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Executive  Summary 


This  report  presents  research  work  performed  under  the  Office  of  Naval  Research  grant  no.  N00014- 
10-1-0455  over  the  four-  year  funding  period  of  2010-2013.  It  includes  contributions  from  Prof. 
Moeness  Amin  (PI),  Dr.  Fauzia  Ahmad  (Research  Professor),  Dr.  Yimin  Zhang  (Research  Professor), 
Dr.  Batu  Chalise  (Postdoctoral  fellow),  Mr.  Pawan  Setlur  (PhD  Student)  and  Mr.  Christian  Debes 
(Visiting  Scholar). 

This  report  deals  with  the  general  problem  of  urban  radar.  It  addresses  major  challenges  in  radar 
imaging  of  simple  and  complex  scenes  inside  enclosed  structures  with  efficient  data  collection 
systems.  The  report  is  divided  into  two  major  sections.  Each  Section  has  its  own  abstract  and 
introduction,  and  also  has  its  own  figure  numbers,  and  references.  The  first  section  includes  a 
comprehensive  discussion  of  key  results  obtained  for  Through-the-Wall  Radar  Imaging  (TWRI)  in 
terms  of  target  detection,  exterior  and  interior  wall  mitigation,  and  multipath  exploitation.  The 
second  section  deals  with  fast  acquisition  and  imaging  using  significantly  reduced  data 
measurements. 

It  has  become  very  evident  that  TWRI  is  emerging  as  a  viable  technology  for  providing  high  quality 
imagery  of  enclosed  structures.  TWRI  makes  use  of  electromagnetic  waves  below  the  S-band  to 
penetrate  through  building  wall  materials.  The  indoor  scene  can  be  illuminated  from  each  antenna, 
and  be  reconstructed  using  the  data  from  the  receive  antennas.  Due  to  the  "see"-through  ability, 
TWRI  has  attracted  much  attention  in  the  last  decade  with  a  variety  of  important  civilian  and  military 
applications.  For  instance,  this  technology  is  employed  for  surveillance  and  detection  of  humans  and 
interior  objects  in  urban  environments,  and  for  search  and  rescue  operations  in  military  situations. 
In  the  first  section,  we  cover  signal  processing  algorithms  that  proved  valuable  in  allowing  proper 
imaging  and  image  recovery  in  the  presence  of  high  clutter,  caused  by  front  walls,  and  multipath, 
caused  by  reflections  from  internal  walls.  Specifically,  we  focus  on  ground-based  imaging  systems 
which  are  relatively  mature,  and  discuss  wall  mitigation  techniques,  multipath  exploitation  methods, 
change  detection  of  moving  target,  and  compressive  sensing  for  fast  data  acquisition. 

Recently,  research  efforts  have  shifted  towards  effective  and  reliable  imaging  under  constraints  on 
aperture  size,  frequency,  and  acquisition  time.  In  this  respect,  scene  reconstructions  are  being 
pursued  with  reduced  data  volume  and  within  the  emerging  compressive  sensing  (CS)  framework.  In 
the  second  section,  we  present  CS  based  scene  reconstruction  techniques  that  address  the  unique 
challenges  associated  with  fast  and  efficient  imaging  in  urban  operations.  Specifically,  we  focus  on 
ground-based  imaging  systems  for  indoor  targets.  We  discuss  CS  based  wall  mitigation,  multipath 
exploitation,  and  change  detection  for  imaging  of  stationary  and  moving  targets  inside  enclosed 
structures. 
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Section  1  -  New  Techniques  for  Improved  Through-the-Wall  Radar 
Imaging:  Theory  and  Applications 


1.1  Introduction 

The  field  of  remote  sensing  has  developed  a  range  of  interesting  imaging  approaches  for  a  variety  of 
applications.  Through-the-wall  sensing  is  a  relatively  new  area  that  addresses  the  desire  to  see  inside 
buildings  for  various  purposes,  including  determining  the  room  layouts,  discerning  the  building 
intent  and  nature  of  activities,  locating  the  occupants,  and  even  identifying  and  classifying  objects 
within  the  building.  Through-the-wall  sensing  is  highly  desired  by  police,  fire  and  rescue,  emergency 
relief  workers,  and  military  operations.  Accurate  sensing  and  imaging  can  allow  a  police  force  to 
obtain  an  accurate  description  of  a  building  in  a  hostage  crisis,  or  allow  firefighters  to  locate  people 
trapped  inside  a  burning  structure.  In  essence,  the  goals  of  through-the-wall  sensing  technology  are 
to  provide  vision  into  otherwise  obscured  areas. 

Each  remote  sensing  application  area  has  driven  different  sensing  modalities  and  imaging  algorithm 
development  based  upon  propagation  characteristics,  sensor  positioning,  and  safety  issues. 
Traditional  optical,  radar,  and  sonar  image  processing  all  begin  with  basic  wave  physics  equations  to 
provide  focusing  to  individual  points.  In  many  radar  applications,  for  example,  data  sampled  from 
many  sensors  are  mathematically  integrated  to  provide  equivalent  focusing  using  free-space 
propagation  assumptions.  Free-space  imaging  is  commonly  seen  in  synthetic  aperture  radar  (SAR) 
techniques  since  atmospheric  distortions  are  often  negligible  and  can  be  safely  ignored  in  first-order 
calculations  [23]. 

Conventional  imaging  approaches  exploit  the  wave  equation  to  compute  the  expected  phase  at  each 
point  in  space  and  time  over  which  the  data  are  collected.  The  complex  phase  front  is  similar  to  the 
spatial  representation  of  the  wavefront  captured  by  a  hologram  [38].  The  complex  returns  can  then 
be  compared  against  the  predicted  returns  from  points  in  the  imaging  target  space  to  focus  on  each 
point  in  that  space.  The  focusing  is  analogous  to  image  reconstruction  in  holography,  where  a  spatial 
pattern  is  projected  back  into  the  originating  target  image  space.  In  true  free-space  conditions,  this 
focusing  approach  represents  a  mathematically  accurate  way  to  perform  imaging.  More 
sophisticated  approaches  extend  beyond  free-space  assumptions  to  allow  for  more  complicated 
propagation  effects,  such  as  adaptive  optics  [78],  atmospheric  correction  for  radar  [23],  and 
matched  field  processing  for  sonar  [79].  Correction  approaches  range  from  simple  wavefront 
calibration  to  more  sophisticated  volumetric  propagation  corrections. 

Free-space  propagation  does  not  apply  to  several  applications  where  transmission  through 
scattering  media  is  encountered,  including  many  modern  imaging  approaches  such  as  geophysical 
sensing,  medical  imaging,  and  more  recently  through-the-wall  imaging  and  sensing.  In  these 
applications,  propagating  signals  diffract  through  a  volume.  As  examples,  geophysical  imaging 
techniques  generally  measure  seismic  propagation  through  the  earth  to  look  for  discontinuities  that 
are  often  indicators  of  oil,  gas,  water,  or  mineral  deposits  [59].  In  medical  imaging,  ultrasound 
tomographic  approaches  account  for  propagation  through  different  tissue  classes  [42]. 

Non-freespace  scattering  applications  are  more  representative  of  the  through-the-wall  sensing 
problem,  albeit  each  has  its  own  distinct  challenges  and  approaches.  In  geophysical  and  medical 
applications,  the  propagation  medium  is  discontinuous  but  still  respectively  fills  the  sensing  volume 
of  the  earth  and  tissue.  To  better  propagate  into  the  volume,  sensors  are  placed  in  direct  contact 
with  the  medium  (e.g.,  seismometers  for  geophysical  sensing,  ultrasound  transducers  for  medical 
imaging).  In  through-the-wall  sensing,  there  are  many  air-material  interfaces  that  dramatically 
change  the  wavefront.  Through-the-wall  sensors  may  be  placed  against  the  front  walls  or  located 
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some  distance  away  from  the  structure.  In  either  case,  attenuation  is  largely  seen  only  in  the 
building  materials  and  contents  rather  than  in  the  large  volumes  of  air  that  occupy  most  of  the  space 
outside  or  within  a  building.  The  rich  through-  the-wall  scattering  environment  makes  volumetric 
tomographic  imaging  approaches  most  relevant  for  through-the-wall  sensing.  Rather  than  using 
free-space  focusing  assumptions,  correcting  for  propagation  effects  may  greatly  improve  the 
imaging  solution  [1-3]. 

Through-the-wall  sensing  is  best  motivated  by  examining  the  applications  primarily  driving  its 
development.  Through-the-wall  sensing  grew  from  the  application  of  ground-penetrating  radar 
systems  to  walls,  with  specific  applications  documented  in  the  literature  since  the  late  1990s 
showing  abilities  to  sense  beyond  a  single  wall  from  near  range  [7,  8,  13,  15  ,17,  31-34,  36,  40,  54, 
67,  68,  70].  Applications  can  be  divided  based  upon  whether  information  is  sought  on  motions 
within  a  structure  or  on  imaging  the  structure  and  its  stationary  contents. 

The  need  to  detect  motion  is  highly  desired  in  situations  like  firefighters  finding  a  child  in  a  burning 
building  or  law  enforcement  officers  locating  hostages  and  their  captors.  Such  applications  can 
resort  to  Doppler  discrimination  of  movement  from  background  clutter  [35,  62,  64,  75].  Motion 
detection  and  localization  can  be  decomposed  into  zero  (0-D),  one  (1-D),  two  (2-D),  or  three- 
dimensional  (3-D)  systems.  A  0-D  system  is  simply  a  motion  detector  and  will  report  any  motion  in 
the  scene.  Such  systems  may  be  useful  to  detect  whether  a  room  is  occupied  or  not,  often  useful  in 
monitoring  applications  such  as  security  systems  or  intruder  detection.  Interior  motion  may  be 
ample  information  for  a  firefighter  to  decide  whether  to  enter  a  room  or  building.  Since  range  or 
angle  is  not  required,  0-D  systems  may  use  continuous  wave  (CW)  tones  as  waveforms.  While  such 
systems  have  exquisite  sensitivity,  it  is  difficult  to  confine  the  sensitivity  to  desired  areas,  so  care 
must  be  taken  to  prevent  undesired  detections  beyond  the  intended  range  or  even  from  the  system 
operator  himself  due  to  multipath  reflections.  0-D  systems  are  not  particularly  useful  in  cases  where 
other  moving  individuals  may  be  within  the  sensor  field  of  view. 

One-dimensional  systems  provide  a  range  to  a  target  but  not  an  angle.  The  extra  dimension  provides 
the  ability  to  separate  and  possibly  discriminate  multiple  targets.  The  range  information  can  help 
bound  whether  a  detection  is  in  the  adjacent  room,  or  perhaps  deeper  within  the  building.  The 
systems  can  obviously  gate  out  detections  from  ranges  beyond  the  ranges  of  interest,  and  may  more 
easily  discriminate  desired  target  motion  from  the  motion  of  the  operator.  The  systems  often  consist 
of  a  single  receiving  antenna  in  addition  to  a  transmitting  antenna  (although  sometimes  a  single 
antenna  can  be  used  for  both  functions).  Motion  of  the  operator  can  be  gated  out  in  range,  or  a 
reference  antenna  pointed  towards  the  radar  operator  can  be  used  for  cancellation.  Ranging  is 
historically  one  of  the  principle  reasons  for  the  development  of  radar  technology,  but  the  short 
stand-off  ranges  of  most  through-the-wall  systems  provide  some  particular  challenges.  Three 
methods  are  often  employed  to  obtain  target  range,  namely,  ultrawideband  pulses  [47,  48,  51,  54, 
80,  83],  dual  frequencies  [9,  10,  62,  91]  or  stepped  CW  [7,  36]. 

Two-dimensional  systems  provide  a  slice  through  the  scene  in  the  range  and  angle  dimensions  [7, 
54,  67].  The  2-D  representation  provides  better  localization  of  the  mover,  at  the  expense  of  a  larger 
antenna  array  whose  elements  are  collinear.  2-D  systems  will  be  subject  to  layover  of  objects, 
meaning  that  objects  out  of  plane  will  appear  to  be  rotated  to  the  imaging  plane,  which  can  lead  to 
distortion  of  the  field  of  view.  3-D  systems  attempt  to  represent  a  volumetric  representation  of  the 
field  of  view  in  the  range,  azimuth,  and  elevation  dimensions  [11].  The  third  dimension  avoids  the 
layover  issue  of  targets  being  projected  into  a  2-D  plane;  at  the  expense  of  a  2-D  array  and  higher 
processing  requirements.  Potentially,  the  third  dimension  can  provide  additional  localization  for 
identification  of  targets.  Height  information  may  allow  discrimination  between  people  and  animals 
such  as  household  pets,  since  radar  cross-section  alone  can  be  unreliable  in  the  presence  of  through- 
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the-wall  multipath  and  other  issues.  The  ultimate  goal  is  higher  resolution  for  even  better  moving 
target  classification.  Both  2-D  and  3-D  processing  techniques  have  generally  used  either 
multilateration  [1,  32]  or  SAR  techniques  [5,  27]  using  either  ultrawideband  pulses  or  stepped  CW 
radars  to  localize  features  in  the  scene. 

Imaging  of  structural  features  and  contents  of  buildings  requires  at  least  2-D  (and  preferably  3-D) 
systems.  It  cannot  rely  on  Doppler  processing  for  separation  of  desired  features,  so  multilateration 
or  SAR  approaches  have  been  the  most  commonly  used  approaches.  The  general  idea  behind 
multilateration  is  to  correlate  range  measurements  from  multiple  sensors  to  specific  points  in  the 
image.  With  sufficient  spatial  diversity  from  a  large  set  of  transmit/receive  combinations,  specific 
reflection  points  will  start  to  integrate  above  the  background  interference.  However,  ambiguities  will 
arise  as  the  number  of  reflection  points  increases  which  can  provide  an  over-determined  system 
relative  to  the  transmit/receive  signal  pairs  which  can  detract  from  image  quality.  SAR-based 
systems  can  be  thought  of  as  a  coherent  extension  of  the  multilateration  concept.  Instead  of 
incoherent  combinations  of  range  returns  from  multiple  transmit/receive  pairs,  coherent  algorithms 
are  used  to  provide  a  complex  matched  filter  to  specific  points  in  the  target  space.  SAR  algorithms 
are  well  established  for  stand-off  imaging  applications  [23].  Stand-off  applications  generally  assume 
free-space  propagation  to  each  point  in  the  target  scene,  although  platform  motion  compensation 
and  atmospheric  effects  are  removed  with  autofocusing  algorithms. 

Most  of  the  SAR  and  multilateration  techniques  usually  neglect  propagation  distortions  such  as 
those  encountered  by  signals  passing  through  walls  and  objects.  Distortions  degrade  the 
performance  and  can  lead  to  ambiguities  in  target  and  wall  localizations.  Free-space  assumptions  no 
longer  apply  after  the  electromagnetic  waves  propagate  through  the  first  wall.  As  a  result,  free- 
space  approximations  may  carry  imaging  systems  through  to  the  first  wall,  but  propagation  effects 
will  then  affect  further  imaging  results.  Shadowing,  attenuation,  multipath,  reflection,  refraction, 
diffraction,  and  dispersion,  all  play  a  role  in  how  the  signals  will  propagate  after  the  first  interface. 
Without  factoring  in  propagation  effects,  imaging  of  contents  within  buildings  will  be  severely 
impacted.  As  such,  image  formation  methods,  array  processing  techniques,  target  detection,  image 
sharpening,  and  clutter  and  multipath  identification  and  rejection  paradigms  must  work  in  concert 
and  be  reexamined  in  view  of  the  nature  and  specificities  of  the  underlying  sensing  problem. 

A  common  occurrence  of  incorrect  localization  is  objects  outside  the  building  being  illuminated  by 
the  reflection  off  the  first  wall  and  subsequently  creating  an  ambiguous  image  visible  inside  the 
building.  Moreover,  the  strong  front  wall  reflection  causes  nearby  indoor  weak  targets  to  go 
undetected.  Multipath  propagation  introduces  ghosts  or  false  targets  in  the  image.  Uncompensated 
refraction  through  walls  can  lead  to  localization  or  focusing  errors,  leading  to  offsets  and  blurring  of 
imaged  targets  [5,  6,  24,  86],  Bragg  scattering  off  repeating  structural  elements,  such  as  rebar  in 
concrete  walls  or  repetitive  voids  in  concrete  block  walls,  can  cause  image  ambiguities  and 
modulation  of  subsequent  wavefronts.  Some  of  the  wall  propagation  effects  can  be  compensated 
for  using  image-focusing  techniques  incorporating  proper  wavefront  corrections  [5,  24].  SAR 
techniques  and  tomographic  algorithms,  specifically  tailored  for  through-the-wall  imaging,  are 
capable  of  making  some  of  the  adjustments  for  wave  propagation  through  solid  materials.  While 
such  approaches  are  well  suited  for  shadowing,  attenuation,  and  refraction  effects,  they  do  not 
account  for  multipath,  Bragg  scattering,  as  well  as  strong  reflections  from  the  front  wall. 

In  this  section,  we  consider  new  algorithmic  that  addresses  some  of  the  aforementioned  unique 
challenges  in  through-the-wall  sensing  operations.  More  specifically,  Section  1.2  deals  with 
techniques  for  the  mitigation  of  the  wall  clutter.  Front  wall  reflections  are  often  stronger  than  target 
reflections,  and  they  tend  to  persist  over  a  long  duration  of  time.  Therefore,  weak  and  close  by 
targets  behind  walls  become  obscured  and  invisible  in  the  image.  Approaches  based  on  both 
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electromagnetic  modeling  and  signal  processing  are  advocated  to  significantly  mitigate  the  front  wall 
clutter.  Section  1.3  presents  an  approach  to  exploit  the  rich  indoor  multipath  environment  for 
improved  target  detection.  It  uses  a  ray  tracing  model  to  implement  a  multipath  exploitation 
algorithm,  which  maps  each  multipath  ghost  to  its  corresponding  true  target  location.  In  doing  so, 
the  algorithm  improves  the  radar  system  performance  by  aiding  in  ameliorating  the  false  positives  in 
the  original  SAR  image  as  well  as  increasing  the  SNR  at  the  target  locations,  culminating  in  enhanced 
behind  the  wall  target  detection  and  localization.  Section  1.4  discusses  a  change  detection  approach 
to  moving  target  indication  for  through-the-wall  applications.  Change  detection  is  used  to  mitigate 
the  heavy  clutter  that  is  caused  by  strong  reflections  from  exterior  and  interior  walls.  Both  coherent 
and  noncoherent  change  detection  techniques  are  examined  and  their  performance  is  compared 
using  real  data  collected  in  a  semi-controlled  laboratory  environment. 

It  is  noted  that  the  section  provides  a  signal  processing  perspective  to  through-the-wall  radar 
imaging.  Outside  the  scope  of  this  section,  there  has  been  significant  work  performed  for  TWRI 
applications  in  the  areas  of  antenna  and  system  design  [16,  44,  65,  83],  wall  attenuation  and 
dispersion  [37,  53,  74],  electromagnetic  modeling  [19,  21,  29,  76],  inverse  scattering  approaches  [26, 
68,  71],  and  polarization  exploitation  [30,  84, 89]. 

The  progress  reported  in  this  section  is  substantial  and  noteworthy.  However,  many  challenging 
scenarios  and  situations  remain  unresolved  using  the  current  techniques  and,  as  such,  further 
research  and  development  is  required.  However,  with  the  advent  of  technology  that  brings  about 
better  hardware  and  improved  system  architectures,  opportunities  for  handling  more  complex 
building  scenarios  will  definitely  increase. 

1.2  Wall  Clutter  Mitigation 

Scattering  from  the  exterior  front  wall  is  typically  stronger  than  that  from  targets  of  interest,  such  as 
a  human  or  a  rifle,  which  have  relatively  small  radar  cross  sections  (RCS).  This  makes  imaging  of 
stationary  targets  behind  walls  difficult  and  challenging.  The  problem  is  further  compounded  when 
targets  are  in  close  vicinity  of  walls  of  either  a  high  dielectric  constant  or  with  layered  structures.  In 
particular,  hollow  cinder  block  walls  contain  an  air-gap  void  within  the  cinder  block  with  disparate 
dielectric  constants.  This  establishes  a  periodic  structure  resonance  cavity  that  traps  electromagnetic 
(EM)  modes.  The  consequence  of  this  layered  composite  structure  on  radar  target  imaging  is  to 
introduce  long  time  constant  relaxations  on  target  detections  in  radar  range  profiles.  Therefore, 
weak  and  close  by  targets  behind  walls  become  obscured  and  sometimes  totally  invisible  in  the 
image.  Thus,  wall  reflections  should  be  suppressed,  or  significantly  mitigated,  prior  to  applying 
image  formation  methods.  One  of  the  simple  but  effective  methods  is  based  on  background 
subtraction.  If  the  received  signals  can  be  approximated  as  the  superposition  of  the  wall  and  the 
target  reflections,  then  subtracting  the  raw  complex  data  without  target  (empty  scene)  from  that 
with  the  target  would  remove  the  wall  contributions  and  eliminate  its  potential  overwhelming 
signature  in  the  image.  Availability  of  the  empty  scene,  however,  is  not  possible  in  many 
applications,  and  one  must  resort  to  other  means  to  deal  with  strong  and  persistent  wall  reflections. 

In  the  past  few  years,  a  number  of  approaches  have  been  proposed  to  mitigate  the  front  wall 
returns  and  effectively  increase  the  signal-to-clutter  ratio.  For  example,  the  wall  reflections  can  be 
gated  out,  the  corresponding  parameters  can  be  estimated,  and  then  used  to  model  and  subtract 
the  wall  contributions  from  the  received  data  [24].  This  wall-dependent  technique  is  effective,  but  its 
performance  is  subject  to  wall  estimation  and  modeling  errors.  Another  method  is  proposed  in  [27], 
which  employs  three  parallel  antenna  arrays  at  different  heights.  The  difference  between  the 
received  signals  at  two  different  arrays  is  used  for  imaging.  Due  to  the  receiver  symmetry  with 
respect  to  the  transmitter,  a  simple  subtraction  of  the  radar  returns  leads  to  wall  reflection 
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attenuation  and  image  improvement.  However,  this  technique  requires  two  additional  arrays  and 
the  effect  of  the  subtraction  operation  on  the  target  reflection  is  unknown  and  cannot  be  controlled. 
In  [29],  the  walls  reflections  are  eliminated  by  operating  the  imaging  radar  in  cross-polarization  as 
the  cross-polarization  returns  from  a  planar  interface  such  as  the  wall  surface  are  theoretically  zero, 
unlike  the  returns  from  humans.  However,  the  cross-polarization  returns  from  the  targets  behind 
walls  are  often  very  weak  compared  to  their  co-polarization  counterparts,  and  as  such,  the  radar 
performance  may  be  limited  by  noise. 

An  eigen-structure  technique  is  developed  in  [77]  to  decompose  the  received  radar  signal  into  three 
subspaces:  clutter  (including  the  front  wall),  target,  and  noise.  Singular  value  decomposition  (SVD  is 
applied  on  the  data  matrix  to  extract  the  target  signatures.  SVD  has  also  been  used  as  a  wall  clutter 
reduction  method  for  TWRI  in  [20,  81].  Another  approach  for  wall  clutter  mitigation  is  based  on 
spatial  filtering.  It  utilizes  the  strong  similarity  between  wall  EM  responses,  as  viewed  by  different 
antennas  along  the  entire  physical  or  synthesized  array  aperture,  is  proposed  in  [87].  The  idea  is  that 
under  monostatic  radar  operation,  the  wall  returns  approximately  assume  equal  values  and  identical 
signal  characteristics  across  the  array  elements,  provided  that  the  extent  of  the  wall  is  much  greater 
than  the  antenna  beamwidth.  On  the  other  hand,  returns  from  targets  with  limited  spatial  extent, 
such  as  humans,  vary  from  sensor  to  sensor.  The  dc  component  corresponding  to  the  constant-type 
return,  which  is  typical  of  walls,  can  be  significantly  reduced  while  preserving  the  target  returns  by 
applying  an  appropriate  spatial  filter  along  the  array  axis.  However,  care  must  be  exercised  in  the 
choice  of  the  spatial  filter  so  that  its  characteristics  cause  minimum  distortion  to  the  target  returns. 
The  spatial-domain  preprocessing  scheme  is  analogous,  in  its  objective,  to  the  moving  target 
indication  (MTI)  clutter  filter  operation  in  the  time  domain. 

In  this  section,  we  present  in  details  two  of  the  aforementioned  methods,  namely,  the  wall 
parameter  estimation  method  and  the  spatial  filtering  technique.  For  the  spatial  filtering  technique, 
we  consider  two  different  types  of  filters,  namely,  moving  average  subtraction  filter  and  the  infinite 
impulse  response  notch  filter,  and  compare  their  effects  on  the  target  image. 

1.2.1  Spatial  Filtering 

1.2. 1.1  Characteristic  of  the  Wall  Reflection 

Assume  a  synthetic  aperture  radar  (SAR)  is  used  in  which  a  single  antenna  at  one  location  transmits 
and  receives  the  radar  signal,  then  moves  to  the  next  location,  and  repeats  the  same  operation 
along  the  axis  parallel  to  the  wall.  Assume  N  antenna  locations.  With  the  wall  reflection,  the  received 
signal  at  the  nth  antenna  due  to  a  single  point  target  of  complex  reflectivity  a  is  given  by, 

z(n,t)  =  sw(t-Tw)  +  os(t-T„)  (1.1) 

where  sw(t )  is  the  signal  reflected  from  the  wall,rw  is  the  two-way  traveling  time  of  the  signal  from 
the  nth  sensor  to  the  wall,  s(t)  is  the  transmitted  signal  convolved  with  the  two-way  transmission 
transfer  function  of  the  wall,  and  rn is  the  two-way  traveling  time  between  the  nth  antenna  and  the 
target.  Note  that  a  is  assumed  to  be  independent  of  frequency  and  tw  does  not  vary  with  the 
sensor  location  since  the  wall  is  parallel  to  the  array.  Furthermore,  if  the  wall  is  homogeneous  and 
much  larger  than  the  beamwidth  of  the  antenna,  the  wall  reflection  sw(t)  will  be  the  same  for  all 
antenna  locations.  This  implies  that  the  first  term  in  (1.1)  assumes  the  same  value  across  the  array 
aperture.  When  the  wall  is  not  parallel  to  the  array  axis,  zw should  be  calculated  for  each  sensor 
location.  Compensations  for  signal  attenuations  due  to  differences  in  the  two-way  traveling 
distances  for  different  antenna  positions  should  be  performed  prior  to  filtering. 
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Unlike  zw,  the  time  delay  zn  in  (1.1)  changes  with  each  antenna  location,  since  the  signal  path  from 
the  antenna  to  the  target  is  different  For  t  =  t0  and  <7  ~  1 ,  the  received  signal  is 


z(n,t0)  =  sw(t0-Tw)  +  s(t0-T„)  (1.2) 

for  n  =  0,1,...9N -1 .  Since  the  time  t  is  fixed,  the  received  signal  is  now  a  function  of  n  via  the 
variable  zn.  We  can  rewrite  (1.2)  as  a  discrete  function  of  n  such  that 

Z/O00  =  vr0  +«/0M  (1-3) 

where  v/q  ~sw(t0  -zw)  and  uto[n\  =  s(t0 -zn).  Since  the  delay  zn  is  not  linear  \nn,utc[n]  for 
w  =  0,l,...,N-l  is  a  nonuniformly  sampled  version  of  s(t) .  Fig.  1(a)  shows  the  signals  received  by 
the  N  antennas  at  a  given  time  instant.  For  example,  suppose  that  the  nc th  sensor  is  the  closest  to 
the  target  (see  Fig.  1(b)).  If  the  wall  diffraction  is  negligible,  then 

Rnc-m  =  y/xf  +  (Sd  +  md )2  =  ^j(xf  +Sd2)  +  2mdSd  +  (md)2  =  ^ R £  +  2 mdSd  +  ( md )2  (1.4) 


where  m  =\n-nc\.\n  most  cases,  the  target's  range  Rn  is  much  larger  than  the  antenna  spacing 
d  .  Then,  R*  »  2 mdSd  +  (md)2 .  Using  Taylor  series  expansion,  we  can  approximate  (1.4)  as, 


R  « R 

nc~m  nr 


1  2  mdSd  +  (md) 

1+ - r- ■ - - 

2  rI 


(1.5) 


Therefore, 


R 


1  2mdSd^jmdf 

2  R 

nc 


(1.6) 


The  received  signal  at  t  =  zn  is  given  by 


z(n,T„c)  =  Sw(T„c  -Tw)+s{T„c  -Tn)*VTn  +5 


ImdSd  +  (md)2 
R  c 


(1.7) 


The  spatial  frequency  transform  of  (1.7)  is  given  by 


z(*»**)=  Ez(«>7»c)e  JKn  *vr  -S(k)  +  S  ( K ) 

c  c  nc  nc 

n- 0 


(1.8) 


where  k  denotes  the  spatial  frequency,  and  ST  ( k )  is  the  target  spatial  signature  at  t  =  t„  such  as 
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N- 1 

\(*)=  1 5 

n=0 


\ 


ImdSd  +  (nid)2 
R  c 


(1.9) 


The  above  analysis  shows  that  separating  wall  reflections  from  target  reflections  amounts  to 
basically  separating  constant  (zero-frequency  signal)  from  nonconstant  valued  signals  across 
antennas,  which  can  be  performed  by  applying  a  proper  spatial  filter. 

1.2. 1.2  Target  Spatial  Signature 

In  the  previous  section,  it  is  shown  that  the  target  signal  across  the  antennas  at  a  fixed  time  (fixed 
down-range)  is  a  nonuniformly  sampled  version  of  the  target  reflection.  In  order  to  minimize  the 
effect  of  the  wall  removal,  which  is  done  by  spatial  filtering,  on  the  target  signal,  it  is  necessary  to 
find  the  target  spatial  frequency  bandwidth  as  viewed  by  the  array  aperture.  This  is  a  critical  factor 
in  designing  the  proposed  spatial  filters  for  the  purpose  of  wall  reflection  removal.  The  target  spatial 
frequency  bandwidth  depends  on  the  target  downrange,  as  shown  in  (1.9),  and  is  also  inversely 
proportional  to  the  extent  of  the  target  signal  return  along  the  array.  At  any  given  time,  the  latter 
depends  on  the  array  length,  target  down  range,  and  the  target  signal  temporal  duration. 

Considering  a  point  target,  which  reflects  the  radar  signal  without  modification,  the  target  spatial 
signature  can  be  analytically  computed.  The  minimum  bandwidth  occurs  when  the  extent  of  the 
target  return  is  maximum.  Two  different  cases  can  be  identified  depending  on  the  target  location  in 
cross-range  relative  to  the  edge  of  the  array.  Fig.  2(a)  shows  the  target  signal  return  in  the  time- 
space  domain  when  the  target  is  located  within  the  array  aperture,  whereas  Fig.  2(b)  shows  the  case 
in  which  the  target  is  outside  the  array  aperture.  In  these  figures,  L  represents  the  array  length,  r 
is  the  signal  temporal  duration,  and  R  is  the  target  range.  Consider  the  case  of  Fig.  2(a),  where  /0  is 
the  cross-range  of  the  target  location  from  the  edge  of  the  array.  Without  loss  of  generality,  assume 
that  l0  is  larger  than  or  equal  to  L  /  2  .  The  maximum  width  of  the  target  return  across  the  array  for 
the  case  of  Fig.  2(a)  is  (Yoon  and  Amin,  2009) 

//  «  ^ 2RSr  +  m\n{^2 RSr  ,L  -/0}  (1.10) 


and  for  the  case  of  Fig.  2(b), 

M  =  TTl  +  2 W8r  -  q2  (1.11) 

In  the  above  equations,  5r  is  the  radar  down-range  resolution  and 

W=TTrTq\  (1.12) 

The  maximum  width  is  achieved  when  the  target  is  located  along  the  centerline  of  the  array 

/W  =min{2A/2  RS^,L}  (1.13) 

Define  /u  /  V2  as  the  effective  width.  Then,  the  maximum  effective  width  is 


/'max.eff  =  min{2  ,  L} 


(1.14) 
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and  the  corresponding  minimum  frequency  bandwidth  is  given  by 


-  (115) 

/^max,eff 

Note  that  kb  mindecreases  as  R  increases.  This  implies  that  the  spatial  filter  should  have  a  narrower 
stopband  to  process  the  more  distant  target. 

1.2.13  Moving  Average  Wall  Removal 

Moving  Average  (MA)  filter  background  subtraction  is  a  noncausal  Finite  impulse  response  (FIR) 
filter,  which  notches  out  the  zero  spatial  frequency  component.  The  MA  subtraction  method  has 
been  effectively  used  for  removal  of  the  ground  reflections  in  ground  penetrating  radar  [18].  The 
effect  of  the  MA  on  the  target  image  can  be  viewed  in  the  context  of  the  corresponding  change  in 
the  point  spread  function  (PSF).  With  the  transmitted  signal  chosen  as  a  stepped-frequency 
waveform  consisting  of  Q  narrowband  signals,  let  z(n,fq)  be  the  signal  received  at  the  «th  antenna 

at  frequency  fq .  Let  z(n,fq)  be  the  signal  after  MA  subtraction.  That  is 


z(n,fq )  =  z{n,fq)~  z(Jq) 


(1.16) 


where  z(fq)  is  the  MAof  z(n,fq)  such  as 


i  m=n+D 

z(fq)  =  ——  I  *(«,/,)  (1-17) 

ZJJ  +  1  m=n-D 

where  2D+1  is  the  filter  length.  When  the  sum  is  taken  over  the  entire  length  of  the  array,  the 
moving  average  becomes  just  an  average  over  the  entire  aperture  and  2D+1=A/.  The  kith  pixel  value 
of  the  image  obtained  by  applying  delay-and-sum  (DS)  beamforming  to  the  received  data  z(n,fq)  is 


given  by 


AM0-1 


b(k,l)=  X  Hz(n,fq)exp{j27fqTn{kJ)} 

n=0q=0 


(1.18) 


where  xn(<k  ^  is  the  two-way  traveling  time,  through  the  air  and  the  wall,  between  the  nth  antenna 
and  the  kith  pixel  location.  The  new  DS  beamforming  image  after  filtering  is  given  by 


b(k,l)  =  lW4)exp{y2^r„>(W)}  =  b(k,l) -  £  sW9)exp {j2nfqtnlkrf  (1.19) 

n- 0  <7=0  n= 0  <7=0 


The  second  term  in  the  above  equation  is  due  to  the  average  subtraction.  This  additional  term 
depends  on  the  signal  waveform,  target  location,  and  the  antenna  positions,  implying  that  the  effect 
of  average  subtraction  varies  according  to  these  parameters.  The  PSF  is  the  DS  image  of  a  point 
source  as  a  function  of  the  down-range  and  the  cross-range.  Fig.  3(a)  shows  the  image  of  the  second 
term  in  (1.19)  when  a  point  target  is  located  6  m  away  from  the  center  of  a  3m  long  antenna  array. 
In  this  case,  the  maximum  value  is  about  0.13  (The  maximum  value  in  PSF  is  one).  The  signal  is  a 
step-frequency  waveform  covering  the  2-3  GHz  band  with  5  MHz  steps.  It  is  noted  that  the 
magnitude  of  this  image  assumes  a  high  value  over  a  large  cross-range  area,  whereas  the  down- 
range  spread  of  the  function  is  small.  Fig.  3(b)  is  the  PSF  without  spatial  filtering.  The  modified  PSF  of 
the  spatial  filtered  signal,  which  is  the  difference  between  Figs.  3(a)  and  3(b),  is  shown  in  Fig.  3(c). 
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We  observe  that  the  size  of  the  mainlobe  of  the  modified  PSF  remains  unchanged  with  about  0.8  dB 
loss  in  the  maximum  value.  However,  the  sidelobes  are  higher  compared  to  the  original  PSF. 
Although  the  average  subtraction  seems  to  work  well  in  this  simulation,  the  limits  of  using  this 
subtraction  method  when  handling  the  distant  targets  can  be  shown  as  follows. 

The  discrete  Fourier  transform  of  z(fq  )  is 


)=  £z(w,/J  exp{-j2xKn/N}=  ))exp{-j2jnai/ N} 


n=0 


n= 0 


N-\ 


(1.20) 


=  Z(Kifq)-z(fq)YJ  exp { -jljucn / N)  =  Z(k,  f  )  -  Nz (/  )S[k] 


n-0 


Note  that 


Z(K,fq)  | 


a:=0 


N-\ 


Y<nJq)=Nt{fq) 

71=0 


Therefore, 


Z(K,fq)  = 


Z(K,fq) 
0 


/c*0 
k  =  0 


(1.21) 


(1.22) 


As  shown  above,  the  subtraction  removes  the  single  spatial  frequency  component  (/c  =  0)  without 
changing  other  components.  The  actual  spatial  frequency  band  that  is  filtered  out  is  always 
MdN  =  \!L,  where  L  is  the  array  length.  When  //maxeff  is  equal  to  L ,  the  target  bandwidth  is  1/Z  . 


In  this  case,  most  of  the  target  signal  return  will  be  removed  by  the  process  of  the  average 
subtraction.  The  range  at  which  /imaxeff  -L  is 


rl 


(1.23) 


For  example,  for  L  =  2  m  and  8r  =  0.075  m  (2  GHz  bandwidth),  RL  is  equal  to  13.33  m,  implying 
that  when  the  target's  down-range  is  larger  than  13.33  m,  the  MA  subtraction  method  will  eliminate 
most  of  the  target  signature,  which  is  undesirable.  Thus,  the  MA  subtraction  is  simple  but  produces 
spatial  filter  characteristics,  which  may  not  be  applicable  to  all  wall  types  and  target  locations. 

1.2. 1.4  Notch  Filtering 

Zero-phase  spatial  MR  filters,  such  as  notch  filters,  are  good  candidates  for  wall  clutter  mitigation. 
The  HR  filter  provides  a  flexible  design.  It  overcomes  the  problem  of  the  fixed  MA  (average)  filter 
characteristics  that  could  be  unsuitable  for  various  wall  and  target  types  and  locations.  Although  the 
MR  filter  spatial  extent  is  truncated  because  of  the  finite  number  of  antennas,  the  filter  is  still  capable 
of  delivering  good  wall  suppression  performance  without  a  significant  change  in  the  target  signal 
when  the  number  of  antennas  is  moderate.  A  simple  MR  notch  filter  that  rejects  zero  frequency 
component  is  given  by 

1  -  jco 

Hnf(oj)  =  ---■■  (1.24) 

1  -ae  J 
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The  positive  parameter  a  (<  1)  determines  the  width  of  the  filter  notch.  Fig.  4  shows  the  frequency 
response  of  the  filter  for  various  values  of  a  and  that  of  MA  subtraction  with  67  antenna  locations. 
Note  the  difference  between  the  notch  filters  and  MA.  The  notch  filter  does  not  have  ripples  and  can 
change  its  stopband  depending  on  the  parameter  a  .  In  the  underlying  problem,  co  is  the  spatial 
frequency  and  the  parameter  a  provides  a  compromise  between  wall  suppression  and  target  signal 
conservation.  This  flexibility  is  an  advantage  of  the  notch  filtering  over  MA  subtraction  method.  As 
shown  in  the  figure,  as  a  moves  closer  to  one,  the  stopband  becomes  narrower  such  that  most  of 
the  signal  except  dc  will  pass  through  the  filtering. 

The  spatial  spectrum  of  wall  reflections  may  have  a  nonzero  width  due  to  unstable  antenna  path  and 
local  inhomogeneities.  In  this  case,  a  wider  filter  stopband  should  be  applied  to  remove  most  of  the 
wall  reflections  around  dc,  but  it  would  also  reduce  the  target  reflection.  Therefore,  the  spatial  filter 
should  be  adjusted  to  the  environment.  It  is  noted  that  the  frequency  spectrum  of  the  target 
reflection  depends  on  the  range  to  the  target,  the  transmitted  waveform,  and  the  distance  between 
antenna  locations.  These  parameters  determine  the  sampling  point  of  the  received  signal  and, 
subsequently,  the  spatial  frequency  characteristics  of  the  target.  If  the  sampled  target  signal  does 
not  have  a  dc  component,  it  will  not  be  severely  affected  by  the  spatial  filter  and,  as  such,  the  target 
image  after  notch  filtering  will  approximately  remain  unaltered.  It  is  noted  that  the  desire  to  have 
the  target  spatial  spectrum  least  attenuated  by  the  spatial  filtering  could  play  an  important  role  in 
designing  signal  waveforms  for  TWRI. 

The  filter  zero-phase  property,  which  is  required  for  a  focused  image  formation,  can  be  achieved  by 
applying  the  notch  filter  twice  in  the  forward  and  backward  directions  [39].  By  two-way  filtering, 
group  delay  can  be  removed.  One  way  to  avoid  the  loss  of  the  target  signal  returns  by  filtering, 
particularly  if  these  returns  vary  slowly  across  the  array  aperture  antennas,  is  to  increase  the  delay 
between  antenna  outputs.  This  is  true  whether  the  MA,  the  notch  filter,  or  any  other  filter  is 
employed.  The  increased  delay  can  be  achieved  by  down  sampling  of  sensors  without  making  the 
sensor  space  too  large  to  cause  aliasing.  When  the  antenna  positions  are  down  sampled  by  M  ,  there 
will  be  A/  different  sets  with  N/M  antennas.  These  M  sets  are  filtered  independently  and  the 
results  are  used  by  the  delay  and  sum  beamforming.  Fig.  5  shows  the  diagram  of  notch  filter 
processing.  The  downsampled  signal  wmt for  n  =  0,\,...,(N  /  M)-\  and  m  =  0,1,.. - 1 


wm  t[n ]  =  z(nM  +  m  -1,0  (1-25) 

The  effective  filter  length  of  the  HR  filter  should  be  considered  when  downsampling  since  the  length 
of  the  downsampled  signal  is  now  N/M .  A  highly  truncated  HR  will  lead  to  undesired  filter 
characteristics. 

1.2. 1.5  Imaging  Results 

A  through-the-wall  SAR  system  was  set  up  in  the  Radar  Imaging  Lab  at  Villanova  University.  A  line 
array  of  length  1.2446m  with  0.0187m  inter-element  spacing  was  synthesized,  parallel  to  a  0.14m 
thick  solid  concrete  wall,  at  a  standoff  distance  of  1.01m.  The  back  and  the  side  walls  of  the  room 
were  covered  with  RF  absorbing  material  to  reduce  clutter.  A  stepped-frequency  signal  covering  1  to 
3  GHz  band  with  2.75  MHz  frequency  step  was  employed.  A  vertical  dihedral  was  used  as  the  target 
and  was  placed  3.04  m  away  on  other  side  of  the  front  wall.  The  size  of  each  face  of  the  dihedral  is 
0.39  m  by  0.28  m.  The  empty  scene  without  the  dihedral  target  present  was  also  measured  for 
comparison. 

First,  we  examine  the  validity  of  the  wall  EM  return  assumption  that  the  spatial  filtering  approach  is 
based  on.  Fig.  6(a)  shows  the  signal  reflected  from  the  concrete  wall.  It  clearly  shows  the  first  few 
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reflections  from  the  wall  and  as  expected,  the  wall  return  assumes  almost  constant  values  across  the 
antenna  positions.  In  order  to  quantify  the  performance  of  the  images,  we  apply  the  target-to- 
clutter  ratio  (TCR)  which  is  commonly  adopted  in  SAR  image  evaluations  [22,  25].  TCR  is  similar  to 
the  MTI  improvement  factor,  except  that  the  latter  is  for  the  time-domain,  whereas  TCR  is  for  the 
image  domain.  TCR  is  calculated  as 


\ 


TCR  =  201og10 


max  |  b(k,l)\ 
{Kl)eAt 

Y}b{k,l)\ 

( k,!)eAc 


(1.26) 


where  At  is  the  target  area,  Ac  is  the  clutter  area,  and  N c  is  the  number  of  pixels  in  the  clutter 
area.  TCR,  in  essence,  is  the  ratio  between  the  maximum  pixel  value  of  the  target  to  the  average 
pixel  value  in  the  clutter  region.  The  clutter  region  is  the  area  where  no  target  is  present  and  the  wall 
reflection  is  severe.  The  clutter  area  is  manually  selected  in  close  vicinity  to  the  wall  where  wall 
reflections  are  most  pronounced.  The  rectangles  depicted  in  the  figures  represent  the  clutter  area 
and  the  target  area  in  our  example.  Fig.  6(b)  shows  the  image  without  any  preprocessing  and  the 
target  is  masked  by  the  wall  response.  Fig.  6(c)  shows  the  result  of  applying  background  subtraction, 
and  the  target  is  clearly  visible.  Fig.  7(a)  and  (b)  demonstrates  that  the  spatial  filtering  approach  is 
effective  in  reducing  the  wall  reflections  without  significantly  compromising  the  target  image.  Fig. 
7(a)  shows  the  DS  image  after  MA  subtraction,  and  Fig.  7(b)  shows  the  DS  image  after  the  notch 
filtering.  TCRs  in  these  figures  are  2.2  dB  (with  wall),  26.2  dB  (background  subtraction),  13.4  dB 
(notch  filtering),  and  15.3  dB  (MA  subtraction).  The  background  subtracted  image  provides  the 
highest  TCR,  while  MA  subtraction  and  notch  filtering  provide  comparable  performance. 

1.2.2  Wall  Parameter  Estimation,  Modelling,  and  Subtraction 

1.2. 2.1  Approach 

In  this  approach,  the  effect  of  the  wall  on  EM  wave  propagation  is  achieved  through  three  steps.  The 
wall  parameters  are  accurately  estimated.  Then,  the  reflected  signal  from  the  wall  is  properly 
modelled.  Finally,  the  modelled  signal  is  subtracted  from  the  measured  signal. 

Building  walls,  such  as  brick,  adobe,  and  poured  concrete  walls,  can  be  modelled  by  homogeneous 
dielectric  slabs.  Using  the  geometric  optics  approach,  the  reflection  coefficient  y  for  vertically  and 
horizontally  polarized  incident  waves  through  the  wall  is  given  by  [24,  45] 


r  = 


*0  i  + 

]  +  R0lRwej2k"d 


*10  --*oi  _ 


l-P  10 

1  +  PlO 


(1.27) 


where  pl0  =  k0x  !kXxi°r  horizontal  polarization  and  pl0  =  £rk0x  /klx  for  vertical  polarization  of  the 
incident  fields.  In  (27),  sr  is  the  relative  dielectric  constant,  d  is  the  thickness  of  the  wall,  and  k0x 
and  £1;care  the  normal  components  of  the  propagation  constants  in  air  and  in  the  dielectric, 
respectively.  The  expressions  for  the  vertical/transverse  magnetic  and  horizontal/transverse  electric 
unit  vectors  are  given  in  [45]. 

The  wall  parameters  can  be  estimated  from  the  time-domain  backscatter  at  a  given  location  [24,  73]. 
Basically,  the  response  of  the  wall  is  gated  from  the  total  backscatter  signal  and  transferred  to  the 
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frequency  domain.  The  mean-squared  error  between  the  calculated  reflection  coefficient  y  of  (1.27) 
and  the  measured  reflection  coefficient  ym  given  by 

1 

s=-^rZ\rMn)-r(fn)\2  (i.28) 

n= 1 

is  computed,  and  its  minimum  is  searched  for  the  wall  thickness  and  permittivity  [24,88].  In  (1.28), 
Nf  is  the  number  of  frequency  points  and  fn  is  the  nth  frequency. 

The  estimated  wall  parameters  are  then  used  to  compute  the  wall  EM  reflection.  This  could  be 
accomplished  either  numerically  using  EM  modelling  software,  or  in  case  of  far-field  conditions, 
analytically  using  the  following  expression  [24] 

GX  ej2k"x° 

0„)=r - m  =0  ,kn)  (1.29) 

8  n  x0 


where  x0  is  the  standoff  distance  of  the  radar  from  the  wall,  A  is  the  wavelength,  G  is  the  antenna 
gain,  kn=27tfn/c  is  the  «th  wave  constant,  c  is  the  speed  of  light  in  free  space,  and  T(di  =  0,kn)  is 
the  reflectivity  matrix  from  the  wall  at  normal  incidence  angle  and  is  given  by, 


r  = 


rv 

o 


o 

ft. 


(130) 


In  (30),  the  subscripts  V  and  'h',  respectively,  denote  the  vertical  and  horizontal  polarization  and 
Xv/h  >s  given  by  (1.27).  This  estimated  wall  reflection  should  be  coherently  subtracted  from  the  total 
received  signal  in  order  to  isolate  the  signature  of  the  target  on  the  other  side  of  the  wall. 

The  isolated  target  return  can  then  be  processed  to  generate  an  image  of  the  scene.  The  effects  of 
the  wall  transmissivity  on  the  target  image  are  the  following:  1)  dislocation  of  the  target  from  its 
actual  position  in  range  and  2)  significant  degradation  in  cross-range  resolution  [5,  24].  To 
compensate  for  the  effects  of  transmission  through  the  wall,  a  compensation  factor  proportional  to 
the  inverse  of  the  squared  of  the  wall  transmittivity  can  be  used  in  the  conventional  free-space 
image  formation  methods.  In  cases  where  transmissivity  is  low  and  the  signal  is  noisy,  division  by  a 
small  noise-affected  number  may  cause  significant  distortion.  In  such  cases,  one  may  use  only  the 
phase  of  transmission  coefficients  as  a  compensation  factor  [24]. 

1.2. 2.2  Imaging  Results 

An  experimental  study  was  conducted  at  the  Radiation  Laboratory  in  University  of  Michigan.  The 
measurement  setup  consists  of  an  HP8753D  vector  network  analyzer,  two  double  ridge  horn 
antennas  with  operational  frequency  range  of  1  to  18  GHz,  an  XY table  along  with  a  control  unit,  and 
a  personal  computer.  The  vertically  aligned  antennas,  one  for  transmission  and  the  other  for 
reception,  are  mounted  on  a  vertical  wooden  rod  (along  the  z-axis),  which  is  attached  to  the  carriage 
of  the  XY  table.  A  2.37m  x  1.77m  wall  composed  of  poured  concrete  blocks  is  made  on  top  of  a  layer 
of  cinder  blocks  inside  the  laboratory.  The  underlying  cinder  block  layer  is  used  to  line  up  the 
antennas  in  height  at  approximately  the  middle  of  the  concrete  block  wall.  Figure  8(a)  shows  the 
side  view  of  the  measurement  setup.  A  small  trihedral  corner  reflector  with  pentagonal  panel 
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geometry  is  used  as  a  point  target  behind  the  wall,  as  shown  in  Fig.  8(b).  The  back  corner  of  the 
trihedral  (scattering  phase  center)  is  at  x  =  -0.71  m  and  at  the  same  height  as  that  of  the  receiver 
antenna,  i.e.,  1.28  m  above  the  ground  plane.  The  transmitting  antenna  is  attached  about  0.25  m 
below  the  receiving  antenna  on  the  wooden  rod.  The  antennas  are  moved  along  a  scan  line  of  length 
95.88  cm,  parallel  to  the  wall,  with  a  spacing  of  2.04  cm.  The  apertures  of  both  antennas  are  about 
0.45  m  away  from  the  wall.  The  frequency  of  operation  is  from  1  to  2.5  GHz  and  the  frequency  step 
is  12.5  MHz.  It  is  noted  that  the  antennas  are  in  the  far-field  region  of  the  target. 

An  image  is  first  constructed  under  assumption  of  free-space  propagation  by  using  the  total  received 
signal,  and  is  shown  in  Fig.  9(a).  Here,  in  addition  to  the  target,  the  wall  is  also  imaged,  primarily  as 
two  parallel  lines  showing  the  front  and  back  boundaries.  Since  backscatter  from  the  wall  is  very 
strong,  the  sidelobes  generated  by  image  formation  spill  over  into  desired  image  domain,  which  is 
manifested  as  multiple  lines  parallel  to  the  wall  surfaces.  Also,  the  target  appears  blurred  and  its 
location  is  biased.  Figure  9(b)  shows  the  image  after  the  estimate,  model,  and  subtract  approach  was 
applied.  The  wall  thickness  and  the  permittivity  were  estimated  using  (1.28)  as  20  cm  and  5.7  +  J0.6, 
respectively.  The  wall  reflection  was  computed  using  (1.29)  and  subtracted  from  the  total  received 
signal,  and  the  resulting  signal  was  used  for  image  formation.  For  comparison,  the  signal  prior  to  wall 
return  subtraction  was  used  for  imaging  wherein  compensation  for  the  effect  of  transmission 
through  the  wall  was  applied  and  is  shown  in  Fig.  9(c).  Although  the  imaged  target  is  at  the  correct 
location  and  is  refocused  in  Fig.  9(c),  existence  of  the  wall  sidelobes  in  the  image  is  still  evident.  By 
comparing  Figs.  9(b)  and  (c),  a  substantial  improvement  in  clutter  rejection  is  observed  after 
subtraction  of  the  modeled  wall  return. 

1.3  Multipath  Exploitation 

The  existence  of  targets  inside  a  room  or  in  an  enclosed  structure  introduces  multipath  in  the  radar 
return,  which  results  in  false  targets  or  ghosts  in  the  radar  images.  These  ghosts  lie  on  or  near  the 
vicinity  of  the  back  and  side  walls,  leading  to  a  cluttered  image  with  several  false  positives.  Without 
a  reasonable  through-the-wall  multipath  model,  it  becomes  difficult  to  associate  a  ghost  to  a 
particular  target.  Identifying  the  nature  of  the  targets  in  the  image  and  tagging  the  ghosts  with  their 
respective  target  according  to  a  developed  multipath  model,  although  important  to  reduce  false 
alarms,  is  not  the  final  goal  of  a  high-performance  through-the-wall  imaging  system.  Since  each 
multipath  provides  some  information  about  the  target,  it  becomes  prudent  to  utilize  the  multipath 
rather  than  ignoring  it,  once  identified.  The  utilization  or  exploitation  of  the  multipath  to  one's 
advantage  represents  a  paradigm  shift  when  compared  to  the  classical  approach  of  either  ignoring 
or  mitigating  it. 

For  through-the-wall  radar  imaging  applications,  the  existence  of  multipath  has  been  recently 
demonstrated  for  stationary  targets  in  [4,  19,  28].  In  [4],  the  authors  use  distributed  fusion  to 
remove  the  false  targets  caused  either  from  multipath  or  target  interactions  for  stationary  scenes 
after  suitable  image  registration.  In  [9],  time  reversal  techniques  are  used  to  alleviate  ghosts  and 
clutter  from  the  target  scene  and  in  [28],  a  synthetic  aperture  radar  (SAR)  based  image  of  a  human 
inside  a  room  is  shown  along  with  possible  multipath  ghosts.  However,  no  more  rigorous  multipath 
modeling  and  analysis  are  presented  in  the  aforementioned  references  [4,  19,  28]  and  the 
references  therein. 

In  this  section,  we  derive  a  model  for  the  multipath  in  an  enclosed  room  of  four  walls.  The  model 
considers  propagation  through  a  front  wall  and  specular  reflections  at  interior  walls.  A  SAR  system  is 
considered  and  stationary  or  slow  moving  targets  are  assumed.  Although  the  multipath  model 
presented  deals  with  walls,  reflections  from  the  ceiling  and  floor  can  be  readily  handled,  similar  to 
the  work  in  [61].  We  demonstrate  analytically  that  the  multipath  as  seen  by  each  sensor  is  displaced, 
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and,  therefore,  we  derive  the  actual  focusing  positions  of  the  multipath  ghosts  in  downrange  and 
crossrange.  The  multipath  model  permits  an  implementation  of  a  multipath  exploitation  algorithm, 
which  associates  as  well  as  maps  each  target  ghost  back  to  its  corresponding  true  target  location.  In 
so  doing,  the  exploitation  algorithm  improves  the  radar  system  performance  by  ameliorating  the 
false  positives  in  the  original  SAR  image  as  well  as  increasing  the  signal-to-clutter  ratio  (SCR)  at  the 
target  locations,  culminating  in  enhanced  behind  the  wall  target  detection  and  localization.  It  is 
noted  that  the  exploitation  algorithm  only  maps  back  target-wall  interactions;  target-target 
interactions  are  left  untreated.  The  multipath  exploitation  algorithm  is  inspired  by  the  work  in  [43], 
wherein  the  shadowed  regions  of  a  target  are  revealed  via  its  multipath  returns.  The  difference, 
however,  is  that  we  are  not  striving  to  reveal  the  shadowed  regions  of  the  target.  In  our  case,  we 
deal  with  targets  with  arbitrary  dielectric  constants,  and  wish  to  associate  and  map  each  multipath 
to  its  true  target  location. 

1.3.1  Image  Formation  Algorithm 

In  this  section,  we  describe  the  through-the-wall  delay-and-sum  beamforming  approach  [5,  11,  12]. 
We  consider  a  SAR  system  in  which  a  single  antenna  at  one  location  transmits  a  wideband  signal  and 
receives  the  radar  return,  and  then  moves  to  the  next  location  and  repeats  the  same  operation 
along  the  axis  parallel  to  the  front  wall.  Assume  N  monostatic  antenna  locations.  The  setup  is  as 

depicted  in  Fig.  10.  Consider  a  single  point  target,  located  atx,  =[xtiyt]T.  At  each  antenna 
location,  the  radar  transmits  a  pulsed  waveform  s(/),  where  *tf  indexes  the  time  within  the  pulse, 
and  measures  the  reflected  signal.  The  target  return  at  the  «th  antenna  location  is  given  by, 


r„  ( t )  =  r„2  (x, )a(x,  )s(t  -  2  z„  (x, ))  (1.31) 

where  a(xt)  represents  the  target  reflectivity,  Tn  (x,)is  the  complex  amplitude  associated  with  the 
one-way  propagation  through  a  dielectric  wall  for  the  wthantenna  location  [63],  and  rn(\t) 
represents  the  one-way  through-the-wall  propagation  delay  from  the  nth  antenna  location  to  the 
target. 

The  scene  of  interest  comprises  several  pixels  indexed  by  the  downrange  and  the  crossrange.  The 
complex  composite  return  from  the  Ath  pixel  location  x*  is  obtained  by  applying  time  delays  to  the 
N  received  signals,  followed  by  weighting  and  summing  the  results.  That  is, 


N- 1 

r(t,  X*  )  =  Z  wnrn  (f  +  2r„  (x*  )) 


n= 0 


(1.32) 


The  signal  r(t,\k)  is  passed  through  a  matched  filter,  with  impulse  response  h(t)  =  s*(-t).  The 
complex  amplitude  assumed  by  the  pixel  \k  in  the  image  /(•)  is  obtained  by  sampling  the  matched 
filter  output  at  time  f  =  0, 


I(*k  )  =  {r(*>  *k  )  *  W)}  1/3=0  (1-33) 

where  denotes  the  convolution  operation.  Equations  (1.32-1.33)  describe  the  standard 
beamforming  approach  in  through-the-wall  radar  imaging.  It  is  noted  that  if  the  imaged  pixel  is  in 
the  vicinity  of  or  at  the  true  target  location,  then  the  complex  amplitude  in  (1.33)  assumes  a  high 
value  as  given  by  the  system's  point  spread  function.  The  process  described  in  (1.32)  and  (1.33)  is 
carried  out  for  all  pixels  in  the  scene  of  interest  to  generate  the  image  of  the  scene. 


14 


The  beamforming  approach  as  described  above  is  pertinent  to  a  point  target  in  a  two-dimensional 
(2D)  scene.  Extensions  to  three-dimensional  (3D)  scenes  and  spatially  extended  targets  are  discussed 
in  [11,  85].  It  is  noted  that,  in  the  above  beamforming  description,  we  have  not  considered  the 
multipath  returns;  neither  have  we  addressed  the  specificities  of  calculating  the  delays,  Tn(-), which 
will  be  treated  in  detail  in  the  following  sections.  Without  loss  of  generality,  we  now  assume  that 
the  weights  wn  =  1,V«  and  the  target  has  unit  reflectivity. 

1.3.2  Multipath  Model 

Consider  a  room  under  surveillance  using  a  SAR  system.  A  priori  knowledge  of  the  room  layout,  i.e., 
wall  locations  and  material  properties,  is  assumed.  The  scene  being  observed  relative  to  the  nth 
sensor  location  is  as  shown  in  Fig.  11.  The  origin  is  marked  as  point  'O'  in  the  figure,  and  the 
standard  convention  for  the  positive  x-  and  y-  axes  is  assumed.  The  nth  sensor  location  is  given  by 

=[-Dxn,0]J  .The  front  wall  has  a  thickness  ^and  dielectric  constant  ^.For  notational 
simplicity,  the  back  and  the  side  walls  are  also  assumed  to  have  sx  as  the  dielectric  constant.  The 
side  walls,  labeled  as  wall-1  and  wall-3,  have  a  length  Dx,  whereas  the  front  wall  and  the  back  wall 

(wall-2)  have  a  length  Z)2.  The  target  is  stationary  and  at  location  A  =  [-xt,yt]7  .The  standoff 
distance  from  the  front  wall  is  constant  for  each  sensor  location  constituting  the  synthetic  aperture, 
and  is  denoted  as  Dy.  In  the  figure,  we  consider  the  direct  path,  referred  to  as  path-A,  and  three 

additional  paths,  namely,  paths-B,  C,  and  D,  which  correspond  to  single-bounce  multipath.  In 
general,  there  exist  other  paths  which  can  contribute  to  the  multipath  returns;  these  include 
multiple  bounces  from  the  back  and  side  walls  as  well  as  paths  that  include  multiple  reflections 
within  each  of  the  four  walls  themselves.  Examples  of  the  former  and  the  latter  are  provided  in  Fig. 
12.  Such  paths  are  defined  as  higher-order  multipath.  Hence,  in  Fig.  11,  we  have  considered  only 
first-order  multipath.  It  is  noted  that  higher-order  multipath  are  in  general  weaker  compared  to  the 
first-order  multipath  due  to  the  secondary  reflections  and  refractions  at  the  various  air-wall 
interfaces  and  high  attenuation  in  the  wall  material.  Therefore,  we  choose  to  exclude  these  paths 
from  the  model. 

The  walls  are  assumed  smooth  with  specular  reflections.  The  smoothness  assumption  is  valid,  since 
for  through-the-wall  radars,  the  wavelength  of  operation  is  much  larger  than  the  roughness  of  the 
walls.  Specular  reflections  are  a  direct  consequence  of  the  wall  smoothness,  and  necessitate  that  the 
angle  of  incidence  be  equal  to  the  angle  of  reflection.  Note  that,  in  general,  the  back  and  side  walls 
may  each  be  of  a  different  material  (interior  or  exterior  grade),  and  thus  may  assume  different 
values  for  the  wall  thickness  and  dielectric  constant.  If  interior  reflections  inside  these  walls  are 
considered,  then  the  thickness  of  these  walls  is  required  [19,  28,  43].  However,  since  we  have 
ignored  such  higher-order  paths,  we  only  require  the  dielectric  constant  of  the  back  and  side  walls  to 
be  known.  As  the  EM  wave  propagates  through  the  front  wall,  it  bends  at  the  medium  discontinuity 
as  dictated  by  Snell's  law,  i.e.,  each  of  the  paths,  as  seen  in  Fig.  11,  has  an  associated  angle  of 
incidence  and  an  angle  of  refraction.  For  example,  the  angles  of  incidence  and  refraction  for  path-B 

are  denoted  as  \f/$  and  respectively.  Similar  nomenclature  follows  for  the  remaining  paths. 

Let  the  reflection  points  on  wall-1,  2,  and  3  be  denoted  by  Bn,Cn,  and  Dn,  with  respective 
position  vectors  B„  =  [0,yBn]T  ■  C„  =  [~xCn  ,£>i  +  Dyf ,  and  D„  =[-D2,yDnf .  It  is  clear  that 
these  position  vectors  are  dependent  on  the  sensor  location.  The  one-way  path  delays  for  the  four 
paths,  with  the  antenna  at  the  «th  location,  are  denoted  by  p  e{A,  B,C9D},  and  can  be 
derived  from  the  geometry  as 
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TAn)  =  tfnfr  sec(^})+(^,  -rf,)sec(^)))/c 
tb']  =  (d\4^\  sec^rB  )  +  Os»  - ) sec(i//^} ) + X,  cosec(^'Bn) )) / c 
*£°  =  sec(^))-sec(^))  +  (2D1  +2D, -y,)sec(W^))/c 

t:(d  =  (d\ -[ih  sec(^D  )  +  (^d„  “  d\ ) sec(ViD  )  +  (D2-x,) cosec(i/^) )) / c 

with  c  being  the  speed  of  light  and  the  coordinates  of  5„,  C„,  and  £)„  given  by 

J'Bn  =yt~xi  cot (y/ff) 

*<:«  =  *»  +  (A  +Dy~yi)  tan(  F*? ) 

=  JV  -  (A  -  */ )  co^(w\b ) ) 

Equations  (1.34)  and  (1.35)  depend  on  the  angles  of  incidence  and  refraction.  Since  the  standoff 
distance  is  a  known  constant,  by  projecting  the  various  paths  to  the  x-axis,  we  obtain  the  following 
equations,  which  are  useful  in  calculating  the  various  angles. 

d\  tan((//^})  +  0,  -rf,)tan(^f  )-£>,„+*,  =0 
d\  tan (Y{$)  +  ( yBn  ~  d\) tan(^)  -  Dxn  =0 

(2 A  +  2D^  -  )  tan(^}c  )  +  </,  (tan(^^ )  -  tan(yj"  ))  “  An  +  */  =  0 

A  tan(^} )  +  (>>£«  -  )  tan((^  jp  )-D2+Dxn=  0 


(1.34) 


(1.35) 


The  angles  of  refraction  can  be  determined  from  Snell's  law  as, 

sin(v'L”)) 


vty  =sin  1 


,  pe{A,B,C,D} 


(1.37) 


Substituting  the  angles  of  refraction  from  (1.37)  in  (1.36),  and  using  (1.35),  we  obtain  a  set  of 
equations  that  can  be  solved  numerically  for  the  angles  of  incidence  by  using  the  Newton  method. 
The  angles  of  refraction  can  then  be  obtained  using  (1.37). 

We  are  now  in  a  position  to  write  the  radar  signal  return  from  the  single  target  scene  as  a 
superposition  of  the  direct  path  and  the  multipath  returns.  For  the  wth  sensor  location,  we  obtain 

r„(0=  E  T2„s(t-2z<pn))  +  £  TpnTqns(t  -tf  (1.38) 

pe{A,B,C,D)  ( p,q)e{A,B,C,D } 

p*q 


where  Tpn  is  the  complex  amplitude  associated  with  reflection  and  transmission  coefficients  for  the 

one-way  propagation  along  path-p,  and  depends  on  the  angles  of  incidence,  the  angles  of  refraction, 
and  dielectric  properties  of  the  walls  [63].  In  (1.38),  the  first  summation  captures  the  two-way 
returns  along  the  direct  path  and  each  multipath.  More  specifically,  the  signal  propagates  along  a 
particular  path -p,  p  e  { A ,  B ,  C,  D},  reaches  the  target  and  retraces  the  same  path  back  to  the  radar, 
i.e.,  path-p.  The  multipath  returns  due  to  the  combination  paths  are  captured  by  the  second 
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summation,  i.e.,  the  wave  initially  travels  to  the  target  via  path -p  and  returns  to  the  radar  through  a 
different  path-q,  q^{A,B,C,D). 

From  Fig.  11  and  (1.35),  we  observe  that  the  coordinates  of  the  points  of  reflection  at  the  back  and 
side  walls  are  sensor  dependent,  implying  that  the  multipath  presents  itself  at  different  locations  to 
the  different  sensors.  Location  of  the  multipath  as  seen  by  each  sensor  and  the  actual  focusing 
position  of  the  resulting  ghost  in  the  image  are  discussed  next. 

1.3.2. 1  Multipath  Locations 

For  simplicity  of  analysis,  consider  the  scenario  in  Fig.  11,  but  without  the  front  wall.  We  first  focus 
on  the  multipath  originating  from  wall-1,  i.e.,  the  signal  travelling  to  the  target  via  path-A,  and 
following  path-B  back  to  the  radar  or  vice  versa.  The  goal  is  to  find  the  location  of  the  multipath 

corresponding  to  the  above  combination  path  as  seen  by  each  of  the  N  sensors.  Let  the  point  P™1 , 
represented  by  the  vector  [-x"1  ,y^  ]'  ,  denote  the  multipath  location  as  seen  by  the  «th  sensor. 
The  superscript  * stresses  that  the  multipath  is  associated  with  the  reflection  at  wall-1  and  the 

subscript ‘t*  indicates  that  the  multipath  is  due  to  the  target  located  at  [-xnyt  ]7  .  Reflecting  path  B 
about  wall-1  yields  an  alternate  radar-target  geometry,  as  shown  in  Fig.  13.  We  obtain  the  virtual 
target  denoted  by  point  A1  with  coordinates  [xt,yt]i  such  that  the  distance  ABn=AxBn,  which 
implies  that 


AR„  +  ABn  +  BnRn  =  AR„  +  AlR„  (1.39) 

In  other  words,  the  combination  path  comprising  of  path-A  and  path-B  has  the  same  length  as  path- 

A  and  the  path  connecting  the  «th  sensor  to  the  virtual  target  at^1.  Due  to  the  monostatic  nature 
of  the  radar  system,  this  combination  path  appears  to  the  «th  sensor  as  the  traditional  two  way 

path  of  length  2RnP”] .  Therefore,  we  have 


R  P  1 

n 1  n 


AR„  +  a}  R 


(1.40) 


Additionally,  we  obtain  the  virtual  radar  at  Rj,  =  [Dxnfi]7  on  the  other  side  of  the  wall,  as  shown  in 
Fig.  13.  Considering  Fig.  13,  it  is  readily  observed  that  the  first  order  multipath  has  the  same  time 
delay  as  a  bistatic  configuration,  comprising  the  radar  and  the  virtual  radar.  That  is,  in  terms  of 
range,  we  have 


R„A  +  R„Al  =  RnA[  +RXnAx 


(1.41) 


where  RnA]  +  RxnAx  represents  the  bistatic  configuration.  In  the  bistatic  case,  the  constant  range 
contour  corresponding  to  the  multipath  is  an  ellipse,  which  has  foci  at  the  radar  and  virtual  radar 
and  passes  through  the  target  and  the  virtual  target  locations.  On  the  other  hand,  the  monostatic 
measurement  scenario  described  by  (1.40)  has  a  circular  constant  range  contour,  centered  at  the 
radar.  Since  the  two  measurements  are  equivalent,  the  location  of  the  multipath  corresponds  to  the 
point  of  intersection  of  the  bistatic  elliptical  and  the  monostatic  circular  constant  range  contours. 
The  equations  for  the  ellipse  and  circle  are  given  by, 
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(1.42) 


Off)2  ,  Off)2  =1 

(RnP„w')2  (R„P„W')2  ~D2X„ 

(*«'  +Dxn)2  +  Off)2  =(R„Pnw')2 


Solving  (42)  for  an  intersection  point,  we  can  readily  see  that  =0  is  the  only  solution.  In  other 
words,  regardless  of  the  location  of  the  target,  its  first-order  multipath  always  falls  on  the  wall.  The 
y-coordinate  of  P™1  can  then  be  derived  as 


>ff  =±i 


KlM  ,  2 

A  Uxn 


(1-43) 


In  (1.43),  the  positive  y-coordinate  is  the  desired  value  as  the  other  solution  lies  behind  the  radar.  It 
is  clear  from  (1.43)  that  the  multipath  location  is  sensor  dependent.  Therefore,  the  locations  of  the 
target  multipath  corresponding  to  the  various  sensor  positions  will  be  displaced  from  one  another. 
In  the  presence  of  the  front  wall,  the  multipath  remains  on  the  wall,  but  the  equations  for  the  ellipse 
and  circle  are  different  than  that  presented  in  (1.42).  Hence,  for  the  scenario  comprising  of  the 
front-wall  and  wall-1,  the  multipath  appears  on  wall-1,  with  its  y-coordinate  given  by, 


y7n  = 


(  ( n )  ,  (w) 

CT\  '  +  CT  n  ' 


cos<rf’>"') 


(1.44) 


where  and  y/^Wl  are,  respectively,  the  angles  of  incidence  and  refraction  for  the  ghost.  The 

solution  in  (1.44)  depends  on  the  angle  y/^W] ,  which  can  be  obtained  by  solving 


d\  tan(^)vVl )  + 


f  (n)  ( n )  'N 

C-l±^pL-  +  dl^  secU^) 

V  J 


sin(  ^?)W')-DX=0 


y,^'  =sin-‘ 


(1.45) 


So  far,  we  have  considered  the  multipath  corresponding  to  wall-1;  there  exist  two  other  multipath 
returns  tagged  to  the  remaining  walls,  namely,  wall-2,  and  wall-3.  We  can  readily  show,  using  similar 
analysis,  that  the  multipath  associated  with  wall-2  appears  on  wall-2,  i.e.,  y™2  =  Dy  +  D{f  at  an  x- 
coordinate  given  by 


w2 

tn 


=  -sin(^")W2)x 


sec(^”)M’2 )  -  tan(^^,”)W2  )cosec(^}")W2 )) 

(n)w  CtP+CtP 

+  Dx  cosec(y/-,  2 ) - — — 

2  y 


(1.46) 


The  angles  in  the  above  equation  are  determined  by  solving  the  equations, 
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(1.47) 


d\  secityf'0"'2  )  +  (£>,-</,) sec(y/j")W2  )  +  dxJe \  sec(^^")M’2 )  - 


cri")  +  ct 


(«) 


=  0 


ir™*  =  sin-1 


sinfr  <”>»,)' 

Pi 


Likewise,  the  multipath  associated  with  wall-3  appears  at  wall-3,  i.e.,  x”3  =- D2 ,  and  at  a  y- 
coordinate  given  by. 


(n)  (n) 

CZ\  4*  CTn ' 


yln  = 


d]  (sec(y4")M'3 )  -  ^  sec(^,”)M'3 )) 


cos(^")M'3 ) 


(1.48) 


The  respective  angles  can  be  obtained  by  solving, 


f  («)  ( n ) 


+  sec(^3) 


sin(^,v”)W3 )-  D2  +  Dx  =0 


(«)w3  •  -1 

y/)tf  3  =sin 


sin(^  ^  ,”)W3 ) 


(1.49) 


From  the  above  equations,  it  is  again  observed  that  the  N  sensors  view  the  multipath,  resulting  from 
the  combination  paths  associated  with  a  particular  wall,  at  different  locations  on  that  wall.  That  is, 
the  multipath  maybe  regarded  as  a  moving  target.  As  a  result,  when  applying  beamforming,  the 
multipath  ghost  appears  at  a  different  pixel  in  the  vicinity  of  the  true  multipath  locations.  We  note 
that  the  multipath  ghost  will  lie  inside  the  room  except  when  the  target  is  near  the  corners  of  the 
room.  In  this  case,  as  the  multipath  corresponding  to  the  N  sensors  may  appear  along  an 
extrapolation  of  the  wall,  the  multipath  ghost  may  appear  outside  the  room.  We  further  note  that 
the  virtual  target  corresponds  to  two-way  propagation  along  the  single-bounce  multipath,  and  is 
readily  seen  to  lie  outside  the  room  perimeter.  The  multipath  focusing  pixel  analysis  for  the 
combination  paths  is  discussed  next. 

1.3. 2. 2  Multipath  Focusing  Analysis 

Consider  Fig.  14,  which  shows  the  multipath  locations  w.r.t  to  walls-1-3,  and  the  focusing  pixels  w.r.t 
to  these  walls.  Consider  the  multipath  locations  associated  with  wall-1;  further  assume  that  the 
focused  multipath  ghost  appears  at  a  pixel  location  given  by 

X?  =  [-Ax,  +  ,  Ay,  +  ]'  (1.50) 

where  xJJ1  ,y™\Y  is  the  true  multipath  location  corresponding  to  the  first  sensor  position. 

Using  a  first  order  Taylor  series  expansion,  which  is  valid  under  conditions  of  a  small  aperture  [82] 
and  when  the  ghost  is  in  the  vicinity  of  the  true  multipath  locations,  we  obtain  the  difference  in 
propagation  path  length  between  the  multipath  ghost  location  and  the  true  multipath  location  w.r.t 
the  first  sensor  position  as 
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(1.51) 


DxX  -  D  tan(y4')Wl )  -  dx  tan(^J)w' )  -  x*1 
Agx  =  Ax,  . . . . . . . 

yj(Dx ,  -  Dy  tan^"* )  -  rf,  tan(^^)w' )  -  *,?  )2  +  (y*  -Dy-dx  )2 

y%-D,-d\ 

+  Ay]  ===== 

V(^1  -  Dy  tan(v41)W|  )-^i  tan(y/^l)w‘ ) - xp  )2  +  (j#  -  Dy  -dx)2 

Following  the  analysis  in  [82],  (1.51)  can  be  expressed  as 

=  Ax]  sin(y4^Wl )  +  Ayx  cos(y^Wl )  (1-52) 

In  general,  for  the  n\h  sensor,  we  have 

A?«‘  =  Ax\  sin(vin)W] )  +  (Ay,  -(; Vm  ~  yn  ) cos(^”)W| ),  n  =  \,...,N  (1.53) 

For  the  multipath  to  focus  at  the  location  xp ,  we  must  have, 

Ag"'  =  0,  forn  =  \,...,N  (1.54) 

This  yields  a  least  squares  (LS)  formulation,  given  by 

A,e,=b,,  e,  =  [Ax,  ,Ay,  ]7,  A,  :=[ap  a"1  ]r 

a,w‘  =[sin(y/,^1)Wl ),..., sin(y/$f)"t)]T,  a"1  =  [cos(^1)Wl ),..., cos(y4A>1  )]r  (1.55) 

•»? 

where  '  o '  denotes  the  Hadamard  or  element-wise  product.  The  solution  of  (1.55)  is  readily  obtained 
by  (AfAjr'b,. 

Now,  considering  the  multipath  w.r.t  wall-2,  we  have  the  following  LS  formulation  for  the  focused 
ghost  pixel. 

A2e2=b2,  e2  =[Ax2,A_y2]r,  A2  :=  [a^2  apf 

a =[sin(^}/)M’2 ),..., sin(ii4A°W2)]7',  a22  =  [cos(^1)W2 ),..., cos(^)W2 )]r  (1.56) 

b2  :=[0,y% -y?x2,...,y% -ytff 

Similarly,  for  the  multipath  from  wall-3,  we  have 

A3e3  =  b3,  e3  =[Ax3,Ay3]r,  A3  :=[ap  apf 

«P  =[sin(^l)W3 ),..., sin(^Ar)M’3)]7',  a23  =  [cos(^1)W3 ),..., cos(^,^K)] 7  (1.57) 

*3t=[0,y%  -yZ,...,y%  -yff 

The  formulations  in  (1.55-1.57)  assume  that  the  sensor  position  increases  from  left  to  right.  In  other 
words,  sensor-1  is  at  the  far  left  of  the  array  whereas  sensor-A/  is  at  the  far  right. 
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1.3.3  Multipath  Exploitation  Algorithm 


Noting  that  the  multipath  ghosts  exist  due  to  the  presence  of  the  target,  we  state  our  objective  as 
follows.  Given  the  beamformed  image  {/(x*)}f=1 ,  our  aim  is  to  associate  each  multipath  ghost  with 
the  respective  target  via  the  model  developed  in  Section  1.3.2.  The  principle  is  captured  in  Fig.  15, 
which  consists  of  two  targets  and  six  false  positives  or  multipath  ghosts.  We  desire  to  associate  and 
map  these  ghosts  to  their  respective  true  target  locations.  The  main  advantages  of  such  an 
association  or  mapping  are  reduction  in  false  positives  in  the  original  beamformed  image,  and  an 
increase  in  the  SCR  at  the  true  target  coordinates.  Note  that  the  first  advantage  is  directly  implied  in 
Fig.  15,  whereas  the  second  advantage  is  explained  as  follows. 

For  ease  of  exposition,  the  technique  is  explained  considering  the  focused  multipath  ghost  from 
wall-1  only.  The  technique  for  exploiting  all  the  ghosts  w.r.t  to  all  three  walls  is  enumerated 

subsequently.  Consider  an  arbitrary  pixel  location  x  =  [x,y]T ,  whose  focused  multipath  ghost  w.r.t 

wall-1  presents  itself  at  location  xW]  =[xW]  ,yW]]T .  Consider  an  intermediate  image,  ^(x),  wherein 

the  association  and  mapping  of  the  focused  multipath  ghosts  is  performed  using  simple  2D 
weighting  functions.  That  is, 

I 1  (x)  €  C  ^  =  X  7(x  A  )0 ,  (x  M 1 ,  x  A ,  cr  2  )0  2  (x,  X  A ,  cr 2 ) 

k=l  (1.58) 

0](xw\xi,<72):=exp(-|xWl  -\k  |2  /cr2)o2(x,xi,cr2):=  1  - exp(-|  x- \k  | 2  /cr2) 

where  <!>](•) and  <I>2(-)are  the  weighting  functions  and  cr2 is  an  arbitrary  variance.  The  weighting 
functions  are  related  to  the  two  dimensional  real  Gaussian  distributions.  A  natural  improvement, 
albeit  not  pursued  here,  is  to  let  the  variance  in  the  y-direction  to  be  equal  to  the  system's  range 
resolution,  the  variance  in  the  x-direction  to  be  equal  to  the  system's  crossrange  resolution,  and 
introducing  a  correlation  which  is  a  function  of  the  actual  ghost  location.  This  would  ensure  that  the 
ensuing  weighting  functions  are  rotated  and  sheared  to  match  the  system's  point  spread  function, 
nevertheless  at  a  higher  computational  cost.  The  role  of  OjQ and  ®2(-)is  explained  as  follows. 
Consider  the  weighting  function  OjQ;  if  the  pixel  x*  is  in  the  vicinity  of  the  focused  multipath 

ghost  pixel  xH  associated  with  the  target  location  x,  then  O^*)  assumes  a  large  value.  On  the 

other  hand,  if  x*  is  not  in  the  vicinity  of  xM],then  OjQ  assumes  a  low  value.  Since  an  exponential 
function  is  used,  the  weights  are  always  between  0  and  1,  the  lower  bound  being  theoretically 

unachievable.  The  vicinity  is  obviously  controlled  by  the  variance  parameter,  cr2.  From  our 
extensive  simulations,  we  found  that  a  lower  variance  must  be  chosen  for  smaller  dimension  rooms. 
For  small  rooms,  it  is  logical  to  assume  that  the  ghosts  and  higher  order  multipath  are  located  quite 
close  to  one  another.  Therefore,  it  becomes  prudent  to  be  conservative  in  selecting  the  variances. 
Large  variances  for  smaller  rooms  lead  to  false  mappings.  Thus,  variance  selection  for  the  problem  at 
hand  must  be  guided  by  the  actual  dimensions  of  the  room  being  imaged.  The  weighting  function 
<D2(-)  rejects  pixel  locations  at  and  in  the  vicinity  of  the  imaged  pixel,  i.e.,  the  location 
corresponding  to  the  assumed  true  target  location  x*  =x.  This  is  required  because,  although  the 
focus  point  is  singular,  but  due  to  the  system's  point  spread  function,  some  of  the  energy  is  spread 
to  the  neighboring  pixels. 

From  our  extensive  simulations,  we  found  that  as  the  original  image  /(•)  is  complex-valued, 
summing  or  integrating  across  the  multipath  pixels  may  result  in  destructive  interference  yielding  a 
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poor  SCR.  This  loss  of  SCR  is  against  our  principle  of  multipath  exploitation.  To  alleviate  this  problem, 
we  use  the  intensity  of  the  original  beamformed  image  in  (1.58)  as 


7,(x)e 3^^  =  XI 7(x*) I2  °i(xW,>x*>c2)<D2(x>x*>°'2) 


(1-59) 


which  implies  an  incoherent  summation.  Now  considering  the  ghosts  w.r.t  the  three  wails,  we  can 
readily  see  that  (1.59)  is  modified  as 


(1.60) 


\Vj=\i\  k= 1 


Image  Z^*)  will  have  no  intensity  at  and  near  the  vicinity  of  the  focused  multipath  pixels.  Hence, 
consider  the  following  composite  image  obtained  by  simple  pixel-wise  multiplication. 


/ /  (x)  =  /i  (x)x  I  /(x)  I2 


(1.61) 


The  effect  of  pixel-wise  multiplication  in  (1.61)  is  as  follows.  The  image  Z^*)  will  have  deep  nulls  at 
the  hypothesized  ghosts  locations,  whereas  a  large  peak  at  the  true  target  location.  Thus,  (1.61)  is 
essentially  a  masking  operation,  nulling  the  ghosts  in  the  original  image  Z(-),  and  simultaneously 
enhancing  the  energy  at  the  true  target  location.  Hence,  the  composite  image  I /(-)  will  alleviate  the 

false  positives,  i.e.,  the  multipath  ghosts  and  increase  the  intensity  of  the  true  target  pixel.  In  the 
approach  described,  the  multipath  ghost  locations  are  readily  obtained  using  the  multipath  focusing 
analysis  in  Section  1.3, 2. 2. 

It  is  noted  that,  in  practice,  different  targets  may  exhibit  significantly  weaker  bistatic  radar  cross 
sections  (RCS)  relative  to  their  monostatic  RCS.  Accordingly,  some  or  maybe  all  of  the  target  ghosts 
may  go  unnoticed.  In  this  case,  the  proposed  multipath  exploitation  approach  will  not  prove 
beneficial  as  there  may  not  be  any  increase  in  the  target  SCR.  Also,  consider  a  scenario  where  one 
target  presents  itself  at  one  of  the  other  target's  focused  multipath  ghost.  Since  multipath  is 
unexploited  in  the  original  beamformed  image,  we  would  declare  the  presence  of  the  target  at  the 
pixel  in  question.  On  the  other  hand,  the  proposed  approach  will  map  this  target  and  of  course  the 
collocated  multipath  ghost  back  to  the  other  target's  location  causing  the  second  target  to  go 
undetected.  The  question  then  arises  whether  one  should  ignore  the  multipath  ghosts  and  consider 
them  as  true  targets  or  instead  exploit  the  multipath  ghosts.  The  former  has  an  unreasonably  high 
false  alarm  rate,  whereas  the  latter  yields  poor  target  detection  for  the  considered  scenario,  which, 
although  rare,  can  occur  in  practice.  However,  noting  that  indoor  targets  of  interest  are  typically  in 
motion,  it  is  reasonable  to  assume  that  a  genuine  target  would  not  persist  on  a  multipath  ghost 
throughout  the  surveillance  period.  Hence,  multipath  exploitation  should  be  preferred.  Another 
scenario  of  interest  is  when  the  direct  path  to  the  target  is  blocked.  In  that  case,  although  the  ghosts 
will  be  mapped  to  the  true  target  coordinates  in  the  image  Ix  (•),  the  masking  operation  in  (1.61)  will 
cause  the  target  to  still  go  undetected  as  the  original  image  ZQ  will  have  a  null  at  the  target  location. 

1.3.4  Experimental  Results 

A  through-the-wall  SAR  system  was  set  up  in  the  Radar  Imaging  Lab  at  Villanova  University.  A 
stepped-frequency  signal  with  696  steps  covering  the  l-3GHz  frequency  band  was  used.  This 
provides  a  range  resolution  of  0.15m  and  a  maximum  unambiguous  range  of  40m.  The  room 
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dimensions  are  Dx  =  5.09m  and  D2  =  3.78m  The  side  walls  and  part  of  the  floor  were  covered  with 
RF  absorbers  and  only  multipath  due  to  the  back  wall  was  considered.  The  front  wall  is  made  of 
0.15m  thick  solid  concrete  blocks  with  a  dielectric  constant  of  7.66.  The  synthetic  aperture  consisted 
of  an  81-element  monostatic  linear  array  with  an  inter-element  spacing  of  2.22cm.  The  standoff 
distance  of  the  array  from  the  front  wall  is  1.53m.  The  experimental  set  up  is  depicted  in  Fig.  16. 

Two  different  experiments  were  conducted.  In  the  first  experiment,  a  0.35m  diameter  metal  sphere, 
labeled  as  Target  V  in  Fig.  16,  was  used  as  a  target.  The  second  experiment  consisted  of  two 
targets,  the  first  being  the  sphere  used  in  the  single  target  case,  while  the  second,  labeled  as  Target 
2'  in  Fig.  16,  was  a  smaller  metallic  sphere  of  0.2032m  diameter.  An  empty  scene  measurement 
without  the  targets  present  in  the  scene  was  also  made  and  was  coherently  subtracted  from  the 
target  scene  data  to  reduce  the  clutter.  The  beamformed  image  corresponding  to  the  single  target 
scene  is  depicted  in  Fig.  17(a);  both  the  target  and  the  ghost  are  clearly  visible  in  the  image.  In 
addition,  other  target-room  interactions,  such  as  those  resulting  from  multiple  bounces  within  the 
front  wall,  can  be  observed.  The  result  after  multipath  exploitation  is  provided  in  Fig.  17(b).  The 
variance  of  the  weighting  functions  in  the  exploitation  algorithm  was  chosen  to  be  equal  to  the 
system's  range  resolution.  We  observe  that  as  expected,  the  ghost  has  been  correctly  mapped  back 
to  the  target,  while  the  other  target-room  interactions  are  left  as  is  in  the  exploited  image.  Figure 
18(a)  shows  the  beamformed  image  corresponding  to  the  two  target  scene.  Both  targets  as  well  as 
the  two  ghosts  with  respect  to  the  back  wall  are  clearly  seen.  The  multipath  exploited  image  is 
provided  in  Fig.  18(b),  in  which  one  can  see  that  the  ghosts  have  been  mapped  back  to  the  target 
locations.  We  note  that,  similar  to  the  single  target  case,  the  other  target-room  and  target-target 
interactions  are  retained  in  Fig.  18(b). 

1.4  Change  Detection  Based  MTi  Approach 

Detection  of  humans  is  one  of  the  most  important  objectives  in  urban  sensing  and  through-the-wall 
radar  technology  [13-15,  28,  46,  50,  58].  Humans  belong  to  the  class  of  animate  objects  which  are 
characterized  by  motion  of  the  limbs,  breathing,  and  heartbeat.  These  features  separate  animate 
and  inanimate  objects  and  allow  the  detection  of  targets  of  interest  to  proceed  based  on  changes  in 
the  phase  of  the  scattered  radar  signals  over  successive  probing  and  data  observations. 

Change  detection  techniques  have  been  recently  used  to  detect  moving  targets  in  the  presence  of 
heavy  clutter  that  is  caused  by  strong  reflections  from  exterior  and  interior  walls.  In  the  case  of 
moving  targets,  the  subtraction  of  two  consecutive  or  non-consecutive  data  frames  or  images  enjoys 
the  same  benefits  as  the  background  subtraction  process  for  the  case  of  stationary  targets.  That  is, 
it  results  in  effective  removal  of  all  wall  returns  and  non-target  clutter.  The  subtraction  operation 
performed  for  either  case  of  stationary  or  moving  targets  is  referred  to  as  Change  Detection  (CD). 
The  length  of  the  time  period  elapsing  between  the  two  datasets  to  be  subtracted  may  differ  for  the 
two  cases.  While  the  stationary  scene  permits  long  time  periods,  the  moving  target  case  necessitates 
short  periods.  When  operating  on  data  frames,  both  operations  can  be  described  by  what  is  known 
as  delay  line  cancellers  (DLC)  [66]. 

It  is  known  that  moving  target  indication  (MTI)  processing  applies  clutter  filters  to  remove  radar 
returns  scattered  from  stationary  objects.  Delay  line  cancellers  can  be  designed  such  that  their 
frequency  responses  place  a  notch  at  DC  and  concurrently  meet  other  desirable  passband  and 
stopband  filter  characteristics  [66].  Doppler  filter  banks  typically  follow  the  delay  line  canceller,  and 
provide  benefits  of  signal  separations,  radial  velocity  measurements,  and  noise  reduction  [66].  It  is 
important  to  note  that  for  urban  sensing  environments,  changes  in  the  backscattered  signal  phase 
due  to  motion  do  not  necessarily  lend  themselves  to  Doppler  frequency  shifts.  The  human  motion 
can  be  abrupt  and  highly  nonstationary,  producing  a  time-dependent  phase  whose  rate  of  change 
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may  fail  to  translate  into  a  single  shift  or  multi-component  sinusoids  that  can  be  captured  by 
different  Doppler  filters.  Instead,  the  corresponding  wide  spectrum  of  human  motions  becomes 
non-localizable  and  can  span  the  entire  radar  frequency  band.  In  lieu  of  Doppler  filters,  time- 
frequency  processing  can  be  applied  to  reveal  the  instantaneous  frequency  signatures,  as  it  is  the 
case  with  Human  gait.  Human  gait  classification  and  biometric  radar  have  been  discussed  in  [55-57, 
72,  75].  However,  apart  from  regularized  motions,  such  as  walking  and  running,  time-frequency 
Doppler  signal  representations  are  very  complex  and  difficult  to  interpret,  especially  when  dealing 
with  non-homogeneous  walls.  Therefore,  the  application  of  Doppler  and  Microdoppler  filters  for 
indoor  target  surveillance  may  not  significantly  add  to  target  motion  detection  that  is  already 
achieved  by  delay  line  cancellers. 

Change  detection  in  through-the-wall  radar  imaging  was  first  discussed  concurrently  in  [59]  and  [52]. 
In  [59],  the  Synchronous  Impulse  Reconstruction  (SIRE)  radar  system,  operating  in  the  frequency 
band  of  300MHz-3GHz  and  employing  a  multi  transmit/receive  design,  was  used  for  MTI.  The  change 
detection  approach  discussed  in  [59]  subtracts  the  down  range  profiles  over  consecutive  data 
frames,  emulating  a  derivative  operator.  It  was  shown  that  the  SIRE  successfully  detects  a  human 
target  moving  within  an  enclosed  structure.  On  the  other  hand,  change  detection  was  applied  in  [52] 
in  the  context  of  background  subtraction  to  detect  stationary  targets  using  data  acquired  during 
interrogations  of  a  scene  at  two  different  time  instants.  Image  intensity  subtraction  was  employed 
and  examples  based  on  ray  tracing  data  covering  the  0.7  to  3.1  GHz  frequency  range  were  provided. 
More  recently,  [69]  applied  change  detection  for  slow-moving  target  detection  in  through-the-wall 
radar  imaging  (TWRI)  by  subtracting  the  image  intensity  corresponding  to  different  data  frames, 
followed  by  a  zero  threshold  operation  to  suppress  the  reference  image.  Examples  of  target 
detection  in  [69]  included  both  EM  modeling  data  and  real  experiments.  Change  detection  using 
first-  and  second-  order  motion  detection  images  for  TWRI  application  were  discussed  in  [41,  90]. 
The  former  has  better  sensitivity  to  motion,  whereas  the  latter  provides  better  signal-to-noise  ratio 
(SNR).  Experimental  data,  reported  in  [41,  90]  using  four  antennas  and  a  signal  bandwidth  of  500 
MHz  to  1  GHz,  demonstrated  that  change  detection  techniques  are  capable  of  detecting  walking 
human  and  simulated  human  breathing  using  machine  motion.  Moving  target  indicator  filtering  was 
also  used  in  [41]  for  detection  of  humans  in  motion  inside  buildings.  An  ultrawideband  frequency- 
modulated  continuous-wave  radar  with  an  extended  frequency  sweep  from  0.5  to  8  GHz  was  used  in 
the  differential  mode  to  track  human  motion  behind  a  brick  wall. 

In  this  section,  we  examine  change  detection  formulation  and  performance  when  clutter  removal  is 
performed  either  coherently  or  noncoherently.  Rather  than  operating  on  successive  pulses,  delay 
line  canceller  equivalences  are  applied  to  different  data  frames  for  each  range  bin  (or  equivalently 
different  complex  amplitude  images  for  each  pixel)  (coherent  CD)  or  different  intensity  images  for 
each  pixel  (noncoherent  CD).  The  frames  can  be  consecutive,  dealing  with  targets  exhibiting  sudden 
short  motions,  or  nonconsecutive  with  relatively  long  time  difference,  for  the  case  in  which  the 
target  changes  its  range  gate  position. 

1.4.1  Signal  Model 

We  develop  the  signal  model  for  wideband  operation  with  multiple  transmitters  and  receivers  under 
coherent  integration.  Sequential  use  of  the  transmitters  with  simultaneous  reception  at  all  receivers 
is  assumed.  It  is  noted  that  time-multiplexing  of  the  transmitters  is  a  viable  option  for  coherent  radar 
operation  in  through-the-wall  radar  sensing  applications.  This  is  due  to  the  following  two  main 
reasons:  a)  Ground-based  urban  imaging  radar  systems  are  typically  vehicle  portable,  which  restricts 
the  available  real  estate  for  deploying  multiple  antennas  on  the  same  platform;  this,  combined  with 
the  important  constraint  of  low  cost,  results  in  deployment  of  far  more  receivers  than  wideband 
transmitters,  and  b)  Indoor  targets  of  interest,  such  as  humans,  move  at  low  velocities.  As  a  result, 
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the  time  required  to  collect  one  data  frame  with  sequential  use  of  the  transmitters  is  small  enough 
so  that  the  target  remains  stationary  during  the  measurement.  The  sequential  transmit  operation  is 
a  salient  feature  of  two  known  through-the-wall  radar  imaging  systems;  one  is  built  by  the  Army 
Research  Lab  [50,  59],  and  the  other  by  the  Defense  Research  and  Development  Canada  [65]. 

Consider  an  imaging  radar  with  an  AT-element  linear  transmit  array  and  an  jV-element  linear  receive 
array.  Let  s(/)  be  the  wideband  signal  used  for  interrogating  the  scene.  The  signal  duration  Ts  is 
assumed  to  be  less  than  the  pulse  repetition  interval  Tr.  Assuming  coherent  integration  of  K  pulses, 
the  coherent  integration  interval  becomes  KTr.  For  the  case  of  a  single  point  target  of  reflectivity 
<jpl  located  at  \p  =  (xp9yp),  the  &th  emitted  pulse  with  the  mth  transmitter  active  is  received  at 
the  «th  receiver  in  the  form 


ZmnO’k)  =  <V(<  “  mKTr  ~  kTr  ~  z p,mn)  (1-62) 

where  rp  mn  is  the  propagation  delay  for  the  signal  to  travel  between  the  mth  transmitter,  the  target 
at  xp,  and  the  nth  receiver.  For  through-the-wall  propagation,  rpmn  will  comprise  the  components 
corresponding  to  traveling  distances  before,  through,  and  after  the  wall  [11].  The  term  mKTr  in 
(1.62)  indicates  the  additional  time  required  for  data  collection  when  the  mth  transmitter  becomes 
active  in  a  sequential  order.  We  assume  that  the  transmitters  transmit  their  pulses  in  contiguous 
time  intervals,  each  of  length  KTr .  It  is  further  assumed  that  we  deal  with  targets  of  constant  radar 
cross  section  (RCS)  over  the  range  of  the  used  frequencies  and  target  aspect  angles.  The  received 
signal  is  passed  through  a  filter  matched  to  the  transmit  waveform  and  the  output  can  be  expressed 
as 


zmn  i‘’k)  =  znw(‘’k)  *  h(t )  =  jJ'  h(u)zmn(t  -  u,k)du 


(1.63) 


where  h{t)  -  With  K  pulses  coherently  integrated,  we  obtain 

K—\ 

*£(')=£  z£it,k)  (1.64) 

*=0 

In  order  to  generate  an  image  of  the  scene  being  interrogated,  the  MN  outputs  corresponding  to  the 
operation  of  M  transmitters  and  N  receivers  are  processed  as  follows.  The  region  of  interest  is 
divided  into  a  finite  number  of  pixels  in  x  and  y,  where  x  and  y  represent  the  crossrange  and 
downrange,  respectively.  The  composite  signal  corresponding  to  the  pixel,  located  at  xq  =  (xqiyq),  is 

obtained  by  summing  time  delayed  versions  of  the  MN  filtered  signals, 

Af-l  AM 

zq  (0  =  Z  'ZZmn(t+Tq,mn)  (1-65) 

m~0  n=Q 

where  rq  mn  is  the  focusing  delay  applied  to  the  matched  filtered  output  of  the  «th  receiver  with  the 

mth  transmitter  active.  It  is  noted  that  additional  weighting  can  be  applied  during  the  summation 
operations  of  (1.65)  to  control  the  sidelobe  level  of  the  transmit-receive  array  pattern  of  the  imaging 
radar  system.  Substituting  (1.62),  (1.63),  and  (1.64)  in  (1.65)  yields 
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(1.66) 


M-\N-\K-\Tr 

Zq(0=  Z  Z  Z  \h(u)zmn(t  +  Tq,mn-U,k)du 

m=0n=0k=0  0 

The  complex  amplitude  image  value  /(x?)for  the  pixel  at  xq  is  obtained  by  sampling  the  signal 
z?(t)attime  t  =  0.  That  is, 


'<x>>=2«<'>L  (1-67) 

The  process  described  by  (1.65)-(1.67)  is  repeated  for  all  pixels  in  the  image  to  generate  the 
composite  image  of  the  scene.  The  general  case  of  multiple  targets  can  be  obtained  by  superposition 
of  target  reflections  [5, 11]. 


1.4.2  Change  Detection 

In  its  simplest  form,  the  MTI  processing  is  achieved  by  using  the  radar  images  corresponding  to  two 
data  frames,  which  may  be  consecutive  or  separated  by  one  or  more  data  frames.  Due  to  time 
multiplexing  of  transmitters,  a  data  frame  length  Tf  is  equal  to  M  coherent  integration  intervals,  i.e. 

Tf-MKTr .  Let  {/^(x^)]^  and  ^t+LTf\\q^ ’_Q  be  the  radar  images  corresponding  to  the  data 

collected  over  the  two  frames,  where  Q  is  the  number  of  pixels  in  each  image  and  L  denotes  the 
number  of  frames  between  the  two  acquisitions.  It  is  noted  that  L  ~  1  represents  the  case  when  the 
two  acquisitions  are  performed  over  consecutive  frames.  The  timing  interval  for  each  data  frame  is 
assumed  to  be  a  fraction  of  a  second  so  that  the  moving  target  appears  stationary  during  data 
collection. 

Change  detection  can  be  performed  in  two  ways,  namely,  coherent  CD  and  noncoherent  CD  [41,  52, 
59,  69,  90].  In  the  former  approach,  complex  amplitude  image  values  corresponding  to  the  two  data 
frames  are  subtracted  on  a  pixel  by  pixel  basis.  The  latter  approach  involves  subtraction  of  image 
pixel  intensity  values  corresponding  to  different  data  frames,  followed  by  a  zero  threshold  operation 
to  suppress  the  reference  image.  For  the  coherent  approach,  the  motion  detection  image  can  be 
expressed  as  follows, 


A(x,  )  =  I('+LT/)(xq  )  -  7(,) (X,  ),  q  =  0,1, . . . ,  Q  - 1 
Using  (66)  and  (67),  we  can  rewrite  (68)  as 


M-\N~\K-\Tr 

A(x?)=  Z  Z  Z  jh(u)Szmn(t  +  Tqmn-u,k)du 

m=0  n= 0  Jc=0  0 


(1.68) 


(1.69) 


where 


(*,k)  =  zmn ( t  -  LMKTr ,k)-zmn(t,k)  =  zmn{t, k )  *[S(t  -  LMKTr  -kTr)-8(t-kTr)\  ( 1. 70) 

and  zmn(t,k) is  defined  in  (1.62).  From  (1.68)-(1.70),  we  observe  that  the  subtraction  of  complex 
amplitude  image  values  for  each  pixel  amounts  to  subtraction  of  the  received  data  corresponding  to 
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the  two  frames  for  each  range  bin,  followed  by  image  formation.  Equation  (1.70)  provides  the 
equivalent  delay  line  canceller  representation  for  change  detection  directly  in  the  data  domain.  A 
block  diagram  of  the  data  domain  CD  algorithm  implementation  is  provided  Fig.  19.  Note  that  the 
windowing  operation  in  the  figure  ensures  application  of  the  algorithm  over  data  frames  of  length 
7y.The  block  diagram  of  Fig.  19  is  a  two-frame  delay  line  canceller,  reminiscent  of  the  commonly 

used  two-pulse  DLC  for  clutter  filter  in  range-Doppler  radars  [66].  More  frames  can  be  used  in  CD 
leading  to  multiple  frames  delay  line  canceller,  which  is  similar  to  three-  or  higher-  pulse  DLC. 
However,  unlike  range-Doppler  radar  systems,  in  which  higher  order  clutter  filters  provide  flexibility 
in  trading  off  clutter  suppression  with  Doppler  frequency  magnitude,  there  are  no  clear  advantages 
of  using  multiple  frames  in  CD  for  the  underlying  problem.  Hunt  suggested  averaging  multiple 
frames,  only  to  emphasize  the  stationary  targets  [90]. 

The  noncoherent  CD  approach  involving  intensity  images  is  depicted  in  Fig.  20.  In  this  case,  the 
intensity  images  corresponding  to  the  two  data  frames  are  first  computed.  Then,  a  difference  image 
is  generated  by  a  pixel  by  pixel  subtraction  of  the  intensity  images  as 


(1.71) 


The  final  motion  detection  image  is  obtained  by  a  zero  threshold  operation 


(1.72) 


This  threshold  results  in  the  target  appearing  in  the  final  motion  detection  image  only  at  the  position 
it  was  occupying  during  acquisition  of  the  second  (test)  data  frame,  since  the  target  image  at 
position  corresponding  to  the  first  (reference)  data  frame  appears  as  a  negative  value  after  intensity 
subtraction  and  is  thus  removed.  It  is  noted  that,  unlike  the  complex  amplitude  image  subtraction, 
there  is  no  equivalent  data  domain  representation  for  the  intensity  image  subtraction  or 
noncoherent  CD  approach  due  to  the  nonlinearity  of  the  modulus  function  in  (1.71).  However,  the 
delay  line  canceller  framework  can  still  be  applied.  The  image  intensity  at  each  pixel  can  be  viewed 
as  a  nonlinear  function  of  the  data  in  the  respective  frame,  and  the  subtraction,  implied  by  the  DLC, 
is  performed  between  two  corresponding  pixels  separated  by  one  or  multiple  frames. 

The  aforementioned  change  detection  schemes,  though  presented  for  motion  detection,  can  also  be 
applied  for  detection  of  stationary  targets  in  the  presence  of  clutter  and  noise  [52].  If  a  target-free 
background  image  of  the  scene  being  interrogated  is  available  (feasible  when  the  radar  system  is 
being  used  for  long  term  surveillance),  it  can  be  subtracted  (complex  amplitude  or  intensity  pixel 
subtraction)  from  an  image  of  the  same  scene  acquired  at  a  later  time  that  may  contain  a  target  of 
interest,  thus  rendering  background  subtraction  as  a  special  case  of  change  detection. 

1.4.3  Experimental  Results 

A  wideband  through-the-wall  imaging  radar  system  was  set  up  in  the  Radar  Imaging  Lab  at  Villanova 
University.  A  stepped-frequency  CW  signal,  consisting  of  101  frequency  steps  of  size  10  MHz, 
covering  the  1.5-2. 5  GHz  band  was  chosen  for  imaging.  The  bandwidth  of  1  GHz  provides  a  range 
resolution  of  15cm  and  the  10  MHz  frequency  step  allows  15m  unambiguous  range,  which  is  roughly 
three  times  the  length  of  the  room  being  imaged.  An  Agilent  network  analyzer,  model  ENA  5071B, 
was  used  for  signal  synthesis  and  data  collection.  A  horn  antenna,  model  ETS-Lindgren  3164-04,  with 
an  operational  bandwidth  from  0.7  to  6  GHz,  was  mounted  on  a  Field  Probe  Scanner  to  synthesize  an 
11-element  uniform  receive  line  array  with  an  inter-element  spacing  of  7.5cm.  Two  horn  antennas 
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(model  ETS-Lindgren  3164-04),  mounted  on  tripods  and  placed  slightly  above  and  on  either  side  of 
the  receive  array,  were  used  as  transmitters.  Each  transmitter  was  spaced  3.75cm  away  from  the 
nearest  receiver.  Thus,  for  the  system  used,  M  ~  2  and  TV  =11.  A  10ft  x  8ft  wall  segment  was 
constructed  utilizing  0.14m  thick  solid  concrete  blocks  with  a  dielectric  constant  of  7.66.  The  receive 
array  was  at  a  standoff  distance  of  1.05m  from  the  wall,  whereas  the  transmitters  were  positioned 
1.33m  in  downrange  from  the  wall. 

In  order  to  illustrate  the  performance  of  change  detection  under  complex  and  intensity  image  delay 
line  cancellers,  we  consider  the  following  experiment.  A  human  walks  away  from  the  front  wall  in  an 
empty  room,  following  a  straight  line  path.  Since  simultaneous  data  cannot  be  collected  at  all 
receive  locations  due  to  the  synthetic  nature  of  the  receive  aperture,  eight  different  target  positions 
were  considered,  separated  by  15cm,  along  the  trajectory,  as  illustrated  in  Fig.  21.  Data  was 
collected  for  all  eight  positions,  with  the  target  stationary  at  each  position,  using  time  multiplexing 
of  the  transmitters.  At  each  position,  coherent  integration  was  performed  over  16  measurements. 
That  is,  K  =  16  for  the  experiment. 

We  first  consider  the  case  where  the  person  moves  from  position  3  to  position  4.  That  is,  the  two 
acquisitions  are  performed  over  consecutive  frames  and  1  =  1.  Figures  22(a)  and  (b)  show  the 
beamformed  images  using  the  data  corresponding  to  positions  3  and  4,  respectively.  Both  images 
show  strong  reflections  from  the  front  wall,  rendering  target  returns  hard  to  detect  in  either  image. 
Figure  23(a)  shows  the  result  of  coherent  change  detection  wherein  the  complex  radar  image  of  Fig. 
22(a)  is  subtracted  from  the  complex  image  of  Fig.  22(b).  We  observe  that  all  stationary  background 
including  the  front  wall  has  been  eliminated  from  the  image,  and  the  target  undergoing  motion  has 
been  detected.  Since  the  target  appears  in  two  consecutive  range  cells,  no  clear  target  displacement 
is  observed  and  the  target  image  appears  blurred.  The  corresponding  noncoherent  CD  results  under 
intensity  image  DLC  are  presented  in  Fig.  23(b).  We  observe  that  the  intensity  image  DLC  removes 
most  of  the  imaged  target  in  the  case  of  consecutive  range  cell  acquisitions.  Moreover,  due  to  the 
nonlinearity  of  the  modulus  function,  the  stationary  background  in  the  image  is  in  general  different 
for  each  data  acquisition,  thereby  causing  the  background  to  be  only  partially  cancelled  out  under 
intensity  image  DLC  [69].  As  a  result,  stronger  artifacts  than  those  of  the  complex  image  DLC  can  be 
clearly  observed.  The  ability  to  detect  the  target  in  the  change  detected  image  is,  therefore,  reduced 
under  intensity  image  DLC  as  compared  to  complex  image  DLC. 

Next,  we  investigate  the  performance  of  the  various  schemes  under  large  target  displacement 
between  consecutive  acquisitions.  Consider  the  case  where  the  person  moves  from  position  1  to 
position  4.  In  this  case,  L  =  3.  Figure  24(a)  shows  the  result  of  coherent  change  detection.  Since  the 
target  undergoes  large  range  migration,  the  moving  target  presents  itself  as  two  targets  in  the 
image.  The  noncoherent  change  detection  results  for  the  large  target  displacement  case  are 
presented  in  Fig.  24(b).  We  observe  that  the  intensity  image  DLC  retains  the  target  at  position  4  and 
removes  the  target  at  position  1.  However,  similar  to  the  small  target  displacement  case,  the 
artifacts  resulting  from  the  background  not  being  totally  cancelled  out  render  target  detection 
difficult. 


1.5  Conclusions 

In  this  section,  we  presented  recent  algorithmic  advances  in  through-the-wall  radar  imaging.  First, 
considering  ground-based  EM  sensing,  we  discussed  two  methods,  namely,  spatial  filtering  and  EM 
modeling  based  solution,  for  mitigating  the  front  wall  return  prior  to  application  of  the  image 
formation  methods.  The  spatial  filtering  technique  builds  on  the  strong  correlation  of  wall  EM 
responses  across  antenna  array  elements  to  reduce  constant-type  return  that  is  typical  of  walls  in 
monostatic  illuminations.  The  spatial  filter,  thus,  allows  the  follow-on  beamforming  to  unmask  and 
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image  behind-the-wall  targets  that  have  limited  spatial  extent,  such  as  humans,  compared  to  walls 
without  the  need  for  o  priori  knowledge  of  the  wall  characteristics.  A  simple  MR  notch  filter  with 
flexible  design  was  compared  with  the  fixed  design  MA  subtraction  filtering,  commonly  used  in  GPR. 
On  the  other  hand,  the  EM  modeling  approach  is  wall-dependent  and  relies  on  accurate  estimation 
of  the  wall  parameters  for  suppression  of  the  wall  return.  Extraction  of  the  wall  parameters  and 
coherent  subtraction  of  the  modeled  wall  reflection  was  shown  to  significantly  improve  the  signal- 
to-  wall-clutter  ratio. 

Second,  an  approach  to  exploit  the  rich  indoor  multipath  environment  for  improved  target  detection 
was  described.  A  ray  tracing  approach  was  used  to  derive  a  multipath  model,  considering  reflections 
inside  an  enclosed  room  comprising  four  homogeneous  walls.  Using  the  model,  it  was  demonstrated 
analytically  that  the  multipath  corresponding  to  each  sensor  appeared  on  the  wall  but  changes 
position  from  one  sensor  to  another.  Hence,  a  least  squares  technique  was  used  to  estimate  its 
actual  focusing  location  in  both  downrange  and  crossrange.  The  model  was  utilized  to  develop  a 
multipath  exploitation  technique  which  associates  multipath  ghosts  with  their  respective  targets  and 
maps  them  to  their  true  target  locations.  This  technique  reduced  the  false  positives  in  the  original 
beamformed  image  as  well  as  increased  the  signal-to-noise  ratio  at  the  true  target  locations. 

Third,  we  discussed  a  change  detection  approach  to  moving  target  indication  for  through-the-wall 
applications.  Change  detection  was  used  to  mitigate  the  heavy  clutter  that  is  caused  by  strong 
reflections  from  exterior  and  interior  walls.  Both  coherent  and  noncoherent  change  detection 
techniques  were  examined  and  their  performance  was  compared  under  both  consecutive  and  non- 
consecutive  acquisitions.  For  non-consecutive  acquisitions,  the  coherent  CD  scheme  showed  two 
sets  of  imaged  targets  corresponding  to  the  positions  of  a  single  target  at  the  two  data  acquisitions, 
whereas  the  noncoherent  CD  retained  only  one  set  of  imaged  targets,  though  with  significantly  more 
artifacts.  For  consecutive  acquisitions,  the  coherent  change  detection  provided  better  performance 
than  the  noncoherent  CD,  which  removed  most  of  the  imaged  target  during  the  zero-thresholding 
step. 
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Figure  1.  Received  signals  using  SAR  data  collection  scheme,  (a)  Array  signal  at  t  =  t0.  (b)  Range 
between  antenna  locations  and  the  target. 
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Figure  2.  Target  return  signal  in  the  time-space  domain,  (a)  Target  is  within  the  array  aperture,  (b) 
Target  is  outside  the  array  aperture 
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Figure  3.  Modification  in  the  PSF  due  to  MA  subtraction,  (a)  Second  term  in  (19).  (b)  PSF  without 
filtering,  (c)  PSF  after  MA  subtraction. 
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Figure  4.  Frequency  response  of  the  notch  filter  for  various  a,  DLC  with  two  pulses,  and  MA 
subtraction  with  67  antennas. 
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Figure  5.  Block  diagram  of  notch  filtering 
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Figure  6.  (a)  Wall  reflections  of  the  solid  concrete  wall  in  time  domain  across  the  antenna  positions, 
(b)  DS  image  with  wall  reflections,  (c)  DS  image  with  background  subtraction. 
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Figure  7.  DS  images  for  solid  concrete  wall,  (a)  MA  subtraction,  (b)  Spatial  notch  filtering. 


Figure  8.  (a)  Side  view  of  the  measurement  setup,  (b)  Trihedral  target. 
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Figure  9.  (a)  Image  of  the  trihedral  behind  wall  using  full  received  data  and  free-space  image 
formation,  (b)  Image  of  trihedral  after  estimated  wall  return  subtraction  and  wall  compensation,  (c) 
Image  of  trihedral  after  wall  compensation. 


Fig.  10.  Through-the-wall  synthetic  aperture  radar  model. 
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Fig.  11.  Multipath  model. 


Fig.  12.  Other  higher  order  paths. 
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Fig.  13.  The  equivalence  between  the  bistatic  and  monostatic  geometries. 
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Fig,  14,  Multipath  locations  and  imaged  location  of  multipath  ghosts. 
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Fig.  15.  Multipath  exploitation  principle. 


Fig.  16.  Scene  Layout 
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Fig.  17.  Single  target  experiment:  (a)  Original,  (b)  after  multipath  exploitation. 
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Fig.  18.  Two  target  experiment:  (a)  Original ,  (b)  after  multipath  exploitation. 
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Figure  19.  Block  diagram  of  the  data  domain  CD  approach,  equivalent  to  complex  amplitude  image 
subtraction  (coherent  CD). 


42 


Figure  20.  Block  diagram  of  the  noncoherent  change  detection  approach. 
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Figure  21.  Layout  of  the  scene  and  target  trajectory. 
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Figure  22.  Beamformed  images  for  the  target  at  position  (a)  3  and  (b)  4  (consecutive  range  cells). 
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Figure  23.  Images  after  change  detection  under  (a)  coherent  CD  and  (b)  noncoherent  CD  for  the 
human  walking  from  position  3  to  position  4. 
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Figure  24.  Images  after  change  detection  under  (a)  coherent  CD  and  (b)  noncoherent  CD  for  the 
human  walking  from  position  1  to  position  4. 


Section  2  -  Compressive  Sensing  for  Through-the-Waii  Radar 
Imaging 

2.1  Introduction 

Through-the-wall  radar  imaging  (TWRI)  is  an  emerging  technology  that  addresses  the  desire  to  see 
inside  buildings  using  electromagnetic  (EM)  waves  for  various  purposes,  including  determining  the 
building  layout,  discerning  the  building  intent  and  nature  of  activities,  locating  and  tracking  the 
occupants,  and  even  identifying  and  classifying  inanimate  objects  of  interest  within  the  building. 
TWRI  is  highly  desirable  for  law  enforcement,  fire  and  rescue,  and  emergency  relief,  and  military 
operations  [1-6]. 


Applications  primarily  driving  TWRI  development  can  be  divided  based  on  whether  information  on 
motions  within  a  structure  or  on  imaging  the  structure  and  its  stationary  contents  is  sought  out.  The 
need  to  detect  motion  is  highly  desirable  to  discern  about  the  building  intent  and  in  many  fire  and 
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hostage  situations.  Discrimination  of  movements  from  background  clutter  can  be  achieved  through 
change  detection  (CD)  or  exploitation  of  Doppler  [7-24].  One-dimensional  motion  detection  and 
localization  systems  employ  a  single  transmitter  and  receiver  and  can  only  provide  range-to-motion, 
whereas  two-  and  three-dimensional  multi-antenna  systems  can  provide  more  accurate  localization 
of  moving  targets.  The  3-D  systems  have  higher  processing  requirements  compared  to  2-D  systems. 
However,  the  third  dimension  provides  height  information,  which  permits  distinguishing  people 
from  animals,  such  as  household  pets.  This  is  important  since  radar  cross  section  alone  for  behind- 
the-wall  targets  can  unreliable. 

Imaging  of  structural  features  and  stationary  targets  inside  buildings  requires  at  least  2-D  and 
preferably  3-D  systems  [25-43].  Because  of  the  lack  of  any  type  of  motion,  these  systems  cannot  rely 
on  Doppler  processing  or  CD  for  target  detection  and  separation.  Synthetic  aperture  radar  (SAR) 
based  approaches  have  been  the  most  commonly  used  algorithms  for  this  purpose.  Most  of  the 
conventional  SAR  techniques  usually  neglect  propagation  distortions  such  as  those  encountered  by 
signals  passing  through  walls  [44].  Distortions  degrade  the  performance  and  can  lead  to  ambiguities 
in  target  and  wall  localizations.  Free-space  assumptions  no  longer  apply  after  the  EM  waves 
propagate  through  the  first  wall.  Without  factoring  in  propagation  effects,  such  as  attenuation, 
reflection,  refraction,  diffraction,  and  dispersion,  imaging  of  contents  within  buildings  will  be 
severely  distorted.  As  such,  image  formation  methods,  array  processing  techniques,  target 
detection,  and  image  sharpening  paradigms  must  work  in  concert  and  be  reexamined  in  view  of  the 
nature  and  specificities  of  the  underlying  sensing  problem. 

In  addition  to  exterior  walls,  the  presence  of  multipath  and  clutter  can  significantly  contaminate  the 
radar  data  leading  to  reduced  system  capabilities  for  imaging  of  building  interiors  and  localization 
and  tracking  of  targets  behind  walls.  The  multiple  reflections  within  the  wall  result  in  wall  residuals 
along  the  range  dimension.  These  wall  reverberations  can  be  stronger  than  target  reflections, 
leading  to  its  masking  and  undetectability,  especially  for  weak  targets  close  to  the  wall  [45]. 
Multipath  stemming  from  multiple  reflections  of  EM  waves  off  the  targets  in  conjunction  with  the 
walls  may  result  in  the  power  being  focused  at  pixels  different  than  those  corresponding  to  the 
target.  This  gives  rise  to  ghosts,  which  can  be  confused  with  the  real  targets  inside  buildings  [46-49]. 
Further,  uncompensated  refraction  through  waits  can  lead  to  localization  or  focusing  errors,  causing 
offsets  and  blurring  of  imaged  targets  [26,  39].  SAR  techniques  and  tomographic  algorithms, 
specifically  tailored  for  TWRI,  are  capable  of  making  some  of  the  adjustments  for  wave  propagation 
through  solid  materials  [26-30,  36-41,  50-57].  While  such  approaches  are  well  suited  for 
shadowing,  attenuation,  and  refraction  effects,  they  do  not  account  for  multipath  as  well  as  strong 
reflections  from  the  front  wall. 

The  problems  caused  by  the  front  wall  reflections  can  be  successfully  tackled  through  wall  clutter 
mitigation  techniques.  Several  approaches  have  been  devised,  which  can  be  categorized  into  those 
based  on  estimating  the  wall  parameters  and  others  incorporating  either  wall  backscattering 
strength  or  invariance  with  antenna  location  [39,  45,  58-61].  In  [39,  58],  a  method  to  extract  the 
dielectric  constant  and  thickness  of  the  non-frequency  dependent  wall  from  the  time-domain 
scattered  field  was  presented.  The  time-domain  response  of  the  wall  was  then  analytically  modeled 
and  removed  from  the  data.  In  [45],  a  spatial  filtering  method  was  applied  to  remove  the  dc 
component  corresponding  to  the  constant-type  radar  return,  typically  associated  with  the  front  wall. 
The  third  method,  presented  in  [59-61],  was  based  not  only  on  the  wall  scattering  invariance  along 
the  array  but  also  on  the  fact  that  wall  reflections  are  relatively  stronger  than  target  reflections.  As  a 
result,  the  wall  subspace  is  usually  captured  in  the  most  dominant  singular  values  when  applying 
singular  value  decomposition  (SVD)  to  the  measured  data  matrix.  The  wall  contribution  can  then  be 
removed  by  orthogonal  subspace  projection. 

Several  methods  have  also  been  devised  for  dealing  with  multipath  ghosts  in  order  to  provide  proper 
representation  of  the  ground  truth.  Earlier  work  attempted  to  mitigate  the  adverse  effects  stemming 
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from  multipath  propagation  [27].  Subsequently,  research  has  been  conducted  to  utilize  the 
additional  information  carried  by  the  multipath  returns.  The  work  in  [49]  considered  multipath 
exploitation  in  TWRI,  assuming  prior  knowledge  of  the  building  layout.  A  scheme  taking  advantage  of 
the  additional  energy  residing  in  the  target  ghosts  was  devised.  An  image  was  first  formed,  the  ghost 
locations  for  each  target  were  calculated,  and  then  the  ghosts  were  mapped  back  onto  the 
corresponding  target.  In  this  way,  the  image  became  ghost-free  with  increased  signal-to-cl utter  ratio 
(SCR). 

More  recently,  the  focus  of  the  TWRI  research  has  shifted  towards  addressing  constraints  on  cost 
and  acquisition  time  in  order  to  achieve  the  ultimate  objective  of  providing  reliable  situational 
awareness  through  high-resolution  imaging  in  a  fast  and  efficient  manner.  This  goal  is  primarily 
challenged  due  to  use  of  wideband  signals  and  large  array  apertures.  Most  radar  imaging  systems 
acquire  samples  in  frequency  (or  time)  and  space,  and  then  apply  compression  to  reduce  the 
amount  of  stored  information.  This  approach  has  three  inherent  inefficiencies.  First,  as  the  demands 
for  high  resolution  and  more  accurate  information  increase,  so  does  the  number  of  data  samples  to 
be  recorded,  stored,  and  subsequently  processed.  Second,  there  are  significant  data  redundancies 
not  exploited  by  the  traditional  sampling  process.  Third,  it  is  wasteful  to  acquire  and  process  data 
samples  that  will  be  discarded  later.  Further,  producing  an  image  of  the  indoor  scene  using  few 
observations  can  be  logistically  important,  as  some  of  the  measurements  in  space  and  frequency  or 
time  can  be  difficult,  unavailable,  or  impossible  to  attain. 


Towards  the  objective  of  providing  timely  actionable  intelligence  in  urban  environments,  the 
emerging  compressive  sensing  (CS)  techniques  have  been  shown  to  yield  reduced  cost  and  efficient 
sensing  operations  that  allow  super-resolution  imaging  of  sparse  behind-the-wall  scenes  [10,  62-76]. 
CS  is  a  very  effective  technique  for  scene  reconstruction  from  a  relatively  small  number  of  data 
samples  without  compromising  the  imaging  quality  [77-89].  In  general,  the  minimum  number  of 
data  samples  or  sampling  rate  that  is  required  for  scene  image  formation  is  governed  by  the  Nyquist 
theorem.  Flowever,  when  the  scene  is  sparse,  CS  provides  very  efficient  sampling,  thereby 
significantly  decreasing  the  required  volume  of  data  collected. 


In  this  section,  we  focus  on  CS  for  TWRI  and  present  norm  reconstruction  techniques  that  address 
the  unique  challenges  associated  with  fast  and  efficient  imaging  in  urban  operations.  Sections  2.2  to 
2.5  deal  with  imaging  of  stationary  scenes,  whereas  moving  target  localization  is  discussed  in 
Sections  2.6  and  2.7.  More  specifically,  Section  2.2  deals  with  CS  based  strategies  for  stepped- 
frequency  based  radar  imaging  of  sparse  stationary  scenes  with  reduced  data  volume  in  spatial  and 
frequency  domains.  Prior  and  complete  removal  of  clutter  is  assumed,  which  renders  the  scene 
sparse.  Section  23  presents  CS  solutions  in  the  presence  of  front  wall  clutter.  Wall  mitigation  in 
conjunction  with  application  of  CS  is  presented  for  the  case  when  the  same  reduced  frequency  set  is 
used  from  all  of  the  employed  antennas.  Section  2.4  considers  imaging  of  the  building  interior 
structures  using  a  CS-based  approach,  which  exploits  prior  information  of  building  construction 
practices  to  form  an  appropriate  sparse  representation  of  the  building  interior  layout.  Section  2.5 
presents  CS  based  multipath  exploitation  technique  to  achieve  good  image  reconstruction  in  rich 
multipath  indoor  environments  from  few  spatial  and  frequency  measurements.  Section  2.6  deals 
with  joint  localization  of  stationary  and  moving  targets  using  CS  based  approaches,  provided  that  the 
indoor  scene  is  sparse  in  both  stationary  and  moving  targets.  Section  2.7  discusses  a  sparsity-based 
CD  approach  to  moving  target  indication  for  TWRI  applications,  and  deals  with  cases  when  the  heavy 
clutter  caused  by  strong  reflections  from  exterior  and  interior  walls  reduces  the  sparsity  of  the 
scene.  Concluding  remarks  are  provided  in  Section  2.8.  It  is  noted  that  for  the  sake  of  not  overly 
complicating  the  notation,  some  symbols  are  used  to  indicate  different  variables  over  different 
sections.  Flowever,  for  those  cases,  these  variables  are  redefined  to  reflect  the  change. 
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2.2  CS  Strategies  in  Frequency  and  Spatial  Domains  For  TWRI 

In  this  section,  we  apply  CS  to  through-the-wall  imaging  of  stationary  scenes,  assuming  prior  and 
complete  removal  of  the  front  wall  clutter  [62,  63].  For  example,  if  the  reference  scene  is  known, 
then  background  subtraction  can  be  performed  for  removal  of  wall  clutter,  thereby  improving  the 
sparsity  of  the  behind-the-wall  stationary  scene.  We  assume  stepped-frequency  based  SAR 
operation.  We  first  present  the  through-the-wall  signal  model,  followed  by  a  description  of  the 
sparsity-based  scene  reconstruction,  highlighting  the  key  equations.  It  is  noted  that  the  problem 
formulation  can  be  modified  in  a  straightforward  manner  for  pulsed  operation  and  multistatic 
systems. 

2.2.1  Through-the-Wall  Signal  Model 


Consider  a  homogeneous  wall  of  thickness  d  and  dielectric  constant  s  located  along  the  x-axis,  and 
the  region  to  be  imaged  located  beyond  the  wall  along  the  positive  z-axis.  Assume  that  an  A/-element 
line  array  of  transceivers  is  located  parallel  to  the  wall  at  a  standoff  distance  z0ff,  as  shown  in  Fig.  1. 

Let  the  nth  transceiver,  located  at  \n  =(xn,-zoff), illuminate  the  scene  with  a  stepped-frequency 
signal  of  M  frequencies,  which  are  equispaced  over  the  desired  bandwidth  0)M_ l  —  C0Q, 

com  =  co0  +  mAco,  m  =  (2.1) 

where  co0  is  the  lowest  frequency  in  the  desired  frequency  band  and  A  CO  is  the  frequency  step  size. 

The  reflections  from  any  targets  in  the  scene  are  measured  only  at  the  same  transceiver  location. 
Assuming  the  scene  contains  P  point  targets  and  the  wall  return  has  been  completely  removed,  the 
output  of  the  nth  transceiver  corresponding  to  the  mth  frequency  is  given  by, 

p-\ 

y(m,n)=  'ZtTpexp(-jamT  )  (2.2) 

P= o 


where  &p  is  the  complex  reflectivity  of  the  pth  target,  and  T  is  the  two-way  traveling  time 

between  the  nth  antenna  and  the  target.  It  is  noted  that  the  complex  amplitude  due  to  free-space 
path  loss,  wall  reflection/transmission  coefficients  and  wall  losses,  is  assumed  to  be  absorbed  into 
the  target  reflectivity.  The  propagation  delay  Tp  n  is  given  by  [27-28,  40] 


 np, air,} 


P,n 


+  - 


21 


np,wall 


+  - 


21, 


np,air,2 


(2.3) 


C  V  c 

where  c  is  the  speed  of  light  in  free-space,  V  =  c/\[s  is  the  speed  through  the  wall,  and  the 
variables  lnp>afr^  lnp,wait’  an^  ^,mr,2rePresent  the  traveling  distances  of  the  signal  before,  through, 


and  beyond  the  wall,  respectively,  from  the  nth  transceiver  to  the  pth  target. 

An  equivalent  matrix-vector  representation  of  the  received  signals  in  (2.2)  can  be  obtained  as 
follows.  Assume  that  the  region  of  interest  is  divided  into  a  finite  number  of  pixels  N.xN  in 
crossrange  and  downrange,  and  the  point  targets  occupy  no  more  than  P  («  Nx  x  Nz )  pixels.  Let 


r(k,l),  k  =  0,1,  —  ? Nx -1,  /  =  0,1,...,JV2  — l,be  a  weighted  indicator  function,  which  takes  the 


value  <Jp  if  the  pth  point  target  exists  at  the  ( k ,  /)-th  pixel;  otherwise,  it  is  zero.  With  the  values 
r{kj )  lexicographically  ordered  into  a  column  vector  r  of  length  NXN „  the  received  signal 
corresponding  to  the  nth  antenna  can  be  expressed  in  matrix-vector  form  as, 


y» = 

where  'F n  is  a  matrix  of  dimensions  MxNxNz,  and  its  mth  row  is  given  by, 


M.= 


00,0 


~j°)mr{NxNz-\),n 


Considering  the  measurement  vector  corresponding  to  all  N  antennas,  defined  as, 


(2.4) 

(2.5) 
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(2.6) 


y= 


T  T  T 

y0  y i  -  y^-i 


the  relationship  between  y  and  r  is  given  by 


y  =  ^r 


where 


VP  = 


IT 


VL ft  \JJj 

T  0  T  1  A  V-l 


'P 


(2.7) 

(2.8) 


The  matrix  ¥  is  a  linear  mapping  between  the  full  data  y  and  the  sparse  vector  r. 


2.2.2  Sparsity-Based  Data  Acquisition  and  Scene  Reconstruction 

The  expression  in  (2.7)  involves  the  full  set  of  measurements  made  at  the  N  array  locations  using  the 
M  frequencies.  For  a  sparse  scene,  it  is  possible  to  recover  r  from  a  reduced  set  of  measurements. 
Consider  y,  which  is  a  vector  of  length  QXQ2  («MV)  consisting  of  elements  chosen  from  y  as 
follows, 

y  =  0y  =  0¥r  (2.9) 

where  O  is  a  QXQ2  xMVmatrix  of  the  form, 

0>  =  kron(5,I0i)-diag{^(O),^(l),...,^(W‘I)}  (2.10) 


In  (2.10),  'kron'  denotes  the  Kronecker  product,  I  is  a  QxxQ2  identity  matrix,  3  is  a  Q2xN 
measurement  matrix  constructed  by  randomly  selecting  Q2  rows  of  an  N  x  identity  matrix,  and 

(p^n\  n  =  0,1,...,#  — 1,  is  a  Qx  xMmeasurement  matrix  constructed  by  randomly  selecting  Qx 
rows  of  an  MxM identity  matrix.  We  note  that  3  determines  the  reduced  antenna  locations, 
whereas  determines  the  reduced  set  of  frequencies  corresponding  to  the  nth  antenna  location. 
The  number  of  measurements  QQ2  required  to  achieve  successful  CS  reconstruction  highly 
depends  on  the  coherence  between  0  and  XF.  For  the  problem  at  hand,  O  is  the  canonical  basis 
and  'P  is  similar  to  the  Fourier  basis,  which  have  been  shown  to  exhibit  maximal  incoherence 
[80].  Given  y,  we  can  recover  r  by  solving  the  following  equation1, 

?  =  argminlrlL  subject  to  y  «  04^  (2.11) 

ii  ii/t 


We  note  that  the  problem  in  (2,11)  can  be  solved  using  convex  relaxation,  greedy  pursuit,  or 
combinatorial  algorithms  [91-96].  In  this  section,  we  consider  Orthogonal  Matching  Pursuit  (OMP), 
which  is  known  to  provide  a  fast  and  easy  to  implement  solution.  Moreover,  OMP  is  better  suited 
when  frequency  measurements  are  used  [95].  It  is  noted  that  the  number  of  iterations  of  the  OMP  is 
usually  associated  with  the  level  of  sparsity  of  the  scene.  In  practice,  this  piece  of  information  is 
often  unavailable  a  priori  and  the  stopping  condition  is  heuristic.  Underestimating  the  sparsity  would 
result  in  the  image  not  being  completely  reconstructed  (underfitting),  while  overestimation  would 
cause  some  of  the  noise  being  treated  as  signal  (overfitting).  Use  of  cross  validation  (CV)  has  been 
also  proposed  to  determine  the  stopping  condition  for  the  greedy  algorithms  [97-99].  CV  is  a 
statistical  technique  that  separates  a  data  set  into  a  training  set  and  a  cross  validation  set.  The 
training  set  is  used  to  detect  the  optimal  stopping  iteration.  There  is,  however,  a  tradeoff  between 
allocating  the  measurements  for  reconstruction  or  CV.  More  details  can  be  found  in  [97,  98]. 


1  Ideally,  minimization  of  the  l0  norm  would  provide  the  sparsest  solution.  Unfortunately,  it  is  NP-hard  to  solve 
the  resulting  minimization  problem.  The  lx  norm  has  been  shown  to  serve  as  a  good  surrogate  for  l0  norm  [90]. 
The  ^-minimization  problem  is  convex,  which  can  be  solved  in  polynomial  time. 
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2.2.3  Illustrative  Results 


A  through-the-wall  wideband  SAR  system  was  set  up  in  the  Radar  Imaging  Lab  at  Villanova 
University.  A  67-element  line  array  with  an  inter-element  spacing  of  0.0187m,  located  along  the  x- 
axis,  was  synthesized  parallel  to  a  0.14m  thick  solid  concrete  wall  of  length  3.05m  and  at  a  standoff 
distance  equal  to  1.24m.  A  stepped-frequency  signal  covering  the  1-3  GHz  frequency  band  with  a 
step  size  of  2.75MHz  was  employed.  Thus,  at  each  scan  position,  the  radar  collects  728  frequency 
measurements.  A  vertical  metal  dihedral  was  used  as  the  target  and  was  placed  at  (0,  4.4)m  on  the 
other  side  of  the  front  wall.  The  size  of  each  face  of  the  dihedral  is  0.39m  by  0.28m.  The  back  and 
the  side  walls  of  the  room  were  covered  with  RF  absorbing  material  to  reduce  clutter.  The  empty 
scene  without  the  dihedral  target  present  was  also  measured  to  enable  background  subtraction  for 
wall  clutter  removal. 

The  region  to  be  imaged  is  chosen  to  be  4.9m  x  5.4m  centered  at  (0,  3.7)m  and  divided  into  33  x  73 
pixels,  respectively.  For  CS,  20%  of  the  frequencies  and  51%  of  the  array  locations  were  used,  which 
collectively  represent  10.2%  of  the  total  data  volume.  Figs.  2(a)  and  2(c)  depict  the  images 

corresponding  to  the  full  dataset  obtained  with  backprojection  and  lx  norm  reconstruction, 
respectively.  Figs.  2(b)  and  2(d)  show  the  images  corresponding  to  the  measured  scene  obtained 
with  backprojection  and  lx  norm  reconstruction,  respectively,  applied  to  the  reduced  background 
subtracted  dataset.  In  Figure  2  and  all  subsequent  figures  in  this  section,  we  plot  the  image  intensity 
with  the  maximum  intensity  value  in  each  image  normalized  to  OdB.  The  true  target  position  is 
indicated  with  a  solid  red  rectangle.  We  observe  that,  with  the  availability  of  the  empty  scene 
measurements,  background  subtraction  renders  the  scene  sparse  and  thus,  CS  based  approach 
generates  an  image  using  reduced  data  where  the  target  can  be  easily  identified.  On  the  other  hand, 
backprojection  applied  to  reduced  dataset  results  in  performance  degradation,  indicated  by  the 
presence  of  many  artifacts  in  the  corresponding  image.  OMP  was  used  to  generate  the  CS  images. 
For  this  particular  example,  the  number  of  OMP  iterations  was  set  to  5. 

2.3  Effects  of  Walls  on  Compressive  Sensing  Solutions 

The  application  of  CS  for  TWRI  as  presented  in  Section  II  assumed  prior  and  complete  removal  of 
front  wall  EM  returns.  Without  this  assumption,  strong  wall  clutter,  which  extends  along  the  range 
dimension,  reduces  the  sparsity  of  the  scene  and,  as  such,  impedes  the  application  of  CS  [71-73]. 
Having  access  to  the  background  scene  is  not  always  possible  in  practical  applications.  In  this  section, 
we  apply  joint  CS  and  wall  mitigation  techniques  using  reduced  data  measurements.  In  essence,  we 
address  wall  clutter  mitigations  in  the  context  of  CS. 

There  are  several  approaches,  which  successfully  mitigate  the  front  wall  contribution  to  the  received 
signal  [39,  45,  58-61],  These  approaches  were  originally  introduced  to  work  on  the  full  data  volume 
and  did  not  account  for  reduced  data  measurements  especially  randomly.  We  examine  the 
performance  of  the  subspace  projection  wall  mitigation  technique  [60]  in  conjunction  with  sparse 
image  reconstruction.  Only  a  small  subset  of  measurements  is  employed  for  both  wall  clutter 
reduction  and  image  formation.  We  consider  the  case  where  the  same  subset  of  frequencies  is  used 
for  each  employed  antenna.  Wall  clutter  mitigation  under  use  of  different  frequencies  across  the 
employed  antennas  is  discussed  in  [68,  73].  It  is  noted  that,  although  not  reported  in  this  section,  the 
spatial  filtering  based  wall  mitigation  scheme  [45]  in  conjunction  with  CS  provides  a  similar 
performance  to  the  subspace  projection  scheme  [73]. 
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2.3.1  Wall  Clutter  Mitigation 


We  first  extend  the  through-the-wall  signal  model  of  (2.2)  to  include  the  front  wall  return.  Without 
the  assumption  of  prior  wall  return  removal,  the  output  of  the  nth  transceiver  corresponding  to  the 
mth  frequency  for  a  scene  of  P  point  targets  is  given  by, 

p- 1 

y(m,n)  =  exp(->mrj  + exp  (rjo>mTpn)  (2.12) 

p=  0 

where  <Jw  is  the  complex  reflectivity  of  the  wall,  and  rw  is  the  two-way  traveling  time  of  the  signal 
from  the  nth  antenna  to  the  wall,  and  is  given  by 

-  (2.13) 

c 

It  is  noted  that  both  the  target  and  wall  reflectivities  in  (2.12)  are  assumed  to  be  independent  of 
frequency  and  aspect  angle.  Many  of  the  walls  and  indoor  targets,  including  humans,  have 
dependency  of  their  reflection  coefficients  on  frequency,  which  could  also  be  a  function  of  angle  and 
polarization.  This  dependency,  if  neglected,  could  be  a  source  of  error.  The  latter,  however,  can  be 
tolerated  for  relatively  limited  aperture  and  bandwidth.  Further  note  that  we  assume  a  simple  scene 
of  P  point  targets  behind  a  front  wall.  The  model  can  be  extended  to  incorporate  returns  from  more 
complex  scenes  involving  multiple  walls  and  room  corners.  These  extensions  are  discussed  in  later 
sections. 


From  (2.12),  we  note  that  rw  does  not  vary  with  the  antenna  location  since  the  array  is  parallel  to 

the  wall.  Furthermore,  as  the  wall  is  homogeneous  and  assumed  to  be  much  larger  than  the 
beamwidth  of  the  antenna,  the  first  term  in  (2.12)  assumes  the  same  value  across  the  array 
aperture.  Unlike  tw  ,  the  time  delay  r n ,  given  by  (2.3),  is  different  for  each  antenna  location,  since 

the  signal  path  from  the  antenna  to  the  target  is  different  from  one  antenna  to  the  other. 


The  signals  received  by  the  N  antennas  at  the  M  frequencies  are  arranged  into  an  M  x  N  matrix,  Y, 


A 


y0 


y„ 


y*-: 


(2.14) 


where  yn  is  the  M  xl  vector  containing  the  stepped-frequency  signal  received  by  the  nth  antenna. 


y„  = 


y(0,n) 


y(m,n) 


y(M-\,n) 


(2.15) 


with  y(m,n)  given  by  (12).  The  eigen-structure  of  the  imaged  scene  is  obtained  by  performing  the 
SVD  of  Y, 

Y  =  VAVh  (2.16) 

where  ‘H'  denotes  the  Hermitian  transpose,  U  and  V  are  unitary  matrices  containing  the  left  and 
right  singular  vectors,  respectively,  and  A  is  a  diagonal  matrix 


(2.17) 


and  Al>Z2>...>AN  are  the  singular  values.  Without  loss  of  generality,  the  number  of 
frequencies  are  assumed  to  exceed  the  number  of  antenna  locations,  i.e.,  M  >  N.  The  subspace 
projection  method  assumes  that  the  wall  returns  and  the  target  reflections  lie  in  different 
subspaces.  Therefore,  the  first  K  dominant  singular  vectors  of  the  Y  matrix  are  used  to  construct  the 
wall  subspace, 
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(2.18) 
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Kail  =Z,U/V/ 
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Methods  for  determining  the  dimensionality  K  of  the  wall  subspace  have  been  reported  in  [59,  60]. 
The  subspace  orthogonal  to  the  wall  subspace  is, 
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wall 


(2.19) 


where  I  is  the  identity  matrix.  To  mitigate  the  wall  returns,  the  data  matrix  Y  is  projected  on  the 
orthogonal  subspace  [60], 


Y  = 


(2.20) 


The  resulting  data  matrix  has  little  or  no  contribution  from  the  front  wall. 


2.3.2  Joint  Wall  Mitigation  and  CS 

Subspace  projection  method  for  wall  clutter  reduction  relies  on  the  fact  that  the  wall  reflections  are 
strong  and  assume  very  close  values  at  the  different  antenna  locations.  When  the  same  set  of 
frequencies  is  employed  for  all  employed  antennas,  the  condition  of  spatial  invariance  of  the  wall 
reflections  is  maintained  [72,  73].  This  permits  direct  application  of  the  subspace  projection  method 
as  a  preprocessing  step  to  the  lx  norm  based  scene  reconstruction  of  (2.11). 


2.3.3  Illustrative  Results 


We  consider  the  same  experimental  setup  as  in  Section  III.C.  Fig.  3(a)  shows  the  result  obtained  with 
lx  norm  reconstruction  using  10.2%  of  the  raw  data  volume  without  background  subtraction.  The 
number  of  OMP  iterations  was  set  to  100.  Comparing  Fig.  3(a)  and  the  corresponding  background 
subtracted  image  of  Fig.  2(d),  it  is  evident  that  in  the  absence  of  access  to  the  background  scene,  the 
wall  mitigation  techniques  must  be  applied,  as  a  preprocessing  step,  prior  to  CS  in  order  to  detect 
the  targets  behind  the  wall. 

First,  we  consider  the  case  when  the  same  set  of  reduced  frequencies  is  used  for  a  reduced  set  of 
antenna  locations.  We  employ  only  10.2%  of  the  data  volume,  i.e.,  20%  of  the  available  frequencies 
and  51%  of  the  antenna  locations.  The  subspace  projection  method  is  applied  to  a  Y  matrix  of 
reduced  dimension  146  x  34 .  The  corresponding  lx  norm  reconstructed  image  obtained  with  OMP  is 
depicted  in  Fig,  3(b),  It  is  clear  that,  even  when  both  spatial  and  frequency  observations  are 
reduced,  the  joint  application  of  wall  clutter  mitigation  and  CS  techniques  successfully  provides  front 
wall  clutter  suppression  and  unmasking  of  the  target. 

2.4  Designated  Dictionary  for  Wall  Detection 

In  this  section,  we  address  the  problem  of  imaging  building  interior  structures  using  a  reduced  set  of 
measurements.  We  consider  interior  walls  as  targets  of  interest  and  attempt  to  reveal  the  building 
interior  layout  based  on  CS  techniques.  We  note  that  construction  practices  suggest  the  exterior  and 
interior  walls  to  be  parallel  or  perpendicular  to  each  other.  This  enables  sparse  scene 
representations  using  a  dictionary  of  possible  wall  orientations  and  locations  [76].  Conventional  CS 
recovery  algorithms  can  then  be  applied  to  reduced  number  of  observations  to  recover  the  positions 
of  various  walls,  which  is  a  primary  goal  in  TWRI. 

2.4.1  Signal  Model  under  Multiple  Parallel  Walls 

Considering  a  monostatic  stepped-frequency  SAR  system  with  N  antenna  positions  located  parallel 
to  the  front  wall,  as  shown  in  Fig.  1,  we  extend  the  signal  model  in  (2.12)  to  include  reflections  from 


51 


multiple  parallel  interior  walls,  in  addition  to  the  returns  from  the  front  wall  and  the  P  point  targets. 
That  is,  the  received  signal  at  the  nth  antenna  location  corresponding  to  the  mth  frequency  can  be 
expressed  as, 

p-\  4,-1 

y(m, n)  =  crw exp(-jcomTw)+  E^exp (-j(DmTpn)+  exp(-/fi>mrW( )  (2.21) 

p= 0  ;=0 

where  Iw  is  the  number  of  interior  walls  parallel  to  the  array  axis,  tw  represents  the  two-way 
traveling  time  of  the  signal  from  the  nth  antenna  to  the  fth  interior  wall  and  crw  is  the  complex 
reflectivity  of  the  fth  interior  wall.  Similar  to  the  front  wall,  the  delays  tw  are  independent  of  the 
variable  n,  as  evident  in  the  subscripts. 

Note  that  the  above  model  contains  contributions  only  from  interior  walls  parallel  to  the  front  wall 
and  the  antenna  array.  This  is  because,  due  to  the  specular  nature  of  the  wall  reflections,  a  SAR 
system  located  parallel  to  the  front  wall  will  only  be  able  to  receive  direct  returns  from  walls  which 
are  parallel  to  the  front  wall.  The  detection  of  perpendicular  walls  is  possible  by  concurrently 
detecting  and  locating  the  canonical  scattering  mechanism  of  corner  features  created  by  the 
junction  of  walls  of  a  room  or  by  having  access  to  another  side  of  the  building.  Extension  of  the 
signal  model  to  incorporate  corner  returns  is  reported  in  [76]. 


Instead  of  the  point-target  based  sensing  matrix  described  in  (2.7),  where  each  antenna  accumulates 
the  contributions  of  all  the  pixels,  we  use  an  alternate  sensing  matrix,  proposed  in  [68],  to  relate  the 
scene  vector,  r,  and  the  observation  vector,  y.  This  matrix  underlines  the  specular  reflections 
produced  by  the  walls.  Due  to  wall  specular  reflections,  and  since  the  array  is  assumed  parallel  to 
the  front  wall  and,  thus,  parallel  to  interior  walls,  the  rays  collected  at  the  nth  antenna  will  be 
produced  by  portions  of  the  walls  that  are  only  in  front  of  this  antenna  (see  Fig.  4(a)).  The  alternate 
matrix,  therefore,  only  considers  the  contributions  of  the  pixels  that  are  located  in  front  of  each 
antenna.  In  so  doing,  the  returns  of  the  walls  located  parallel  to  the  array  axis  are  emphasized.  As 
such,  it  is  most  suited  to  the  specific  building  structure  imaging  problem,  wherein  the  signal  returns 
are  mainly  caused  by  EM  reflections  of  exterior  and  interior  walls.  The  alternate  linear  model  can  be 
expressed  as 

y  =  Tr  (2.22) 


where 


*-hr 


with  x¥n  defined  as, 


kl  = 


’[(do),*]6 


~jcomT(,  0,0) 


’Ktf.-uw.-fti.]* 


~Jcomr(Nx-\,Nz-\) 


(2.23) 

(2.24) 


In  (2.24),  Tk  /  is  the  two-way  signal  propagation  time  associated  with  the  downrange  of  the  ( k ,  /)th 


pixel,  and  the  function  works  as  an  indicator  function  in  the  following  way, 

_Jl  if  the  (k,  /)th  pixel  is  in  front  of  the  «th  antenna 
[0  otherwise  2"i 

That  is,  if  xk  andx„  represent  the  crossrange  coordinates  of  the  (k,  /)th  pixel  and  the  nth  antenna 
location,  respectively,  and  dx  is  the  crossrange  sampling  step,  then  7)  =1  provided  that 

xk-dx/2  <x„<xk  +dx/2  (see  Fig.  4(b)). 
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2.4.2  Sparsifying  Dictionary  for  Wall  Detection 


Since  the  number  of  parallel  walls  is  typically  much  smaller  compared  to  the  downrange  extent  of 
the  building,  the  decomposition  of  the  image  into  parallel  walls  can  be  considered  as  sparse.  Note 
that  although  other  indoor  targets,  such  as  furniture  and  humans,  may  be  present,  their  projections 
onto  the  horizontal  lines  are  expected  to  be  negligible  compared  to  those  of  the  walls. 

In  order  to  obtain  a  linear  matrix-vector  relation  between  the  scene  and  the  horizontal  projections, 
we  define  a  sparsifying  matrix  R  composed  of  possible  wall  locations.  Specifically,  each  column  of  the 
dictionary  R  represents  an  image  containing  a  single  wall  of  length  lx  pixels,  located  at  a  specific 
crossrange  and  at  a  specific  downrange  in  the  image.  Consider  the  crossrange  to  be  divided  into  k 
non-overlapping  blocks  of  lx  pixels  each  (see  Fig.  5(a)),  and  the  downrange  division  defined  by  the 
pixel  grid.  The  number  of  blocks  Nc  is  determined  by  the  value  of  lx,  which  is  the  minimum 
expected  wall  length  in  the  scene.  Therefore,  the  dimension  of  R  is  NXNZ  x  NCNZ,  where  the 
product  NCNZ  denotes  the  number  of  possible  wall  locations.  Figure  5(b)  shows  a  simplified 

scheme  of  the  sparsifying  dictionary  generation.  The  projection  associated  with  each  wall  location  is 
given  by, 

0  2>(M)  (2.26) 

lx  keB[b ] 

where  B[b\  indicates  the  bth  crossrange  block  and  h-  1.2  Defining 

mg  =  [g(1)(0)  -  gW)(0)  g(1)(i)  -  gw)( i)  -  gm(N:- 1)  -  *<*•>(*, -i)J 

the  linear  system  of  equations  relating  the  observed  data  y  and  the  sparse  vector  g  is  given  by, 

y  =  ^Rg  (2.28) 

In  practice  and  by  the  virtue  of  collecting  signal  reflections  corresponding  to  the  zero  aspect  angle, 
any  interior  wall  outside  the  synthetic  array  extent  will  not  be  visible  to  the  system.  Finally,  the  CS 
image  in  this  case  is  obtained  by  first  recovering  the  projection  vector  g  using  lx  norm  minimization 
with  a  reduced  set  of  measurements  and  then  forming  the  product  Rg. 

It  is  noted  that  we  are  implicitly  assuming  that  the  extents  of  the  walls  in  the  scene  are  integer 
multiples  of  the  block  of  lx  pixels.  In  case  this  condition  is  not  satisfied,  the  maximum  error  in 

determining  the  wall  extent  will  be  at  most  equal  to  the  chosen  block  size.  Note  that  incorporation 
of  the  corner  effects  will  help  resolve  this  issue,  since  the  localization  of  corners  will  identify  the  wall 
extent  [76]. 

2.4.3  Illustrative  Results 

A  through-the-wall  SAR  system  was  set  up  in  the  Radar  Imaging  Lab,  Villanova  University.  A  stepped- 
frequency  signal  consisting  of  335  frequencies  covering  the  1  to  2  GHz  frequency  band  was  used  for 
interrogating  the  scene.  A  monostatic  synthetic  aperture  array,  consisting  of  71-element  locations 
with  an  inter-element  spacing  of  2.2cm,  was  employed.  The  scene  consisted  of  two  parallel  plywood 
walls,  each  2.25cm  thick,  1.83m  wide,  and  2.43m  high.  Both  walls  were  centered  at  Om  in 
crossrange.  The  first  and  the  second  walls  were  located  at  respective  distances  of  3.25m  and  5.1m 
from  the  antenna  baseline.  Figure  6(a)  depicts  the  geometry  of  the  experimental  scene. 

The  region  to  be  imaged  is  chosen  to  be  5.65m  (crossrange)  x  4.45m  (downrange),  centered  at  (0, 
4.23)m,  and  is  divided  into  128  x  128  pixels.  For  the  CS  approach,  we  use  a  uniform  subset  of  only  84 
frequencies  at  each  of  the  18  uniformly  spaced  antenna  locations,  which  represent  6.4%  of  the  full 
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data  volume.  The  CS  reconstructed  image  is  shown  in  Fig.  6(b).  We  note  that  the  proposed  algorithm 
was  able  to  reconstruct  both  walls.  However,  it  can  be  observed  in  Fig.  6(b)  that  ghost  walls  appear 
immediately  behind  each  true  wall  position.  These  ghosts  are  attributed  to  the  dihedral-type 
reflections  from  the  wall-floor  junctions. 

2.5  CS  and  Multipath  Exploitation 

In  this  section,  we  consider  the  problem  of  multipath  in  view  of  the  requirements  of  fast  data 
acquisition  and  reduced  measurements.  Multipath  ghosts  may  cast  a  sparse  scene  as  a  populated 
scene,  and  at  minimum  will  render  the  scene  less  sparse,  degrading  the  performance  of  CS  based 
reconstruction.  A  CS  method  that  directly  incorporates  multipath  exploitation  into  sparse  signal 
reconstruction  for  imaging  of  stationary  scenes  with  a  stepped-frequency  monostatic  SAR  is 
presented.  Assuming  prior  knowledge  of  the  building  layout,  the  propagation  delays  corresponding 
to  different  multipath  returns  for  each  assumed  target  position  are  calculated,  and  the  multipath 
returns  associated  with  reflections  from  the  same  wall  are  grouped  together  and  represented  by  one 
measurement  matrix.  This  allows  CS  solutions  to  focus  the  returns  on  the  true  target  positions 
without  ghosting.  Although  not  considered  in  this  section,  it  is  noted  that  the  clutter  due  to  front 
wall  reverberations  can  be  mitigated  by  adapting  a  similar  multipath  formulation,  which  maps  back 
multiple  reflections  within  the  wall  after  separating  wall  and  target  returns  [100]. 

2.5.1  Multipath  Propagation  Model 

We  refer  to  the  signal  that  propagates  from  the  antenna  through  the  front  wall  to  the  target  and 
back  to  the  antenna  as  the  direct  target  return.  Multipath  propagation  corresponds  to  indirect 
paths,  which  involve  reflections  at  one  or  more  interior  walls  by  which  the  signal  may  reach  the 
target.  Multipath  can  also  occur  due  to  reflections  from  the  floor  and  ceiling  and  interactions  among 
different  targets.  In  considering  wall  reflections  and  assuming  diffuse  target  scattering,  there  are 
two  typical  cases  for  multipath.  In  the  first  case,  the  wave  traverses  a  path  that  consists  of  two  parts; 
one  part  is  the  propagation  path  to  the  target  and  back  to  the  receiver,  and  the  other  part  is  a  round 
trip  path  from  the  target  to  an  interior  wall.  As  the  signal  weakens  at  each  secondary  wall  reflection, 
this  case  can  usually  be  neglected.  Furthermore,  except  when  the  target  is  close  to  an  interior  wall, 
the  corresponding  propagation  delay  is  high,  and,  most  likely,  would  be  equivalent  to  the  direct-path 
delay  of  a  target  that  lies  outside  the  perimeter  of  the  room  being  imaged.  Thus,  if  necessary,  this 
type  of  multipath  can  be  gated  out.  The  second  case  is  a  bistatic  scattering  scenario,  where  the 
signal  propagation  on  transmit  and  receive  takes  place  along  different  paths.  This  is  the  dominant 
case  of  multipath  with  one  of  the  paths  being  the  direct  propagation,  to  or  from  the  target,  and  the 
other  involving  a  secondary  reflection  at  an  interior  wall. 

Other  higher-order  multipath  returns  are  possible  as  well.  Signals  reaching  the  target  can  undergo 
multiple  reflections  within  the  front  wall.  We  refer  to  such  signals  as  wall  ringing  multipath.  Also,  the 
reflection  at  the  interior  wall  can  occur  at  the  outer  wall-air  interface.  This  will  result,  however,  in 
additional  attenuation  and,  therefore,  can  be  neglected.  In  order  to  derive  the  multipath  signal 
model,  we  assume  perfect  knowledge  of  the  front  wall,  i.e.  location,  thickness,  and  dielectric 
constant,  as  well  as  the  location  of  the  interior  walls. 

2.5. 1.1  Interior  Wall  Multipath 

Consider  the  antenna-target  geometry  illustrated  in  Fig.  7(a),  where  the  front  wall  has  been  ignored 
for  simplicity.  The  pth  target  is  located  at  xp  =  ( xp,zp ),  and  the  interior  wall  is  parallel  to  the  z-axis 

and  located  at  x  =  xw.  Multipath  propagation  consists  of  the  forward  propagation  from  the  nth 
antenna  to  the  target  along  the  path  ;z?"and  the  return  from  the  target  via  a  reflection  at  the 
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interior  wall  along  the  path  p '.  Assuming  specular  reflection  at  the  wall  interface;  we  observe  from 
Fig.  7(a)  that  reflecting  the  return  path  about  the  interior  wall  yields  an  alternative  antenna-target 
geometry.  We  obtain  a  virtual  target  located  at  x'  =  ( 2xw  -xp9zp)9  and  the  delay  associated  with 

path  p'  is  the  same  as  that  of  the  path  p9  from  the  virtual  target  to  the  antenna.  This 
correspondence  simplifies  the  calculation  of  the  one-way  propagation  delay  r^pJ  associated  with 
path  p It  is  noted  that  this  principle  can  be  used  for  multipath  via  any  interior  wall. 

From  the  position  of  the  virtual  target  of  an  assumed  target  location;  we  can  calculate  the 
propagation  delay  along  path  P’  as  follows.  Under  the  assumption  of  free  space  propagation;  the 
delay  can  be  simply  calculated  as  the  Euclidean  distance  from  the  virtual  target  to  the  receiver 
divided  by  the  propagation  speed  of  the  wave.  In  the  TWRI  scenario;  however;  the  wave  has  to  pass 
through  the  front  wall  on  its  way  from  the  virtual  target  to  the  receiver.  As  the  front  wall  parameters 
are  assumed  to  be  known;  the  delay  can  be  readily  calculated  from  geometric  considerations  using 
Snell's  law  [28]. 


2.5. 1.2  Wall  Ringing  Multipath 


The  effect  of  wall  ringing  on  the  target  image  can  be  delineated  through  Fig.  7(b);  which  depicts  the 
wall  and  the  incident,  reflected,  and  refracted  waves.  The  distance  between  the  target  and  the  array 
element  in  crossrange  direction,  Ax,  can  be  expressed  as 

Ax  =  (Az-d)  tan  9air  +  d(  1  +  2iw)  tan  dwall  (2.29) 

where  A z  is  the  distance  between  target  and  array  element  in  downrange  direction,  and 
0air  and  0waU  are  the  angles  in  the  air  and  in  the  wall  medium,  respectively.  The  integer  iw  denotes 


the  number  of  internal  reflections  within  the  wall.  The  case  iw  =  0  describes  the  direct  path  as 


derived  in  [28].  From  Snell's  law, 


sin^>  =4~s 

Sln  Owall 


(2.30) 


Equations  (2.29)  and  (2.30)  form  a  nonlinear  system  of  equations  that  can  be  solved  numerically  for 
the  unknown  angles,  e.g.,  using  the  Newton  method.  Having  the  solution  for  the  incidence  and 
refraction  angles,  we  can  express  the  one-way  propagation  delay  associated  with  the  wall  ringing 
multipath  as  [101] 


T  = 


(A z-d)  ^  V ~sd(\  +  2iw ) 

CCOS0air  CCOS0wal, 


(2.31) 


2.5.2  Received  Signal  Model 

Having  described  the  two  principal  multipath  mechanisms  in  TWRI,  namely  the  interior  wall  and  wall 
ringing  types  of  multipath,  we  are  now  in  a  position  to  develop  a  multipath  model  for  the  received 
signal.  We  assume  that  the  front  wall  returns  have  been  suppressed  and  the  measured  data  contains 
only  the  target  returns.  The  case  with  the  wall  returns  present  in  the  measurements  is  discussed  in 
[100]. 

Each  path  p  from  the  transmitter  to  a  target  and  back  to  receiver  can  be  divided  into  two  parts,  p 9 
and  p"9  where  p"  denotes  the  partial  path  from  the  transmitter  to  the  scattering  target  and  pf  is 
the  return  path  back  to  the  receiver.  For  each  target-transceiver  combination,  there  exist  a  number 
of  partial  paths  due  to  the  interior  wall  and  wall  ringing  multipath  phenomena.  Let 
=0,1, ...,^-1,  and  P"2 ,  i2  =  09l9...9R2 -1,  denote  the  feasible  partial  paths.  Any 
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combination  of  pf  and  p?  results  in  a  round-trip  path  Pi9i  =  0,1,. — 1.  We  can  establish  a 

function  that  maps  the  index  /  of  the  round-trip  path  to  a  pair  of  indices  of  the  partial  paths, 
i  i-»  ( ,  i2  ).  Hence,  we  can  determine  the  maximum  number  R  <  R}R2  of  possible  paths  for  each 
target-transceiver  pair.  Note  that,  in  practice,^  «  RXR29  as  some  round-trip  paths  may  be  equal 
due  to  symmetry  while  some  others  could  be  strongly  attenuated,  and,  thereby  can  be  neglected. 
We  follow  the  convention  that  P0  refers  to  the  direct  round-trip  path. 


The  round-trip  delay  of  the  signal  along  path  pit  consisting  of  the  partial  parts  p?  and  p"9  can  be 
calculated  as 


T  n'  -  +  r 


.(0 
'  pyn 


'  P,n 


(232) 


We  also  associate  a  complex  amplitude  w{?  for  each  possible  path  corresponding  to  the  pth  target, 
with  the  direct  path,  which  is  typically  the  strongest  in  TWRI,  having  =  1 . 


Without  loss  of  generality,  we  assume  the  same  number  of  propagation  paths  for  each  target.  The 
unavailability  of  a  path  for  a  particular  target  is  reflected  by  a  corresponding  path  amplitude  of  zero. 
The  received  signal  at  the  nth  antenna  due  to  the  mth  frequency  can,  therefore,  be  expressed  as 

y(m, #?)  =  E  E y^ptrf  exp (-j(OmT^n )  (2.33) 

i=0  p= 0 

As  the  bistatic  radar  cross  section  (RCS)  of  a  target  could  be  different  from  its  monostatic  RCS,  the 
target  reflectivity  is  considered  to  be  dependent  on  the  propagation  path.  For  convenience,  the  path 
amplitude  (33)  can  be  absorbed  into  the  target  reflectivity  leading  to 

y(m,n)  =  £  E>J}  exp (2.34) 

i—0  p- 0 

Note  that  (2.34)  is  a  generalization  of  the  non-multipath  propagation  model  (2.2).  If  the  number  of 
propagation  paths  is  set  to  1,  then  the  two  models  are  equivalent. 

The  matrix-vector  form  for  the  received  signal  under  multipath  propagation  is  given  by 

y  =  ^(0)r(0)  +  ^(I)r(1)  + . . .  +  (2.35) 

where 

r(0  _  L(0  r(0  f 

r  -  [r00  ...  j 

[^/(0]J9  =exp(-y'fl)mr^),  m  =  smodM,n  =  [s/M]  (2.36) 

5  =  0,1 . MN-\,q  =  0,\ . ^Ar2-1 

The  term  q  =  0,1,... , NXNZ  - 1,  takes  the  value  <j^  if  the  pth  point  target  exists  at  the  pth 
pixel;  otherwise,  it  is  zero.  Finally,  the  reduced  measurement  vector  y  can  be  obtained  from  (2.35) 
as  y  =  Oy,  where  the  QXQ2  xMVmatrix  Ois  defined  in  (2.10). 

2.5.3  Sparse  Scene  Reconstruction  with  Multipath  Exploitation 

Within  the  CS  framework,  we  aim  at  undoing  the  ghosts,  i.e.,  inverting  the  multipath  measurement 
model  and  achieving  a  reconstruction,  wherein  only  the  true  targets  remain. 

In  practice,  any  prior  knowledge  about  the  exact  relationship  between  the  various  sub-images  r^of 
the  sparse  scene  is  either  limited  or  nonexistent.  However,  we  know  with  certainty  that  the  sub¬ 
images  r(0\rn\...,r(*_i)describe  the  same  underlying  scene.  That  is,  the  support  of  the  R  images 
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is  the  same,  or  at  least  approximately  the  same.  The  common  structure  property  of  the  sparse  scene 
suggests  the  application  of  a  group  sparse  reconstruction. 


All  unknown  vectors  in  (2.35)  can  be  stacked  to  form  a  tall  vector  of  length  NXN=R 

iT 


r  = 


.(0)7 


A  V 


.(*-i)7 


The  reduced  measurement  vector  y  can  then  be  expressed  as 

y  =  Br 

where  B  =  [oT(0)  Ot(1)  •••  W**-0]  has  dimensions  QiQ2  xNxNzR. 


(2.37) 

(2.38) 


We  proceed  to  reconstruct  the  images  r  from  y  under  measurement  model  (2.38).  It  has  been 
shown  that  a  group  sparse  reconstruction  can  be  obtained  by  a  mixed  lx  -/2norm  regularization 
[102-105].  Thus,  we  solve 


?  =  argminl||y-B?|g+a||?||2jl 


where  a  is  the  so-called  regularization  parameter  and 


2,1  * 

<7=0 


r(0)  (1) 
:q  ’  rq  ’ 


,(*-!) 


ri  -*s 

II 2  q= 0 


(2.39) 


(2.40) 


'  1=0 


is  the  mixed  l±  —  l2  norm.  As  defined  in  Eq.  (40),  the  mixed  lt  -  l2  norm  behaves  like  an  norm  on 
the  vector  ,  q  =  0,1, ---,NXNZ  -  1,  and  therefore,  induces  group  sparsity. 


In  other  words,  each  1^]71|  ,  and  equivalently  each  —  ^]r,  is 


'q  >'q  >  >'q 

encouraged  to  be  set  to  zero.  On  the  other  hand,  within  the  groups,  the  l2  norm  does  not  promote 
sparsity  [106].  The  convex  optimization  problem  (2.39)  can  be  solved  using  SparSA  [102],  YALL  group 
[103],  or  other  available  schemes  [105, 107]. 


Once  a  solution  r  is  obtained,  the  sub-images  can  be  noncoherently  combined  to  form  an  overall 
image  with  improved  signal-to-noise-and-clutter  ratio  (SCNR),  with  the  elements  of  the  composite 
image  rGS  defined  as 


IrjV®,. 

q  5  q  5 

5  q  J  | 

<7  =  0, 


NxN:-\. 


(2.41) 


2.5.4  Illustrative  Results 

An  experiment  was  conducted  in  a  semi-controlled  environment  at  the  Radar  Imaging  Lab,  Villanova 
University.  A  single  aluminum  pipe  (61  cm  long,  7.6  cm  diameter)  was  placed  upright  on  a  1.2  m  high 
foam  pedestal  at  3.67  m  downrange  and  0.31  m  crossrange,  as  shown  in  Fig.  8.  A  77-element 
uniform  linear  monostatic  array  with  an  inter-element  spacing  of  1.9  cm  was  used  for  imaging.  The 
origin  of  the  coordinate  system  is  chosen  to  be  at  the  center  of  the  array.  The  0.2  m  thick  concrete 
front  wall  was  located  parallel  to  the  array  at  2.44  m  downrange.  The  left  sidewall  was  at  a 
crossrange  of  -1.83  m,  whereas  the  back  wall  was  at  6.37  m  downrange  (See  Fig.  8).  Also,  there  was 
a  protruding  corner  on  the  right  at  3.4  m  crossrange  and  4.57  m  downrange.  A  stepped-frequency 
signal,  consisting  of  801  equally  spaced  frequency  steps  covering  the  1  to  3  GFIz  band  was  employed. 
The  left  and  right  side  walls  were  covered  with  RF  absorbing  material,  but  the  protruding  right 
corner  and  the  back  wall  were  left  uncovered. 
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We  consider  background  subtracted  data  to  focus  only  on  target  multipath.  Figure  9(a)  depicts  the 
backprojection  image  using  all  available  data.  Apparently,  only  the  multipath  ghosts  due  to  the  back 
wall,  and  the  protruding  corner  in  the  back  right  are  visible.  Hence,  we  only  consider  these  two 
multipath  propagation  cases  for  the  group  sparse  CS  scheme.  We  use  25%  of  the  array  elements  and 
50%  of  the  frequencies.  The  corresponding  CS  reconstruction  is  shown  in  Fig.  9(b).  The  multipath 
ghosts  have  been  clearly  suppressed. 

2.6  CS-Based  Change  Detection  for  Moving  Target  Localization 


In  this  section,  we  consider  sparsity-driven  CD  for  human  motion  indication  in  TWRI  applications.  CD 
can  be  used  in  lieu  of  Doppler  processing,  wherein  motion  detection  is  accomplished  by  subtraction 
of  data  frames  acquired  over  successive  probing  of  the  scene.  In  so  doing,  CD  mitigates  the  heavy 
clutter  that  is  caused  by  strong  reflections  from  exterior  and  interior  walls  and  also  removes 
stationary  objects  present  in  the  enclosed  structure,  thereby  rendering  a  densely  populated  scene 
sparse  [7,  9-10].  As  a  result,  it  becomes  possible  to  exploit  CS  techniques  for  achieving  reduction  in 
the  data  volume.  We  assume  a  multistatic  imaging  system  with  physical  transmit  and  receive 
apertures  and  a  wideband  transmit  pulse.  We  establish  an  appropriate  CD  model  for  translational 
motion  that  permits  formulation  of  linear  modeling  with  sensing  matrices,  so  as  to  apply  CS  for 
scene  reconstruction.  Other  types  of  human  motions  involving  sudden  short  movements  of  the 
limbs,  head,  and/or  torso  are  discussed  in  [70]. 


2.6.1  Signal  Model 


Consider  wideband  radar  operation  with  M  transmitters  and  N  receivers.  A  sequential  multiplexing 
of  the  transmitters  with  simultaneous  reception  at  multiple  receivers  is  assumed.  As  such,  a  signal 
model  can  be  developed  based  on  single  active  transmitters.  We  note  that  the  timing  interval  for 
each  data  frame  is  assumed  to  be  a  fraction  of  a  second  so  that  the  moving  target  appears  stationary 
during  each  data  collection  interval. 


Let  sT  (/)  be  the  wideband  baseband  signal  used  for  interrogating  the  scene.  For  the  case  of  a  single 
point  target  with  reflectivity  a p ,  located  at  \p  =  (xp,zp)  behind  a  wall,  the  pulse  emitted  by  the 
mth  transmitter  with  phase  center  at  =  {xm-zofr)  is  received  at  the  nth  receiver  with  phase 
center  at  xrn  =  {xm,-z0jy) in  the  form 

y mn  (0  =  (0  +  bmn  (0,  amn  (t)  =  apsT(t-  Tp  mn  )  exp(  -jcocrp  mn  )  (2.42) 

where  coc  is  the  carrier  frequency,  r p  mn  is  the  propagation  delay  for  the  signal  to  travel  between 

the  mth  transmitter,  the  target  at  and  the  nth  receiver,  and  bmn(t)  represents  the  contribution 

of  the  stationary  background  at  the  nth  receiver  with  the  mth  transmitter  active.  The  delay  rpmn 

consists  of  the  components  corresponding  to  traveling  distances  before,  through,  and  after  the  wall, 
similar  to  (3). 

In  its  simplest  form,  CD  is  achieved  by  coherent  subtraction  of  the  data  corresponding  to  two  data 
frames,  which  may  be  consecutive  or  separated  by  one  or  more  data  frames.  This  subtraction 
operation  is  performed  for  each  range  bin.  CD  results  in  the  set  of  difference  signals, 


58 


symn  (0  =  y{™]  (0  -  ^  (0  =  °t+X)  (0  -  «£  (0  (2-43) 

where  L  denotes  the  number  of  frames  between  the  two  time  acquisitions.  The  component  of  the 
radar  return  from  the  stationary  background  is  the  same  over  the  two  time  intervals,  and  is  thus 
removed  from  the  difference  signal.  Using  (42)  and  (43),  the  (m,  n)-th  difference  signal  can  be 
expressed  as, 

8ymn  (t)  —a  psT(t—T  )  exp(  -  ja>cx )-apsT(t  -r%„  )  exp( -jmcr(p]mn )  (2.44) 

where  r^mn  and  are  the  respective  two-way  propagation  delays  for  the  signal  to  travel 

between  the  mth  transmitter,  the  target,  and  the  nth  receiver,  during  the  first  and  the  second  data 
acquisitions,  respectively. 


2.6.2  Sparsity-Driven  Change  Detection  under  Translational  Motion 


Consider  the  difference  signal  in  (2.44)  for  the  case  where  the  target  is  undergoing  translational 
motion.  Two  nonconsecutive  data  frames  with  relatively  long  time  difference  are  used,  i .e.L  »  1 
[108].  In  this  case,  the  target  will  change  its  range  gate  position  during  the  time  elapsed  between  the 
two  data  acquisitions.  As  seen  from  (2.44),  the  moving  target  will  present  itself  as  two  targets,  one 
corresponding  to  the  target  position  during  the  first  time  interval  and  the  other  corresponding  to  the 
target  location  during  the  second  data  frame.  It  is  noted  that  the  imaged  target  at  the  reference 
position  corresponding  to  the  first  data  frame  cannot  be  suppressed  for  the  coherent  CD  approach. 

On  the  other  hand,  the  noncoherent  CD  approach  that  deals  with  differences  of  image  magnitudes 
corresponding  to  the  two  data  frames,  allows  suppression  of  the  reference  image  through  a  zero 
thresholding  operation  [23].  However,  as  the  noncoherent  approach  requires  the  scene 
reconstruction  to  be  performed  prior  to  CD,  it  is  not  a  feasible  option  for  sparsity-based  imaging, 
which  relies  on  coherent  CD  to  render  the  scene  sparse.  Therefore,  we  rewrite  (2.44)  as, 

2 

fymniO  =  l>°isT  (t  -  zj  mn ) exp(-j<X)cTj  mn )  (2.45) 
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(2.46) 


If  we  sample  the  difference  signal  8ymn(t)  at  times  {/*}£=(j  to  obtain  the  ATxl  vector  Aymn  and 
form  the  concatenated  NxN:x  1  scene  reflectivity  vector  r,  then  using  the  developed  signal 
model  in  (2.45),  we  obtain  the  linear  system  of  equations, 

^  mn='¥mnr  (2-47) 

The  qth  column  of  ^mri  consists  of  the  received  signal  corresponding  to  a  target  at  pixel  \q  and  the 
kth  element  of  the  qth  column  can  be  written  as  [70,  83] 


= 


sr({k  -Tq,mn)exp(-JacTq,mn) 


Jqymn  || 2 


,k q  =  0,1  ,...,nxn:  -1 


(2.48) 


where  Tg  mn  is  the  two-way  signal  traveling  time  from  the  mth  transmitter  to  the  qth  pixel  and  back 


v  q,mn 

to  the  nth  receiver.  Note  that  the  kth  element  of  the  vector  s^^is  sT(tk  ~rqmn ),  which  implies 
that  the  denominator  in  the  R.H.S.  of  (2.48)  is  the  energy  in  the  time  signal.  Therefore,  each  column 
of  has  unit  norm.  Further  note  that  if  there  is  a  target  at  the  qth  pixel,  the  value  of  the  qth 
element  of  r  should  be  oq\  otherwise,  it  is  zero. 
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The  CD  model  described  in  (2.47-2.48)  permits  the  scene  reconstruction  within  the  CS  framework. 
We  measure  a  J  («  K)  dimensional  vector  of  elements  randomly  chosen  from  Ay  mn.  The  new 
measurements  can  be  expressed  as 

=  9mA  ymn  =  Vmn'VmnV  (2-49) 

where  <pmn  is  a  J  x  K  measurement  matrix.  Several  types  of  measurement  matrices  have  been 

reported  in  the  literature  [[83],  [86],  [109]  and  the  references  therein].  To  name  a  few,  a 
measurement  matrix  whose  elements  are  drawn  from  a  Gaussian  distribution,  a  measurement 
matrix  having  random  ±1  entries  with  probability  of  0.5,  or  a  random  matrix  whose  entries  can  be 
constructed  by  randomly  selecting  rows  of  a  K  x  K  identity  matrix.  It  was  shown  in  [83]  that  the 
measurement  matrix  with  random  ±1  elements  requires  the  least  amount  of  compressive 
measurements  for  the  same  radar  imaging  performance,  and  permits  a  relatively  straight  forward 
data  acquisition  implementation.  Therefore,  we  choose  to  use  such  a  measurement  matrix  in  image 
reconstructions. 


Given  Aymn  for  m  =  0,1,. . .,M -1,  w  =  0,lv..,Af-l,  we  can  recover  r  by  solving  the  following 
equation, 

r  =  argminllrlL  subject  to  Ay  (2.50) 

r  11  M/i 

where 


^  ~  [Too  *  •  •  ^(m-ixam)  ]r  >  ®  “  diag((p00 ,  (p0l , . . . ,  <P(m-\xn-\)  ) 
Ay  =  [Ayoo  ^Yoi 


(2.51) 


Equations  (2.50,  2.51)  represent  one  strategy  that  can  be  adopted  for  sparsity-based  CD  approach, 
wherein  a  reduced  number  of  time  samples  are  chosen  randomly  for  all  the  transmitter-receiver 
pairs  constituting  the  array  apertures.  The  above  two  equations  can  also  be  extended  so  that  the 
reduction  in  data  measurements  includes  both  spatial  and  time  samples.  The  latter  strategy  is  not 
considered  in  this  section. 


2.6.3  Illustrative  Results 


A  through-the-wall  wideband  pulsed  radar  system  was  used  for  data  collection  in  the  Radar  Imaging 
Lab  at  Villanova  University.  The  system  uses  a  0.7ns  Gaussian  pulse  for  scene  interrogation.  The 
pulse  is  up-converted  to  3  GHz  for  transmission  and  down-converted  to  baseband  through  in-phase 
and  quadrature  demodulation  on  reception.  The  system  operational  bandwidth  from  1.5  -  4.5  GHz 
provides  a  range  resolution  of  5cm.  The  peak  transmit  power  is  25dBm.  Transmission  is  through  a 
single  horn  antenna,  which  is  mounted  on  a  tripod.  An  8-element  line  array  with  an  inter-element 
spacing  of  0.06m,  is  used  as  the  receiver  and  is  placed  to  the  right  of  the  transmit  antenna.  The 
center-to-center  separation  between  the  transmitter  and  the  leftmost  receive  antenna  is  0.28m,  as 
shown  in  Fig.  10.  A  3.65m  x  2,6m  wall  segment  was  constructed  utilizing  1cm  thick  cement  board  on 
a  2-by-4  wood  stud  frame.  The  transmit  antenna  and  the  receive  array  were  at  a  standoff  distance  of 
1.19m  from  the  wall.  The  system  refresh  rate  is  100Hz. 

In  the  experiment,  a  person  walked  away  from  the  wall  in  an  empty  room  (the  back  and  the  side 
walls  were  covered  with  RF  absorbing  material)  along  a  straight  line  path.  The  path  is  located  0.5m 
to  the  right  of  the  center  of  the  scene,  as  shown  in  Fig.  10.  The  data  collection  started  with  the 
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target  at  position  1  and  ended  after  the  target  reached  position  3,  with  the  target  pausing  at  each 
position  along  the  trajectory  for  a  second.  Consider  the  data  frames  corresponding  to  the  target  at 
position  2  and  position  3.  Each  frame  consists  of  20  pulses,  which  are  coherently  integrated  to 
improve  the  signal-to-noise  ratio.  The  imaging  region  (target  space)  is  chosen  to  be  3m  x  3m, 
centered  at  (0.5m,  4m),  and  divided  into  61  x  61  grid  points  in  crossrange  and  downrange,  resulting 
in  3721  unknowns.  The  space-time  response  of  the  target  space  consists  of  8  x  1536  space-time 
measurements.  For  sparsity-based  CD,  only  5%  of  the  1536  time  samples  are  randomly  selected  at 
each  of  the  8  receive  antenna  locations,  resulting  in  8  x  77  space-time  measured  data.  Figure  11 
depicts  the  corresponding  result.  We  observe  that,  as  the  human  changed  its  range  gate  position 
during  the  time  elapsed  between  the  two  acquisitions,  it  presents  itself  as  two  targets  in  the  image, 
and  is  correctly  localized  at  both  of  its  positions. 


2.7  CS  General  Formulation  for  Stationary  and  Moving  Targets 

As  seen  in  the  previous  sections,  the  presence  of  the  front  wall  renders  the  target  detection  problem 
very  difficult  and  challenging,  and  has  an  adverse  effect  on  the  scene  reconstruction  performance 
when  employing  CS.  Different  strategies  have  been  devised  for  suppression  of  the  wall  clutter  to 
enable  target  detection  behind  walls.  Change  detection  enables  detection  and  localization  of 
moving  targets.  Clutter  cancellation  filtering  provides  another  option  [87,  110].  However,  along  with 
the  wall  clutter,  both  of  these  methods  also  suppress  the  returns  from  the  stationary  targets  of 
interest  in  the  scene,  and  as  such,  allow  subsequent  application  of  CS  to  recover  only  the  moving 
targets.  Wall  clutter  mitigation  methods  can  be  applied  to  remove  the  wall  and  enable  joint 
detection  of  stationary  and  moving  targets.  However,  these  methods  assume  monostatic  operation 
with  the  array  located  parallel  to  the  front  wall,  and  exploit  the  strength  and  invariance  of  the  wall 
return  across  the  array  under  such  a  deployment  for  mitigating  the  wall  return.  As  such,  they  may 
not  perform  as  well  under  other  situations. 

For  multistatic  imaging  radar  systems  using  ultra-wideband  (UWB)  pulses,  an  alternate  option  is  to 
employ  time  gating,  in  lieu  of  the  aforementioned  clutter  cancellation  methods.  The  compact 
temporal  support  of  the  signal  renders  time  gating  a  viable  option  for  suppressing  the  wall  returns. 

This  enhances  the  SCR  and  maintains  the  sparsity  of  the  scene,  thereby  permitting  the  application  of 
CS  techniques  for  simultaneous  localization  of  stationary  and  moving  targets  with  few  observations 
[74]. 

2.7.1  Signal  Model 

Consider  the  scene  layout  depicted  in  Fig.  12.  Note  that  although  the  M-element  transmit  and  N- 
element  receive  arrays  are  assumed  to  be  parallel  to  the  front  wall  for  notational  simplicity,  this  is 
not  a  requirement.  Let  7J.be  the  pulse  repetition  interval.  Consider  a  coherent  processing  interval  of 
I  pulses  per  transmitter  and  a  single  point  target  moving  slowly  away  from  the  origin  with  constant 
horizontal  and  vertical  velocity  components  (vxpiv2p),  as  depicted  in  Fig.  12.  Let  the  target 

position  be  xp  =  ( xp,zp )  at  time  t  =  0.  Assume  that  the  timing  interval  for  sequencing  through 
the  transmitters  is  short  enough  so  that  the  target  appears  stationary  during  each  data  collection 
interval  of  length  ITr .  This  implies  that  the  target  position  corresponding  to  the  fth  pulse  is  given  by 

xp(i)  =  (xp  +vxpiITr,zp+v.piITr )  (2.52) 

The  baseband  target  return  measured  by  the  nth  receiver  corresponding  to  the  ith  pulse  emitted  by 
the  mth  transmitter  is  given  by  [74] 

y mm  (0  =  (XpsT  ( t  -  iITr  -  mTr  -  Tp  mn  (/)))  exp {-jO)cTp  mn  (/))  (2.53) 
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where  r  (f)  is  the  propagation  delay  for  the  rth  pulse  to  travel  from  the  mth  transmitter  to  the 

target  at  xp  (/),  and  back  to  the  nth  receiver.  In  the  presence  of  P  point  targets,  the  received  signal 

component  corresponding  to  the  targets  will  be  a  superposition  of  the  individual  target  returns  in 
(2.53)  with  p  —  0,1, ...,/*  —  1 .  Interactions  between  the  targets  and  multipath  returns  are  ignored  in 
this  model.  Note  that  any  stationary  targets  behind  the  wall  are  included  in  this  model  and  would 
correspond  to  the  motion  parameter  pair  {vxp,vzp)  =  (0,0).  Further  note  that  the  slowly  moving 
targets  are  assumed  to  remain  within  the  same  range  cell  over  the  coherent  processing  interval. 


On  the  other  hand,  as  the  wall  is  a  specular  reflector,  the  baseband  wall  return  received  at  the  nth 
receiver  corresponding  to  the  /th  pulse  emitted  by  the  mth  transmitter  can  be  expressed  as 

yZ‘!  (0  =  °wsT  (t  -  ilTr  -  mTr  -  rw  m„  )) exp( ~j(OcTwnm  )  +  (/)  (2.54) 

where  zwmn  is  the  propagation  delay  from  the  mth  transmitter  to  the  wall  and  back  to  the  nth 

receiver,  and  represents  the  wall  reverberations  of  decaying  amplitudes  resulting  from 

multiple  reflections  within  the  wall  (see  Fig.  13).  The  propagation  delay  rw  mn  is  given  by  [111] 


where 


'w,mn 


)2+*2 


off 


+ 


,)2+z. 


Off 


k  w,mn 


(2.55) 


(2.56) 


is  the  point  of  reflection  on  the  wall  corresponding  to  the  mth  transmitter  and  the  nth  receiver,  as 
shown  in  Fig.  13.  Note  that,  as  the  wall  is  stationary,  the  delay  rw  mn  does  not  vary  from  one  pulse  to 


the  next.  Therefore,  the  expression  in  (2.54)  assumes  the  same  value  for  /  =  0,1,...,/  —  1. 


Combining  (2.53)  and  (2.54),  the  total  baseband  signal  received  by  the  nth  receiver,  corresponding 
to  the  /th  pulse  with  the  mth  transmitter  active,  is  given  by 

y’rnm  (0 = yZn!  (0  +  iVl  (0  (2-57) 

P= 0 

By  gating  out  the  wall  return  in  the  time  domain,  we  gain  access  to  the  sparse  behind-the-wall  scene 
of  a  few  stationary  and  moving  targets  of  interest.  Therefore,  the  time-gated  received  signal 
contains  only  contributions  from  the  P  targets  behind  the  wall  as  well  as  any  residuals  of  the  wall  not 
removed  or  fully  mitigated  by  gating.  In  this  section,  we  assume  that  wall  clutter  is  effectively 
suppressed  by  gating.  Therefore,  using  (2.57),  we  obtain 

ymm(0=Y.yL(‘)  (2-ss) 

P= 0 


2.7.2  Linear  Model  Formulation  and  CS  reconstruction 

With  the  observed  scene  divided  into  Nx  xNz  pixels  in  crossrange  and  downrange,  consider  Nv 
and  Nv *  discrete  values  of  the  expected  horizontal  and  vertical  velocities,  respectively.  Therefore, 
an  image  with  NxxNz  pixels  in  crossrange  and  downrange  is  associated  with  each  considered 

horizontal  and  vertical  velocity  pair,  resulting  in  a  four-dimensional  target  space.  Note  that  the 
considered  velocities  contain  the  (0,  0)  velocity  pair  to  include  stationary  targets. 
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Sampling  the  received  signal  ymni (t)  at  times  {^}*=o  /  we  obtain  a  A:  x  1  vector  ymni.  For  the /th 
velocity  pair  we  vectorize  the  corresponding  crossrange  vs.  downrange  image  into  an 

NxNzx  1  scene  reflectivity  vector  i*(v*/>vz/)-The  vector  r(vx/,v2/)  is  a  weighted  indicator  vector 
defining  the  scene  reflectivity  corresponding  to  the  /th  considered  velocity  pair,  i.e.,  if  there  is  a 
target  at  the  spatial  grid  point  (x,  z)  with  motion  parameters  (vxj9vzj)9  then  the  value  of  the 
corresponding  element  of  r(vxl9vzl)  should  be  nonzero;  otherwise,  it  is  zero. 

Using  the  developed  signal  model  in  (2.53)  and  (2.58),  we  obtain  the  linear  system  of  equations, 

Ymm  ='¥mni(vxl,vzl)r(vxl,vzl),  l  =  0,\,...,NVxNV'-l  (2.59) 

where  the  matrix  (vx! , vz/ )  is  of  dimension  KxNxNz.  The  qth  column  of  xFm/J/(vJf/,vz/) 
consists  of  the  received  signal  corresponding  to  a  target  at  pixel  \q  with  motion  parameters 
(vxI,v:i),  and  the  /'th  element  of  the  qth  column  can  be  written  as 

['¥mni(vxi>vzi)h,q  =  sT(tk -ilTr -mTr -Tqmn(i))exp(-jo)cTqm„(i)),  q  =  0,1 . NXNZ  -1 

(2.60) 

where  rqmn{i)  is  the  two-way  signal  traveling  time,  corresponding  to  (vx/,vz/),  from  the  mth 
transmitter  to  the  qth  spatial  grid  point  and  back  to  the  nth  receiver  for  the  /th  pulse. 

Stacking  the  received  signal  samples  corresponding  to  /  pulses  from  all  MN  transmitting  and 
receiving  element  pairs,  we  obtain  the  MNIK  x  1  measurement  vector  y  as 

y  =  x¥(vxl,vzl)r(vxl,vzi),  l  =  0,\,...,(NVxNVi -1)  (2.61) 

where 

Vxt  9  V=l  )  =  [%00  (V^/  5  Vzl  )»  *  •  *  5  ^(M-lXAMXZ-l)  (V^/  5  Vzl  )]  *  (2.62) 

Finally,  forming  the  MNIKxN XN:NV*NV  matrix  4*  as 

'F  =  '¥(vX(nVxnVx-\)9 vz(jyVjcjyv^-j))]«/  (2.63) 

we  obtain  the  linear  matrix  equation 

y^'Fr  (2.64) 

with  r  being  the  concatenation  of  target  reflectivity  vectors  corresponding  to  every  possible 
considered  velocity  combination. 


The  model  described  in  (64)  permits  the  scene  reconstruction  within  the  CS  framework.  We  measure 
a  J  <  MNIK  dimensional  vector  of  elements  randomly  chosen  from  y.  The  reduced  set  of 
measurements  can  be  expressed  as 

y  =  OxFr  (2.65) 


where  O  is  a  J  x  MNIK  measurement  matrix.  For  measurement  reduction  simultaneously  along 
the  spatial,  slow  time,  and  fast  time  dimensions,  the  specific  structure  of  the  matrix  <D  is  given  by 


®  =  kron  (® , , I  )  •  kron  (®  2 ,  )  •  kron  (®  3 , 1  J{MN )  •  diag {<D 4(0) ,  O 4‘ , . . . , O 


(i) 


} 

(2.66) 


where  Iq  is  an  identity  matrix  with  the  subscript  indicating  its  dimensions,  and 
and  J2  denote  the  reduced  number  of  transmit  elements,  receive  elements,  pulses,  and  fast  time 
samples,  respectively,  with  the  total  number  of  reduced  measurements  J  =  J2-  The  matrix 

Oj  is  an  M{xM  matrix,  02  is  an  N{xN  matrix,  03  is  a  J2XI  matrix,  and  each  of  the  04 
matrices  is  a  matrix  for  determining  the  reduced  number  of  transmitting  elements,  receiving 
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elements,  pulses  and  fast  time  samples,  respectively.  Each  of  the  three  matrices  and 03 

consists  of  randomly  selected  rows  of  an  identity  matrix.  These  choices  of  reduced  matrix 
dimensions  amount  to  selection  of  subsets  of  existing  available  degrees  of  freedom  offered  by  the 
fully  deployed  imaging  system.  Any  other  matrix  structure  does  not  yield  to  any  hardware  simplicity 
or  saving  in  acquisition  time.  On  the  other  hand,  three  different  choices,  discussed  in  Section  VLB, 
are  available  for  compressive  acquisition  of  each  pulse  in  fast  time. 

Given  the  reduced  measurement  vector  y  in  (2.65),  we  can  recover  r  by  solving  the  following 
equation, 

r  =  argmin||r||  subject  to  (2.67) 

r  *1 

We  note  that  the  reconstructed  vector  can  be  rearranged  into  Nv  Nv  matrices  of  dimensions 
NxxNz  in  order  to  depict  the  estimated  target  reflectivity  for  different  vertical  and  horizontal 
velocity  combinations.  Note  that  i)  Stationary  targets  will  be  localized  for  the  (0,0)  velocity  pair,  and 
ii)  Two  targets  located  at  the  same  spatial  location  but  moving  with  different  velocities  will  be 
distinguished  and  their  corresponding  reflectivity  and  motion  parameters  will  be  estimated. 

2.7.3  Illustrative  Results 

A  real  data  collection  experiment  was  conducted  in  the  Radar  Imaging  Laboratory,  Villanova 
University.  The  system  and  signal  parameters  are  the  same  as  described  in  Section  VI.C.  The  origin  of 
the  coordinate  system  was  chosen  to  be  at  the  center  of  the  receive  array.  The  scene  behind  the 
wall  consisted  of  one  stationary  target  and  one  moving  target,  as  shown  in  Fig.  14.  A  metal  sphere  of 
0.3  m  diameter,  placed  on  a  1  m  high  Styrofoam  pedestal,  was  used  as  the  stationary  target.  The 
pedestal  was  located  1.25  m  behind  the  wall,  centered  at  (0.49  m,  2.45  m).  A  person  walked  towards 
the  front  wall  at  a  speed  of  0.7  m/s  approximately  along  a  straight  line  path,  which  is  located  0.2  m 
to  the  right  of  the  transmitter.  The  back  and  the  right  side  wall  in  the  region  behind  the  front  wall 
were  covered  with  RF  absorbing  material,  whereas  the  8  in  thick  concrete  side  wall  on  the  left  and 
the  floor  were  uncovered.  A  coherent  processing  interval  of  15  pulses  was  selected. 

The  image  region  is  chosen  to  be  4  m  x  6  m,  centered  at  (-0.31  m,  3  m),  and  divided  into  41  x  36 
pixels  in  crossrange  and  downrange.  As  the  human  moves  directly  towards  the  radar,  we  only 
consider  varying  vertical  velocity  from  -1.4  m/s  to  0  m/s,  with  a  step  size  of  0.7  m/s,  resulting  in 
three  velocity  pixels.  The  space-slow  time-fast  time  response  of  the  scene  consists  of  8  x  15  x  2872 
measurements.  First,  we  reconstruct  the  scene  without  time  gating  the  wall  response.  Only  33.3%  of 
the  15  pulses  and  13.9%  of  the  fast-time  samples  are  randomly  selected  for  each  of  the  8  receive 
elements,  resulting  in  8  x  5  x  400  space-slow  time-fast  time  measured  data.  This  is  equivalent  to 
4.6%  of  the  total  data  volume.  Figure  15  depicts  the  CS  based  result,  corresponding  to  the  three 
velocity  bins,  obtained  with  the  number  of  OMP  iterations  set  to  50.  We  observe  from  Figs.  15(a) 
and  15(b)  that  both  the  stationary  sphere  and  the  moving  person  cannot  be  localized.  The  reason 
behind  this  failure  is  two-fold:  1)  The  front  wall  is  a  strong  extended  target  and  as  such,  most  of  the 
degrees  of  freedom  of  the  reconstruction  process  are  used  up  for  the  wall,  and  2)  The  low  SCR,  due 
to  the  much  weaker  returns  from  the  moving  and  stationary  targets  compared  to  the  front  wall 
reflections,  causes  the  targets  to  be  not  reconstructed  with  the  residual  degrees  of  freedom  of  the 
OMP.  These  results  confirm  that  the  performance  of  the  sparse  reconstruction  scheme  is  hindered 
by  the  presence  of  the  front  wall. 

After  removal  of  the  front  wall  return  from  the  received  signals  through  time  gating,  the  space-slow 
time-fast  time  data  includes  8  x  15  x  2048  measurements.  For  CS,  we  used  all  eight  receivers, 
randomly  selected  5  pulses  (33.3%  of  15)  and  chose  400  Gaussian  random  measurements  (19.5%  of 
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2048)  in  fast  time,  which  amounts  to  using  6.5%  of  the  total  data  volume.  The  number  of  OMP 
iterations  was  set  to  4.  Figures  16(a),  16(b),  and  16(c)  are  the  respective  images  corresponding  to  the 
0  m/s,  -0.7  m/s,  and  -1.4  m/s  velocities.  It  is  apparent  that  with  the  wall  gated  out,  both  the 
stationary  and  moving  targets  have  been  correctly  localized  even  with  the  reduced  set  of 
measurements. 

2.8  Conclusion 

In  this  section,  we  presented  a  review  of  important  approaches  for  sparse  behind-the-wall  scene 
reconstruction  using  CS.  These  approaches  address  the  unique  challenges  associated  with  fast  and 
efficient  imaging  in  urban  operations.  First,  considering  stepped-frequency  SAR  operation,  we 
presented  a  linear  matrix  modeling  formulation,  which  enabled  application  of  sparsity  based 
reconstruction  of  a  scene  of  stationary  targets  using  a  significantly  reduced  data  volume.  Access  to 
background  scene  without  the  targets  of  interest  was  assumed  to  render  the  scene  sparse  upon 
coherent  subtraction.  Subsequent  sparse  reconstruction  using  a  much  reduced  data  volume  was 
shown  to  successfully  detect  and  accurately  localize  the  targets. 

Second,  assuming  no  prior  access  to  a  background  scene,  we  examined  the  performance  of  joint 
mitigation  of  the  wall  backscattering  and  sparse  scene  reconstruction  in  TWRI  applications.  We 
focused  on  subspace  projections  approach,  which  is  a  leading  method  for  combating  wall  clutter. 
Using  real  data  collected  with  a  stepped-frequency  radar,  we  demonstrated  that  the  subspace 
projection  method  maintains  proper  performance  when  acting  on  reduced  data  measurements. 

Third,  a  sparsity-based  approach  for  imaging  of  interior  building  structure  was  presented.  The 
technique  made  use  of  the  prior  information  about  building  construction  practices  of  interior  walls 
to  both  devise  an  appropriate  linear  model  and  design  a  sparsifying  dictionary  based  on  the 
expected  wall  alignment  relative  to  the  radar's  scan  direction.  The  scheme  was  shown  to  provide 
reliable  determination  of  building  layouts,  while  achieving  substantial  reduction  in  data  volume. 

Fourth,  we  described  a  group  sparse  reconstruction  method  to  exploit  the  rich  indoor  multipath 
environment  for  improved  target  detection  under  efficient  data  collection.  A  ray  tracing  approach 
was  used  to  derive  a  multipath  model,  considering  reflections  not  only  due  to  targets  interactions 
with  interior  walls,  but  also  the  multipath  propagation  resulting  from  ringing  within  the  front  wall. 
Using  stepped-frequency  radar  data,  it  was  shown  that  this  technique  successfully  reconstructed  the 
ground  truth  without  multipath  ghosts,  and  also  increased  the  SCR  at  the  true  target  locations. 

Fifth,  we  detected  and  localized  moving  humans  behind  walls  and  inside  enclosed  structures  using 
an  approach  that  combines  sparsity-driven  radar  imaging  and  change  detection.  Removal  of 
stationary  background  via  CD  resulted  in  a  sparse  scene  of  moving  targets,  whereby  CS  schemes 
could  exploit  full  benefits  of  sparsity-driven  imaging.  An  appropriate  CD  linear  model  was  developed 
that  allowed  scene  reconstruction  within  the  CS  framework.  Using  pulsed  radar  operation,  it  was 
demonstrated  that  a  sizable  reduction  in  the  data  volume  is  provided  by  CS  without  degradation  in 
system  performance. 

Finally,  we  presented  a  CS  based  technique  for  joint  localization  of  stationary  and  moving  targets  in 
TWRI  applications.  The  front  wall  returns  were  suppressed  through  time  gating,  which  was  made 
possible  by  the  short  temporal  support  characteristic  of  the  UWB  transmit  waveform.  The  SCR 
enhancement  as  a  result  of  time  gating  permitted  the  application  of  CS  techniques  for  scene 
reconstruction  with  few  observations.  We  established  an  appropriate  signal  model  that  enabled 
formulation  of  linear  modeling  with  sensing  matrices  for  reconstruction  of  the  downrange- 
crossrange-velocity  space.  Results  based  on  real  data  experiments  demonstrated  that  joint 
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localization  of  stationary  and  moving  targets  can  be  achieved  via  sparse  regularization  using  a 
reduced  set  of  measurements  without  any  degradation  in  system  performance. 
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FIGURE  CAPTIONS 

Figure  1.  Geometry  on  transmit  of  the  equivalent  two-dimensional  problem. 

Figure  2.  Imaging  results  after  background  subtraction,  (a)  Backprojection  image  using  full  data  (b) 
Backprojection  image  using  10%  data  volume,  (c)  CS  reconstructed  image  using  full  data,  (d)  CS 
reconstructed  image  using  10%  of  the  data. 

Figure  3.  CS  based  imaging  result  (a)  using  full  data  volume  without  background  subtraction,  (b) 
using  10%  data  volume  with  the  same  frequency  set  at  each  antenna. 

Figure  4.  (a)  Specular  reflections  produced  by  walls,  (b)  Indicator  function. 

Figure  5.  (a)  Crossrange  division  into  blocks  of  lx  pixels,  (b)  Sparsifying  dictionary  generation. 

Figure  6.  (a)  Scene  geometry,  (b)  reconstructed  image. 

Figure  7.  (a)  Multipath  propagation  via  reflection  at  an  interior  wall,  (b)  Wall  ringing  propagation 
with  iw= 1  internal  bounces. 

Figure  8.  Scene  Layout. 

Figure  9.  (a)  Backprojection  image  with  full  data  volume;  (b)  Group  sparse  reconstruction  with  25% 
of  the  antenna  elements  and  50%  of  the  frequencies. 

Figure  10.  Scene  Layout  for  the  target  undergoing  translational  motion. 

Figure  11.  Sparsity-based  CD  image  using  5%  of  the  data  volume. 

Figure  12.  Geometry  on  transmit  and  receive. 

Figure  13.  Wall  reverberations. 

Figure  14.  The  configuration  of  the  experiment. 

Figure  15.  Imaging  result  for  both  stationary  and  moving  targets  without  time  gating,  (a)  CS 
reconstructed  image  o(0,0),  (b)  CS  reconstructed  image  o(0,-0.7),  (c)  CS  reconstructed  image  o( 0,- 
1.4). 

Figure  16.  Imaging  result  for  both  stationary  and  moving  targets  after  time  gating:  (a)  CS 
reconstructed  image  o(0,0),  (b)  CS  reconstructed  image  o(0,-0.7),  (c)  CS  reconstructed  image  o(0,- 
1.4). 
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Figure  1.  Geometry  on  transmit  of  the  equivalent  two-dimensional  problem. 
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Figure  2.  Imaging  results  after  background  subtraction,  (a)  Backprojection  image  using  full  data  (b) 
Backprojection  image  using  10%  data  volume,  (c)  CS  reconstructed  image  using  full  data,  (d)  CS 
reconstructed  image  using  10%  of  the  data. 
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Figure  3.  CS  based  imaging  result  (a)  using  full  data  volume  without  background  subtraction,  (b) 
using  10%  data  volume  with  the  same  frequency  set  at  each  antenna. 


Figure  4.  (a)  Specular  reflections  produced  by  walls,  (b)  Indicator  function. 
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Figure  5.  (a)  Crossrange  division  into  blocks  of  /x  pixels,  (b)  Sparsifying  dictionary  generation. 


(b) 

Figure  6.  (a)  Scene  geometry,  (b)  reconstructed  image. 
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Figure  7.  (a)  Multipath  propagation  via  reflection  at  an  interior  wall,  (b)  Wall  ringing  propagation 
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Figure  8.  Scene  Layout. 


(a)  (b) 

Figure  9.  (a)  Backprojection  image  with  full  data  volume;  (b)  Group  sparse  reconstruction  with  25% 
of  the  antenna  elements  and  50%  of  the  frequencies. 
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Figure  10.  Scene  Layout  for  the  target  undergoing  translational  motion. 
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Figure  11.  Sparsity-based  CD  image  using  5%  of  the  data  volume. 
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Figure  12.  Geometry  on  transmit  and  receive. 
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Air 


Figure  13.  Wall  reverberations. 
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(b)  (c) 

Figure  15.  Imaging  result  for  both  stationary  and  moving  targets  without  time  gating,  (a)  CS 
reconstructed  image  cr(0,0),  (b)  CS  reconstructed  image  cr(0,-0.7),(c)  CS  reconstructed  image 
<t(0 -1.4). 
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Figure  16.  Imaging  result  for  both  stationary  and  moving  targets  after  time  gating:  (a)  CS 
reconstructed  image  cr(0,0);  (b)  CS  reconstructed  image  <r(0,-0.7),  (c)  CS  reconstructed  image 
£7(0, -1.4) 
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